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ABSTRACT 


This  report  describes  some  of  the  issues  involved  in  using  stepped- 
frequency  radar  pulse  trains  to  obtain  high  resolution  target  range 
profiles.  Radar  returns  from  stepped-f requency  pulse  trains  may  be  pro¬ 
cessed  either  coherently  or  noncoherently.  Coherent  processing  consists  of 
taking  the  Fourier  transform  of  the  N  coherently -detected  (complex)  pulse 
•  returns  from  each  range  cell,  where  N  is  the  number  of  pulses  in  the  train. 

The  output  of  the  transform  is  the  high  resolution  range  profile  of  the 
scatterers  within  that  range  cell  (i.e.,  within  the  span  of  a  single  pulse 
width).  Noncoherent  processing  consists  of  taking  the  Fourier  transform  of 
the  squared^magnitudes  of  the  pulse  returns  from  a  given  range  cell  and 
yields  the  auto-correlation  function  of  this  range'  profile.  In  either 
case,  the  resulting  range  resolution  (of  both  the  profile  and  itt’^auto- 
correlation  function)  is  determined  by  the  total  bandwidth  (i.e.,  the  fre¬ 
quency  spread)  of  the  pulse  train. 

The  stability  requirements  of  the  variable  frequency  source  used  to 
generate  ths  stepped  frequencies  of  the  pulse  train  are  determined  for  both 
the  coherent  and  noncoherent  processing  cases,  In  both  cases, r when  the 

•  l  —  •  .  , 

Fourier  transform  process  is  Implemented  using  a  DFT,  a  basic  requirement 
is  that  the  frequency  spacing  of  the  pulses  in  the  train  be  uniform. 

However,  the  frequency  stability  tolerance  required  for  coherent  processing, 
is  found  to  be  two  to  three  orders-of-magnltude  more  stringent  than  i a  the 
noncoherent  case.  The  effect  of  radar  platform* and/or  target  motion  on  the 
measured  range  profiles  and  range  profile  auto-correlation  functions  is  also 


ill 


determined.  In  general,  noncoherent  processing  is  found  to  be  much  more 
robust  than  coherent  processing  with  respect  to  the  effects  of  such  target 
and/or  radar  platform  motion. 

Finally  the  target-clutter  contrast  provided  by  the  resulting  range 
profiles  and  range  profile  auto-correlation  functions  is  determined.  The 
target-to-clutter  contrast  provided  by  the  range  profile  auto-correlation 
function  is  found  to  be  significantly  poorer  than  that  of  the  corresponding 
high-resolution  range  profile. 
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1.0  INTRODUCTION 


It  is  well  known  that  the  range  resolution  of  a  radar  is  inversely 
proportional  to  the  bandwidth  of  the  signal  waveform.  A  wide  bandwidth 
signal  waveform  may  be  realized  either  in  the  form  of  a  single  short  CW 
pulse,  a  single  swept-f requency  "chirp”  pulse,  or  a  train  of  pulses,  each 
of  which  is  transmitted  ?t  a  different, frequency.  However  in  most  prac- 
tical  radar  applications  involving  waveform  bandwidths  in  excess  of  20  MHz, 
only  the  latter  tvo  of  the  above  three  approaches  are  generally  used 
because  of  the  limited  energy  that  can  be  transmitted  in  a  very  short  CW 
pulse. 

The  frequency-diverse  pulse  train  may  be  considered  as  a  sequence  of 
samples  of  a  hypothetical,  single  swept-f requency  pulse,  the  time  duration 
and  frequency  spread  of  which  are  equal  to  those  of  the  pulse  train.  One 
advantage  of  using  a  stepped-f requency  pulse  train  as  a  high-resolution 
radar  waveform.  In  lieu  of  a  single  swept-f requency  pulse  having  the  same 
bandwidth,  is  that  with  the  pulse-train,  the  instantaneous  bandwidth  of  the 
radar  receiver  components  (i.e.,  mixers,  IF  amplifier,  etc.)  need  only  be 
as  great  as  that  of  a  single  pulse,  rather  than  as  large  as  the  frequency 
spread  over  the  entire  pulse  train.  In  addition,  coherent  frequency- 
synthesisers  for  generating  stepped-f requency  pulse  trains  are  generally 
more  robust  and  less  expensive  to  Implement  than  ultra-linear  wide-band 
chirp  waveforms  having  the  same  time-bandwidth  product.  This  issue  is  par¬ 
ticularly  germane  to  the  Implementation  of  wideband  waveforms  in  missile 
seeker  radars. 

Coherent  processing  of  the  returns  from  a  stepped-f requency  pulse  train 
is  similar  to  digital  pulse  compression  of  the  samples  of  a  hypothetical 
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continuous  chirp  waveform  because  of  the  discrete  nature  of  both  processes, 
and  the  associated  sampling  considerations.  As  in  the  case  of  normal  pulse 
compression,  it  yields  a  high  resolution  range  profile  of  the  scetterers 
contained  within  the  range  span  of  a  single  pulsewidth. 

When  a  high-bandwidth  waveform  is  implemented  in  the  form  of  a 
single  linearly  swept-f requency  "chirp“  pulse,  the  phase  at  every  point  in 
the  pulse  is  a  known,  continuous  quadratic  function  of  time.  A  reasonable 
questicn  is  whether  or  not  the  phases  of  the  individual  pulses  in  a 
frequency-diverse  pulse  train  also  need  to  be  related  in  such  a  known  way 
in  order  to  obtain  the  high  resolution  range  profile,  or  whether  it  suf¬ 
fices  for  the  radar  to, be  coherent  only  at  each  individual  frequency  of  the 
pulse  train,  independently?  The  Importance  of  this  issue  lies  in  the  fact 
that  if  the  latter  condition  of  independent  coherence  at  each  individual 
frequency  in  the  pulse  train  suffices,  then  i:  is  possible  to  generate  the 
N  frequencies  required  merely  by  using  N  different,  stable,  but  indepen¬ 
dent,  frequency  sources  such  as  crystal  oscillators.  Such  an  arrangement 
provides  the  require!  coherency  only  at  each  of  the  N  frequencies  indivi¬ 
dually,  since  t  '.e  phase  relationships  between  each  of  the  N  different 
stable,  frequency  sources  will  be  completely  arbitrary. 

1.1  Coherent  Processing  Overview 

It  is  shown  herein  that  such  an  arrangement  is  permissible  with  regard 
to  coherent  pro  :essing  of  the  pulse  returns..  Known  phase  relationships 
between  the  different  pulses  in  a  otepped-f requency  pulse  train  do  not 
appear  to  be  required  in  order  to  obtain  a  high  resolution  range  profile. 

If  this  were  not  true,  then  the  different  frequencies  of  the  pulses  in  the 


train  would  all  have  to  be  generated  from  a  common  coherent  source  such  as 
a  coherent  frequency-synthesizer.  However,  as  will  be  discussed  subsequently, 
it  is  often  mere  convenient,  for  other  reasons,  to  use  a  coherent  frequency 
synthesizer,  rather  than  N  independent  crystal  oscillators,  to  generate  the 
N  different  frequencies  of  the  pulse  train.  Therefore,  as  a  practical 
matter,  the  fact  that  N  independent  crystal  oscillators  could  be  used  is 
primarily  of  academic  interest. 

The  fact  that  high  range  resolution  can  be  obtained  by  coherently  pro¬ 
cessing  the  returns  from  a  frequency -stepped  pulse  train,  without  requiring  ' 
that  the  different  frequencies  in  the  pulse  train  be  generated  by  a 
coherent  transmitter,  appears  to  have  been  first  described  in  1S68,  in  the 
paper  "High  Range  Resolution  by  Means  of  Pulse-Pulse  Frequency  Shifting", 

Ref.  (1]«  in  that  paper  it  was  pointed  out  that  the  only  requirement  for 
generating  a  high-resolution  range  profile,  using  a  frequency-stepped  pulse 
train,  is  that  the  radar  have  a  coherent-on-receive  capability;  in  prin¬ 
ciple,  the  frequency-diverse  'ransaltted  pulses  can  be  generated  using  a 
non-coherent  source  such  as  a  variable  frequency  magnetron. 

However,  coherent  processing  of  the  returns  from  a  stepped-f requancy 
pulse  train  does  require  compensation  of  any  non-linear  differential  phase 
shifts  (l.e.  at  zero  range)  between  the  different  carrier  frequencies  that 
may  occur  in  the  front  end  of  the  radar  due  to  component  bandwidth  limita¬ 
tions.  In  addition,  any  target-doppler  induced  phase  shifts  in  the  pulse 
returns  must  also  be  compensated  for. 

In  the  case  of  coherent  processing  it  is  shown  that  the  DFT  Is  a 
matched  filter  for  the  returns  from  a  frequency-diverse  pulse  train,  under 
the  following  conditions. 
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Each  of  the  pulse  returns  is  coherently  detected  by  using  the 
transmitted  carrier  for  a  given  pulse  as  the  coherent  reference  for 
that  pulse.  Thus,  a  different  coherent  reference  is  used  for  each 
pulse  in  the  train. 

2.  The  spacing  beiveen  the  different  frequencies  at  which  the  pulses  in 
the  train  are  transmitted  is  uniform. 

3.  There  is  no  relative  radial  motion  between  the  target  and  the  radar 
during  the  transmission  of  the  pulse  train. 

i 

4.,  There  are  no  non-linear  differential  phase  shifts  through  the  radar 
system  between  the  different  carrier  frequencies. 

Performing  a  DPT  on  the  N  frequency  diverse  pulse  returns  from  each 
range  cell  is  similar  to  digital  pulse  compression. 

Because  of  the  fact  that  the  DFT  is  a  matched  filter  for  the  returns 
from  the  pulse  train  only  in  the  case  of  zero  radial  velocity,  the  effect 
of  non-zero  radial  velocity  is  also  investigated.  Two  effects  of 
uncompensated  non-zero  radial  velocity  on  the  resulting  high  resolution 
range  profile  are  noted.  The  first  is  the  well-known  rangs-doppler 
coupling  effect  that  causes  the  target's  range-profile  to  be  shifted  in 
range  by  an  amount  proportloual  to  the  product  of  the  relative  radial  velo¬ 
city  and  the  radar  carrier  frequency.  The  second  effect  is  an  attenuation 
and  dispersion  of  the  high  resolution  range  profile  produced  at  the  output 
of  the  OFT  due  to  the  filter  mismatch  between  the  DFT  and  the  doppler* 
shifted  pulse,  returns.  This  can  be  considerably  more  serious  than  the 
range  shifting  of  the  profile  because  the  resulting  dispersion  degrades  the 
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desired  resolution.  These  effects  are  shown  to  be  a  function  of  the  dimen¬ 
sionless  parameter,  V-,  that  is  equal  to  the  ratio  of  the  relative  radial 
motion,  NVT,  of  the  target  during  the  transmission  of  the  pulse  train, 
to  the  ideal  value  of  range  resolution  provided  by  the  pulse  train  band¬ 
width,  NAf.  Thus  P NVT/Ar,  where  V  is  the  radial  velocity,  T  is 
the  pulse  repetition  interval  and  Ar  is  the  range  resolution  corresponding 
to  the  pulse  train  bandwidth  NAf.  Serious  degradation  in  resolution 
of  the  range  profile  is  shown  to  occur  when  P>3.  These  effects  can,  of 
course,  be  eliminated  by  measuring  the  relative  radial  velocity  of  the 
target  and  applying  motion-compensation  phase  corrections  to  the  individual 
complex  pulse  returns.  The  most  practical' method  of  minimizing  P,  for  a 
given  radial  velocity  and  waveform  bandwidth,  short  of  applying  pulse-pulse 
muclon  compensation,  appears  to  be  to  reduce  the  duration  of  the  pulse  train 
by  increasing  the  prf. 

1.2  Noncoherent  Processing  Overview 

If  the  individual  frequency  sources  used  to  generate  the  N  frequencies 
of  the  pulse  train  are  not  all  phase  stable  and/or  the  radar  does  not  have 
a  coherent-on-receive  capability,  it  will  not  be  possible  to  generate  the 
high  resolution  range  profile  of  the  scatterers  encompassed  in  the  width  of 
a  single  pulse.  Such  might  be  the  case,  for  example,  if  the  N  frequencies 
were  generated  using  a  digitally-controlled  VCO  (voltage-controlled 
oscillator).  However,  in  this  event  it  is  still  possible  to  obtain  the 
auto-correlation  of  the  high  resolution  range  profile  by  noncoherently  pro¬ 
cessing  the  resulting  frequency-diverse  pulse  returns.  Knowledge  of  a 
radar  target's  range-profile  auto-correlation  function  is  sometimes  useful 
in  estimating  the  range  extent  of  that  target. 
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Whereas  coherent  processing  consists  of  computing  the  Fourier  trans¬ 


form  of  the  complex  pulse  returns  (phase  and  amplitude!  or  in-phase  and 
quadrature  voltages),  noncoherent  processing  consists  of  computing  the 
Fourier  transform  of  the  received  power  (or  voltage  amplitude  squared)  of 
the  frequency-diverse  pulse  returns.  The  term  "noncoherent  processing" 
derives  from  the  fact  that  coherent  detection  of  the  pulse  returns  is  not 
required  in  this  case  -  envelope  detection  suffices  to  measure  the  received 
power  of  the  returns.  However,  the  required  computational  effort  is  nearly 
the  same  in  both  cases. 

Although  only  the  «uto-cor relation  of  the  high-resolution  range  pro¬ 
file,  rather  than  the  range  profile  Itself,  can  be  obtained  in  this  way, 
noncoherent  processing  does  have  the  advantage  that  it  is  not  affected  by 
either  non-linear  differential  phase  shifts  over  the  frequency  span  of  the 
pulse  train,  or  by  target  motion.  Hence,  motion  compensation  is  not 
required,' and  the  process  is  therefore  somewhat  more  robust  than  coherent 
processing. 

1.3  Considerations  Affecting  Number  of  Frequencies  Required  (Number 
of  Pulses)  and  Frequency  Spacing 

When  the  Fourier  transform  is  implemented  in  terms  of  a  DFT,  an 
implicit  requirement  Is  that  the  frequency  spacing  between  the  M  frequen¬ 
cies  of  the  pulses  in  the  train  be  uniform.  This  requirement  applies 
regardless  of  whether  coherent  or  noncoherent  procwASlng  is  used,  and  is 
most  easily  met  by  using  a  coherent  frequency  synthesiser  as  the  variable 
frequency  source. 

Because  the  frequency-diverse  pulses  are  effectively  frequency  samples 
of  a  continuous,  frequency-swept  waveform,  the  limitations  of  the  tyqulst 
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sampling  tneorem  apply.  Since  the  sampling  is  performed  in  frequency 
rather  than  in  time,  the  frequency  spacing  between  the  different  pulses  in 
the  train  must  be  fine  enough  so  that  ambiguities  due  to  aliasing  of  the 
computed  high  resolution  range  profiles,  or  range  profile  auto-correlations, 

.  do  not  become  a  problem.  It  is  shown  that  the  maximum  frequency  spacing 
between  pulses  in  the  train  should  be  no  greater  than  one-half  the  single¬ 
pulse  bandwidth  (i.e.  one-half  the  reciprocal  of  the  pulse  width).  This 
limit  applies  to  both  coherent  and  noncoherent  processing  of  the  returns; 
however  for  somewhat  different  reasons  in  each  case.  In  the  former  case  it 
is  required  to  avoid  ambiguities  in  the  interpretation  of  the  resulting 
range  profile;  in  the  latter,  to  prevent  aliasing  of  the  corresponding  range 
profile  auto-correlation  function.  The  number  of  different  frequencies 
required  (or  equivalently,  the  number  of  pulses  in  the  tram)  is  determined 
by  the  ratio  of  the  pulse  train  bandwidth  (the  reciprocal  of  the  desired 
resolution)  divided  by  the  frequency  spacing  between  pulses.  Because  of 
the  constraint  on  frequency  spacing  described  above,  the  minimum  number  of 
frequencies  required  fs  equal  to  twice  the  time-bandwidth  product  obtained 
by  multiplying  the  pulse  train  bandwidth  by  the  widith  of  a  single  pulse. 

As  a  result,  frequency  diverse  pulse  trains  used  to  obtain  high  range  reso¬ 
lution  are  usually  at  least  32  pulses  long,  and  use  a  frequency  separation 
between  pulses  no  more  than  one  half, the  single  pulse  bandwidth. 

In  either  case,  the  combination  of  the  requirements  for  a  relatively 
large  number  of  frequencies  and  for  uniform  frequency  spacing  between  them 
is  best  met  by  use  of  i  digitally-controlled  coherent  frequency  synthe¬ 
sizer.  Such  a  device  ilso  provides  the  phase  stability  at  each  frequency 


that  is  required  for  coherent  processing.  Precise  uniformity  of  the 
spacing  between  adjacent  frequencies  is  often  difficult  to  achieve  even 
with  a  coherent  frequency  synthesizer  -  particularly  one  having  a  fast 
switching  time.  Consequently,  the  effect  of  non-uniformities  in  frequency 
spacing  on  processing  are  determined. 

In  theory ,  coherent  processing  to  obtain  the  range  profile  appears  to 
be  no  more  complicated  than  incoherent  processing  to  obtain  the  range  pro- 
file  auto-correlation.  However  in  practice,  the  successful  implementation 
of  coherent  processing  requires  both  accurate  calibration  and  compensation 
of  differential  phase  shifts  through  the  system  at  the  different  carrier 
frequencies,  extremely  tight  bounds  on  the  uniformity  of  the  frequency  spacing 
between  adjacent  pulses  In  the  train,  as  well  as  measurement  and  compen¬ 
sation  of  any  target  doppler  velocity,  although  the  latter  effect  can  often 
be  neglected  if  the  target  doppler  velocities  are  sufficiently  low.  On  the 
other  hand,  noncoherent  processing  is  not  affected  by  either  differential 
phase  shifts  or  relative  target  motion  and  the  bounds  on  the  allowable  fre¬ 
quency  errors  are  about  1000  time?  larger  than  for  coherent  processing. 
However,  it  yields  only  the  range  profile  auto-correlation  function  rather 
than  the  range  profile  itself.  Therefore,  there  is  a  trade-off  between 
these  two  techniques.  In  effect,  the  additional  costs  paid  for  obtaining 
the  range  profile  rather  than  its  autocorrelation  function,  are  accurate 
calibration  and  compensation  of  phase  shifts  that  result  from  either 
frequency-dependent  delays  in  the  radar  system  and/or  target  motion,  and 
tighter  bounds  on  the  allowable  non-unlformlty  In  the  spacing  of  the  fre¬ 
quencies  output  by  the  frequency  synthesizer. 


1.4.  Organization  of  Report 

A  mathematical  model  of  the  returns  received  from  a  distributed  target 
as  a  function  of  frequency,  and  the  processing  of  those  returns  to  obtain 
either  the  target's  high-resolution  range  profile,  or  the  auto-correlation 
of  that  profile,  is  presented  in  chapter  2.  An  expression  for  the  complex 
voltage  returned  from  a  distributed  target  as  a  function  of  transmitter 
frequency,  range,  pulse  shape  and  the  target's  range  scattering  profile  Is 
derived  in  Section  2.1.  Using  that  expression,  the  generation  of  the 
target's  high  resolution  range  profile  obtained  by  coherently  processing 
the  frequency  diverse  pulse  returns  from  the  range  cell  containing  that 
target,  is  described  in  Section  2.2.  Section  2.2.1  describes  how  the 
required  processing  of  the  coherently  detected  returns  is  Implemented  using 
a  DPT.  Section  2.2. 1.1  describes  the  problem  of  ambiguous  "wrap-around" 
of  the  resulting  high  resolution  range  profile  and  how  this  can  be  avoided 
by  limiting  the  value  of  the  inter-pulse  frequency,  spacing,  Af  to  a  value 
less  than  one-half  the  single-pulse  bandwidth,  and  by  properly  choosing  the 
timing  of  the  A/0  sampling  strobes.  The  effects  of  non-  uniform  frequency 
spacing  on  coherent  processing  are  described  in  Section  2.2.1. 2.  The  bound 
on  the  allowable  deviation  in  frequency  from  uniform  spacing  is  shown  to  be 
inversely  proportional  to  the  radar  prl.  Section  2.2.2  describes  coherent 
processing  of  the  stepped-f requency  pulse  returns  from  a  matched  filtering 
point  of  view,  as  well  as  the  differences  between  this  process  and  the  appli¬ 
cation  of  DFT's  in  the  compression  of  conventional  linear  FM  pulses.  The 
effects  of  non-zero  relative  radial  velocity  between  the  radar  and  the 
target  on  the  resulting  high  resolution  range  profile  are  described  in 


Section  2.2.3. 


Section  2.3  gives  a  similar  description  of  how  the  stepped-frequency 
pulse  returns  may  be  noncoherently  processed  to  obtain  the  auto-correlation 
function  of  this  high-resolution  range  profile.  It  is  shown  that  com¬ 
putation  of  the  resulting  range  profile  autocorrelation  function  is  not 
affected  either  by  target  velocity  or  by  non-linaar  differential  phase 
shifts  through  the  system,  and  that  the  bounds  on  the  allowable  non- 
uniformity  in  frequency  spacing  of  the  pulses  in  the  train  are  two  or  three 
orders  of  magnitude  greater  than  in  the  case  of  coherent  processing. 

Chapter  3  describes  the  effect  of  coherent  and  noncoherent  processing 
of  frequency-diverse  pulse  returns  on  the  amplitude  contrast  between  the 
target  and  clutter  responses  in  the  resulting  high  resolution  profiles. 

This  is  an  important  issue  because  the  detectability  of  targets  in  clutter 
is  based  upon  the  amplitude  contrast  between  the  return  from  a  target  and 
that  from  the  surrounding  clutter. 

In  the  unresolved  case,  the  principal  measure  of  this  contrast  is  the 
ratio  of  the  power  return  from  a  target  plus  the  clutter  within  a  single 
resolution  cell,  to  that  from  clutter  only  within  the  same  size  resolution 
cell.  .This  ratio  is  simply  (T+C)/C  and  is  usually  expressed  in  dB. 

In  the  resolved  case,  this  contract  is  defined  as  the  amplitude  ratio. 
In  dB,  of  the  target  response  to  that  of  the  mean  clutter  background  in  the 
high  resolution  profile  that  is  obtained  by  coherent  or  noncoherent  pro¬ 
cessing  of  the  frequency-diverse  pulse  returns.  The  difference  between  the 
target-clutter  contrast  in  the  high  resolution  profiles  obtained  after  pro¬ 
cessing,  and  the  original  unresolved  contrast  (T+C)/C,  is  defined  as  the 
contrast  enhancement  ratio,  E. 


Since  target  contrast  generally  depends  upon, resolution,  it  is  not 
surprising  to  find  that  E  is  primarily  a  function  of  the  resolution  improve¬ 
ment  factor,  I,  which  is  defined  as  the  ratio  of  the  original  single  pulse 
range  resolution  to  that  of  the  resulting  high  resolution  profile.  Note 
that  I  is  also  equal  to  the  time-bandwidth  product  of  the  pulse-train  band¬ 
width  and  the  width  of  a  single  pulse. 

Section  3*1  descsrlbes  the  effect  of  coherent  processing  on  target 
contrast.  Since  (T+C)/C  is  a  general  measure  of  target  contrast,  and  the 
clutter  return  C  is  linearly  proportional  to  the  area  of  the  resolution 
cell,  decreasing  the  size  of  the  resolution  cell  Increases  the  contrast 
between  targer  and  clutter  returns.  Consequently  for  coherent  processing, 
the  contrast  enhancement  factor  for  the  high-resolution  range  profile 
increases  in  proportion  to  the  Improvement  in  resolution  until  tbe  resolu¬ 
tion  becomes  fine  enough  to  start  separating  (i.e. ,  resolving)  the  indivi¬ 
dual  scatterer3  of  the  target.  At  that  point,  E  continues  to  increase  with 
I  but  at  a  somewhat  lower  rate.  Asymptotically,  E  -  I/M,  wher*'  M  is  the 
number  of  resolved  target  scatterers  in  the  high  resolution  profile. 

The  effect  of  noncoherent  processing  on  the  contrast  between  the 
target  responses  and  the  clutter  background  of  the  resulting  range  profile 
auto-correlation  function  is  described  in  Section  3.2.  It  Is  seen1  that, 
contrary  to  coherent  processing,  the  target-to-clutter  contrast  in  the 
resulting  range  profile  auto-correlation  function  is  often  less  than  the 
original  unresolved  target  contrast  (T+C)/C.  Thus,  the  contrast  enhan¬ 
cement  factor  for  noncoherent  processing  is  often  less  chan  unity  (or  nega¬ 
tive  in  dB).  Nevertheless,  E  does  Increase  with  increasing  resolution,  but 
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at  a  ruch  slower  rate  than  for  coherent  processing,  and  in  a  fairly  compli¬ 
cated  way.  In  this  case,  in  addition  to  being  a  function  of  I  and  M,  the 
contrast  enhancement  ratio  E  is  also  a  function  of  the  unresolved  t&rget- 
to-clutter  rnfic,  T/C. 

Finally,  some  computer  simulations  of  actual  high-resolution  range 
profiles  and  range  profile  autocorrelation  functions  that  are  obtained 
using  stepped-f requency  pulse  trains,  are  presented  in  Chapter  4.  The  pur¬ 
pose  of  these  simulations  is  to  illustrate  some  of  the  more  important 
attributes  of  these  profiles  that  were  theoretically  predicted  in  Chapters 
2  and  3.  Specifically,  these  simulation  results  illustrate  the  effect  of 
resolution  improvement  on  target  contrast  enhancement  and  the  effect  of 
uncompensated  radial  velocity  on  the  resulting  high  resolution  range- 
profiles  and  range-profile  autocorrelation  functions. 

Simulated  high  resolution  range-profiles  and  range-profile- 
autocorrelation  functions  are  generated  for  the  following  cases: 

1.  Clutter  only 

2.  A  stationary  two-point-scatterer  target  in  clutter  at  various 
target-to-clutter  ratios 

3.  A  stationary  three-point-scatterer  target  in  clutter 

4.  A  stationary  two-point-scatterer  target  in  clutter,  observed  by  a 
moving  radar 

5.  A  two-point-scatterer  target  moving  through  clutter  and  observed 
by  a  stationary  radar. 


Since  Che  clutter  model  is  random,  the  simulation  is  Monte-Carlo  and 
the  resulting  high  resolution  profiles  are  themselves  random  functions. 

Hence  the  simulation  results  are  presented  in  terms  of  the  mean  and  the 
mean-plus-three-times  the  standard  deviation  of  the  high  resolution  range 
profile  and  range  profile  autocorrelation  functions.  In  all  cases,  the 
results  of  these  simulations  confirm  the  profile  attributes  that  are 
theoretically  predicted  in  Chapters  2  and  3. 

The  simulation  results  clearly  illustrate  the  superior  target  contrast 
enhancement  provided  by  coherent  processing  relative  to  that  provided  by 
noncoherent  processing.  Specifically,  in  the  case  of  coherent  processing, 
the  target  contrast  with  respect  to  the  clutter  in  the  high  resolution 
range  profile  is  seen  to  increase  in  direct  proportion  to  tne  improvement 
in  resolution  (i.e,  pulse  train  bandwidth).  However,  in  the  case  of  nonco¬ 
herent  processing,  the  target-clutter  contrast  in  the  resulting  range  pro¬ 
file  auto-correlation  function  is  not  only  significantly  less  than  in  the 
case  coherent  processing,  but  also  increases  only  in  proportion  to  the 
square-root  of  the  resolution  improvement  ratio. 

Despite  the  lower  contrast  enhancement  attained  using  non-coherent 
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processing,  the  moving  target/moving  radar  simulations  illustrate  the 
relative  robustness  of, the  resulting  range-profile  auto-eorrelation  function 
with  respect  to  uncompensated  non-sero  radial  velocities  or  other  sources 
of  phase  error. 


2.0  ANALYSIS 


2.1  Overview  of  Analysis 

In  this  chapter,  a  mathematical  model  of  the  radar  pulse  returns 
received  from  a  distributed  target  is  derived  as  a  function  of  radar  fre¬ 
quency,  sampling  time,  pulse  shape,  and  the  target's  complex  range  scattering 
profile.  The  distributed  target  is  assumed  to  be  unresolved  -  i.e.  the 
range  extent  of  target's  range  scattering  profile  is  assumed  to  be  less 
than  the  radar  pulse  width.  This  model  is  then  used  tc  develop  the  pro¬ 
cessing  required  to  obtain  the  target's  high  resolution  range-profile  and 
range-profile  auto-correlatioa, function  from  the  radar  returns  received  in 
response  to  transmission  of  a  frequency-agile  -pulse  train. 

The  approach  taken  in  this  analysis  is  perhaps  unconventional  in  that 
it  emphasizes  the  basis  physics  of  the  processes  involved  rather  than 
abstract  mathematical  concepts.  The  approach  is  as  follows.  An  expression 
is  derived  for  the  complex  return  from  a  radar  target  as  a  function  of 
range,  carrier  frequency  and  A/D  sampling  time.  When  the  transmitted  pulse 
carrier  Is  used  as  the  coherent  reference  signal,  the  phase  of  the  detected 
radar  return  from  the  target  is  equal  to  the  product  of  that  carrier  fre¬ 
quency  and  the  round  trip  delay  time  corresponding  to  the  target  range. ’ 

The  A/D  sampling  time  at  which  the  complex  target  return  is  measured  prin¬ 
cipally  affects  only  the  amount  of  attenuation  of  the  measured  target 
return  by  the  shape  of  the  radar  pulse.  However,  as  will  be  shown  in 
Section  2.2. 1.1,  in  the  case  of  coherent  processing,  the  timing  of  the  A/D 
sampling  strobes  must  be  synchronized  to  the  fundamental  freque.  ..y  source 
(l.e.,  synthesizer)  used  to  generate  the  stepped  carrier  frequencies 


of  che  pulse  Crain,  in  order  Co  avoid  aliasing  of  Che  high  resoluCioo  range 
profiles. 

The  resulting  expression  for  Che  complex  return  from  a  distributed 
target  as  a  function  of  frequency  is  then  shown  to  be  the  Fourier  transform 
of  the  target  range  profile,  weighted  by  the  pulse  shape.  Consequently, 
the  discrete  Fourier  transform  of  the  range  profile  is  obtained  by 
assuring  the  complex  return  from  the  target  at  a  number  of  different  fre¬ 
quencies.  Therefore  the  desired  range  profile  can,  in  theory,  be 
reconstructed  by  taking  the  Inverse  transform  of  the  set  of:  complex  returns 
measured  at  these  different  frequencies.  When  using  a  frequency  agile  pulse 
train,  these  measurements  of  the  complex  target'  returns  are  made  at  a 
finite  number  of  discrete  frequencies  equal  to  the  number  of  pulses  in  the 
train.  Therefore,  as  a  practical  matter,  the  required  inverse  transform  is 
implemented  using  the  DFT.  Indeed,  it  Is  shown  that  in  the  case  of  zero 
target  radial  velocity,  DFT  is  the  matched  filter  for  this  pulse  train. 

2.1.1  Derivation  of  Complex  Returns  from  an  Arbitrary 
Distributed  Target 

Let  the  output  of  a  variable  frequency  source  in  the  exciter  of  a 
radar,  for  frequency  f^  be  given  by: 

u(t)  *  e-^2***6  +  Sk  *  +k(t»  (2-1) 

where  3^  is  an  arbitrary  constant  phese  angle  and  ^1.:)  represents  the 
phase  drift  of  the  oscillator  with  respect  to  the  constant  frequency  ffc,  at. 
a  function  of  time .  The  time  origin  corresponds  to  zero  range  -  the  time 
of  transmission  of  a  pulse. 


c 


Let  the  same  frequency  source  also  serve  as  the  receiver  local  oscillator 
signal  for  the  return  signal  received  frcm  a  target  in  response  to  the  trans¬ 
mission  of  u(t).  Hence,  this  L.O.  reference  signal  may  also  be  written  as: 


r(t)  .  e J ( 2u ffct  +  9 ^  +  i|>ic(t)) 


(2-2) 


Note  that  — —  represents  the  frequency-drift  of  the  local  oscillator  over 

the  tine  interval,  t,  between  the  time  of  transmission  of  the  pulse  and  the 

reception  of  the  return  signal  from  the  target.  By  definition;  -^(0)  “  0 
<fyk(0) 

and  — —  -  0.  In  the  case  of  a  stable  frequency  source,  such  as  a 

crystal  oscillator  i|;k(t)  *  0  for  all  t.  The  signal  reflected  from  a  point 

target  is  the  original  transmitted  signal,  aft),  delayed  by  the  ro  .nd  trip 

time,  x,  and  modulated  by  the  targets  complex  reflection  coefficient. 

y(t)  *  At  u(trx)eJ*c  (2“3a> 

where  At  and  are  the  amplitude  and  phase,  respectively,  of  the  target's 

complex  reflection  coefficient  and  where  x  is  the  round-trip  time  to  the 

point  target,  x  is  related  to  the  range  r  by: 


-  £ 

c 


(2-3b) 


After  mixing  with  the  L.O.  and  reference  signals  the  received  signal  from 
the  target  becomes: 

v(t)  -  AteJ*c.u(t-x)r*(t)  -  ♦  ♦k(t))  (2-4) 

Note  that  the  arbitrary  phase  angle  6g  of  the  transmitted  signal  does  not 
affect  the  output  of  the  mixer  since  the  same  frequency  source  is  used  to 
generate  both  the  transmitted  signal  and  local  oscillator  reference. 
However,  the  signal  v(t)  in  eqn  (2-4)  Is  affected  by  the  phase  drift  ^(t) 
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of  Che  exciter  away  froa  Che  fixed  frequency  f^,  as  indicated  by  Che 
difference  term:  ^(t)  -  ^j-(t-x). 

The  above  analysis  represents  the  complex  return  signal  received  froa 
a  point  target,  at  range  r,  at  the  output  of  the  receiver  mixer.  When  the 
transmitted  signal  is  a  finite  duration  pulse,  the  effect  of  the  pulse 
shape  on  the  received  signal  is  illustrated  in  Pig.  2-1.  The  amplitude  of 
Che  received  return  is  attenuated  by  the  pulse  shape  by  an  amount  S(ta-ix0)t 
where  S  is  the  pulse  shape  function,  t8  is  the  A/D  saapling  time,  and  Xq  is 
the  range  delay  to  the  point  target. 

Nov  consider  the  return  froa  a  distributed  target,  in  response  to  a 
finite  duration  pulse,  as  illustrated  in  Fig.  2-2.  The  dimensions  of  both 
the  target  range  profile  ahd  the  pulse  shape  are  given  in  units  of  tiae. 

In  the  sequel ,  the  designation  E(t)  will  be  used  to  represent  the  complex 
received  signal  before  mixing  with  the  L.O.  and  coherent  reference  signals, 
while  V(t),  will  represent  Che  complex  received  signal  after  mixing  (l.e. 
the  output  of  the  coherent  detector).  As  can  be  seen  froa  Fig.  2-2,  a 
sauple  of  the  received  signal  (before  mixing  with  the  L.O.  reference)  at 
sample  time  t,  due  to  a  differential  scattering  element  p(x)dx  located  at 
range  delay  (l.e.  round  trip  time),  x,  is  thus; 

dE(f  ,t)  -  p(x)s(t-x)e^2"f,t(t"*)  +  ♦  0k)dx  (2-5) 

In  eqn  (2-5),  p(x)  represents  the  complex  scattering  profile  of  the  target, 
s(x)  represents  the  real  pulse  shape,  and  a^2wf|c(t-x)  ♦ir.tt-x)  +  9k) 
represents  the  phase  of  the  return  froa  the  scattering  element  p(x)dx.  As 
before,  0^  is  an  arbitrary  phase  angle.  Consequently,  the  received  voltage 
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A 

|  136795  NJ 

POINT  TARGET 

AT  RANGE  Xq 

■ 

*0 

RANGE  DELAY  (Time),  * 


RANGE  DELAY  TIME,  ? 
ts  =  A/D  SAMPLING  INSTANT 


MAGNITUDE  GF  TARGET  RESPONSE  AT  SAMPLING  TIME  ts 
|E(f,«s>|  =  A  5* 

Fig.  2-1  Single  pulse  return  from  «  polnt-scatterer  target  of  amplitude  A 
and  located  at  range-delay  Xo>  measured  by  the  radar  receiver  at  sampling 
time,  tM* 
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RANGE  MAGNITUDE  PROFILE  OF  DISTRIBUTED  TARGET 


|i3«;9«  n) 


RANGE  DELAY  TIME,  t 


!,  =  A/D  SAMPLING  TIME 

TOTAL  RETURN  FROM  DISTRIBUTED  TARGET  p(x) 
AT  SAMPLING  TIME  t. 


E<fk.g  =  J>x)  S  *  H‘*dx 

-ac 

E(fK.gs  vA(x))eJ«)S(lt•x1^•K^B,k(,*■‘,,^»^,»'*•,  *  HklAx 


Fig.  2-2  Formulation  of  the  complex  return  from  a  distributed  target  as 
the  convolution  of  the  target's  complex  range  scattering  profile  and  the 
radar  pulse  shape. 
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at  sample  time,  t,  due  Co  a  distributed  target,  p(x),  is  the  coherent  sum 
of  all  of  the  different  scattering  elements  in  the  target,  and  is  repre¬ 
sented  by  the  following  convolution  integral: 

E(fk,t)  -  /  p(x)s(t-x)  ^ 2irf|c(t— x>  +  ^(t-x)  +  3k>dx  (2-6) 

Hence,  for  a  given  scattering  profile,  p(x),  the  received  voltage  E(f,t) 
depends  not  only  upon  the  pulse  shape  s(*)  and  on  the  time  of  the  sampling 

Instant,  t,  but  also  upon  the  frequency  of  the  transmitted  pulse,  fk.  This 

frequency  is  the  frequency  of  the  transmitted  signal  at  the  time  of  pulse 
transmssion;  t-0,  and  in  accordance  with  eqn  (2-1)  is  equal  to 

Now  suppose  that  E(fk,x)  is  to  be  determined  for  a  number  of  discrete 
frequencies,  fk.  Tf  these  frequencies  are  generated  Independently  of  each 
other,  then  the  arbitrary  initial  phase  angle,  ek,  will  also  be  unique  for 
each  frequency.  Thus  the  complex  pulse  return  at  frequency  fk  is 

E(fk,t>  -  eJ<2,rfkt  +  p(x)s(t-x)e*j(2irfkx  “  *k(t-*»dx  (2-6a) 

Eqn  ( 2-6a)  can  be  further  simplified  because  ^i(t-x)  ■  0.  An  explanation  of 
why  this  is  so  follows:  let  the  nominal  width  of  the  pulse,  whose  shape  is 
defined  by  s(x)  be  t.  Then  for  all  practical  purposes,  s(t-x)  *  0  for 
J  t-x  |  >  2t.  Furthermore,  since  by  definition  t|»(0)  •  0  and  i|>(0)  -  0,  it 

also  follows  that  iji(t-x)  -  0  when  j  t-x  |  <  2r»  Hence,  the  above  may  be 

approximated,  as  follows: 


fk  + 


<tyk 

dt 


t-0 


“  fk; 


since 


#k 

dt 


t-0 


-  0,  by  definition. 
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t+2i 

E(fk,t)  -  ej<2nfkt  +  6k)  jT  p(x)s(t-x)e“J<2ltfkX)dx  (2-6b) 

C-2t 

or  equivalently  fcince  s(t-x)  ”  0  for  j  t-xj  >2t: 

E(£k,t>  -  .•><2*fl‘t  +  6k)  f  p(x)s(t-x)e”^2ir^kX^dx  (2-6c) 

Then  after  mixing  this  signal  with  the  L.O.  reference  signal  of  eqn  (2-2)  to 
base-band  complex  video,  the  complex  voltage  at  the  output  of  the  detector ' 
becomes: 

+• 

V(fk,t)  ■  e'^k(t>  J  p(x)s(t-x)^2irfkX  dx  (2-7) 

Note  that  the  mixing  operation  removes  the  frequency  fk  and  the 
corresponding  initial  phase  angle  ek,  but  Introduces  the  phase  change  <|>k(t) 
corresponding  to  the  frequency  drift  of  the  L*  0.  during  the  round-trip 
signal  delay,  t. 

Now  define:  Q(fk,t)  =  J  p(x)s(t-x)e  ^2xfk*dx  (2-8) 

Thus  eqn  (2-7)  may  be  written  as:  V(fk,t)  •  e^^k^C^Q(fk,t)  (2-?8a) 

Eqn  (2-8)  may  be  Interpreted  as  the  definition  of  the  following  Fourier 
transform  pair, 

Q(f,f)  < — >  q(x,t)  =  p(x)s(t— x)  (2-9) 

where  the  time-domain  variable  is  x,  and  t  is  a  parameter  representing  the 
sampling  time  instant.  The  function  q(x,t)  represents  the  scattering  pro¬ 
file  of  the  target,  weighted  by  the  pulse  shape,  s(x). 


p(x) s(t-x)  -  f  ejl('(f,t>V(f,t)ej2fffxdf  (2-10a) 

where  Che  phase  drift  funcClon  ^(t)  has  been  wriCCen  as  ip(f,t)  since  Che 
phase  drift  will  generally  also  depend  upon  the  transmit tea  frequency,  f. 
Recall  from  the  definition  in  eqn  (2-7)  that  V(f,t)  is  the  complex  output 
voltage  of  the  coherent  detector  in  response  to  a  pulse  transmitted  at  fre¬ 
quency  f.  The  quantity  ip(f,t)  represents  the  combination  of  the  phase 
drift  of  the  master  oscillator  between  the  time  at  which  the  pulse  was 
transmitted  at  frequency  f,  and  the  time,  t,  at  which  a  return  from  the 
target  was  received  in  response  to  that  pulse,  and  any'  non-linear  differen¬ 
tial  phase  delay  through  the  transmitter  and  microwave  components  as  a 
function  of  carrier  frequency.  The  latter  differential  phase  delays  are 
only  a  function  of  frequency  and  are  independent  of  time,  and  can  be  repre¬ 
sented  by  ^»(f  ,o)  -  \|»0(f )  a  o.  For  the  present,  the  range  to  the  target  is 
assumed  to  remain  constant.  The  effect  of  relative  target  motion  on 
coherent  processing  will  be  examined  subsequently  In  Section  2.2.3. 

Actually,  when  using  a  frequency-diverse  pulse  train,  the  measurements 
of  V(f,t)  are  not  made  over  a  continuum  of  frequencies,  but  rather  at  a  set 
of  discrete  frequencies,  ffc.  Hence  the  continuous  transform  in  (2-10a)  is 
actually  approximated  In  discrete  form,  follows: 


N-i 

p(x)a(t-x)  -  q(x,t)  »  l  ejlt,k(t)  V(fk,t)eJ2ltfltX  Af  (2-11) 

k-0 

where  ^(t)  =  \|i(fk,t) 

The  summation  process  defined  by  eqn.  (2-11)  effectively  constitutes  a  special 
case  of  digital  pulse  compression.  V(f]ltt)  Is  the  complex  voltage  measured 
at  the  output  of  the  coherent  detector  at  each  frequency  k-0  . ...N-l. 
These  measurements  must  be  corrected  for  the  phase  drift,  ^(t),  of  the 
local  oscillator  during  the  round-trip  time,  t,  to  the  target,  relative  to 
the  phase  of  a  hypothetical  constant  frequency  signal  of  frequency  f^  at 
which  the  given  pulse  was  transmitted.  If  the  source  of  each  of  the  N  fre¬ 
quencies,  ffc,  k  -  1....N  are  stable,  and  there  are  no  non-linear  differen¬ 
tial  phase  delays  as  a  function  of  frequency,  then  \f*^(t)  -  0  for  al.l  k  and 
t  and  no  corrections  to  the  measured  complex  return  V(f]c,t)  are  required. 
However  If  the  frequency  sources  are  not  stable,  or  the  system  does  not 
have  a  linear  phase  response  as  a  function  of  frequency,  then  coherent  pro¬ 
cessing  (eqn.  11)  to  obtain  the  range  profile  q(x,t)  -  p(x)«(t~x)  is  only 
possible  If  the  N  phase  drift  functions  ^^(t)  are  known  a  priori  as  a  func¬ 
tion,  of  tlme/range  for  each  of  the  N  frequencies  f^.  Such  functions  could 
be  stored  in  tabular  fora  or  approximated  by  a  polynomial.  Unfortunately, 
these  functions  will  generally  be  different  for  each  of  the  If  frequencies, 

A  measure  of  oscillator  stability  is  the  maximum  value  of  $k(t)  over 
the  expected  soan  of  target  ranges:  r  •  t  ^  .  If  J  $k(t)  | £.  for  all 
t  <  — r(max),  then  the  correction  represented  by  e^k^  can  probably  be 
neglected*  This  will  usually  be  the  case  for  most  crystal  controlled  fre¬ 
quency  sources.  However,  if  this  crlterlor  is  not  sat,  coherent  processing 
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will  only  be  possible  if  ^(t)  can  be  measured  and  calibrated  as  a  function 
of  delay  t  for  each  frequency  f^,  and  the  required  correction  applied  in 
eqn  (11)  as  a  function  of  k  and  t.  Otherwise,  one  cannot  obtain  the 
target’s  high-resolution  range  profile. 

2.2.1  DFT  Implementation  of  Coherent  Processing 
The  following  expression  for  the  high  resolution  range  profile  of  a 
target,  was  derived  in  the  preceding  Section  (eqn  (2-11))  for  the  case  of 
no  non-linear  differential  phase  shifts  as  a  function  of  frequency  (i.e. 
^k(c)  *  0),  and  zero  radial  velocity. 

a— l 

q(x,t)  -  l  V(f k,t)  exp(J2xf)cx)  Af  (2-12) 

k-0 

where:  x  is  the  range  to  the  target  in  units  of  round  trip  delay 

■  time 

t  is  the  time  instant,  relative  to  pulse  transmit  time,  at 

which  the  returns  are  sampled 

V(ffc,t)  is  the  complex  signal  voltage  sample  measured  at 
time,  t,  in  response  to  a  pulse  transmitted  at 
frequency  ffc. 

fit  is  the  carrier  frequency  at  which  the  kth  pulse  in  the 
train  is  transmitted 

q(x,t)  is  the  complex,  pulse-shaped  weighted,  high  resolu¬ 
tion  range  profile  of  the  distributed  target  as  a 
function  of  range  x.  The  resolution  of  this  pro¬ 
file  is  equal  to  the  reciprocal  of  the  waveform 
bandwidth,  NAf , 

N  is  the  number  of  pulses  in  the  train. 

Af  is  the  frequency  spacing  between  pulses  in  H* 

It  should  be  noted  that  the  reference  waveform  used  for  the  coherent 
detector  . to  obtain  V(f^,t)  at  each  frequency  is  the  pulse  carrier 
transmitted  at  that  frequency.  Thus  a  different  reference  signal  is  used 
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for  each  pulse  return.  Therefore,  the  different  carrier  frequencies  used 
must  each  be  stable  but  they  do  not  necessarily  have  to  have  known  relative 
phase  relationships  relative  to  one  another. 

In  order  to  use  a  DFT  to  implement  eqn.  (2-12),  the  frequency  spacing 
between  the  pulses  in  the  train  must  be  uniform. 


Thus:  fjt  -ffc-i  -Af  for  all  k 


Or:  fk  -f0  +kAf  ;  k  -0,1, 


.N  -  1 


Also  write  the  range  delay  x  as:  x  -Xq  +nAx 


(2-13) 


(2-14) 


Where  Xq  Is  an  arbitrary  constant  range  to  be  determined.  Ax  is  the  range 
resolution  corresponding  to  one  DFT  bin  width,  and  n  is  the  DFT  bln  index 
n  -0,1, . N— 1. 

Substituting  (2-13)  and  (2-14)  into  (.2-12)  gives: 

N— 1 

q(x,t)  -  exp(j2irf0(xo+uAx))  £  V(fit,t)exp( j2xkAf(xo+nAx))  Af  (2-15) 

k-0 

Choose  the  arbitrary  constant  ^  to  be  the  maximum  integral  multiple  of 
1/A  f  that  is  just  less  than  the  minimum  range  to  the  distributed  target 


whose  profile  is  to  be  determined. 


Xq  -M/Af  <xmin;  where  M  is  an  integer 


(2-16) 


Substituting  (2-16)  into  (2-15)  gives  the  following  result  for  the 
magnitude  of  the  target's  high-resolution  range  profile. 


q(x,t)|  -|  l  V(fit,t)exp(j2rknAfAx)  Af 
k-0 


(2-17) 


Note  that  . the  resolution  (bin  width)  in  range  delay,  Ax,  is  equal  .to  the 


reciprocal  of  the  waveform  bandwidth: 
Ax  -1/ (NAf) 


(2-18) 


Substituting  (2-18)  into  (2-17)  gives  the  desired  result  for  the 
realization  cf  eqn  (2-12)  by  means  of  an  N-point  DFT. 

N"*l 

|  q(x,  t)j  •  |  l  V(fjt,t)exp(  j2xkn/N)  Af  |  (2-19) 

k»o 

where  the  range  delay  x  is  defined  by  eqn.  (2-14)  and  (2-16),  and  the  DFT  bin 
width  is  given  by  (2-18).  The  magnltu4e  of  the  output  in  bin  n  of  the  DFT  is 
therefore  the  amplitude  of  the  individual  target  scatterer(s)  located  at 
range  x  -Xq  +nAx,  where  Xq  is  defined  by  eqn.  (2-16).  Thua  the  magnitude  of 
the  DFT  output  represents  the  range  profile  of  the  target  over  a  range 
extent  of  c/2(NAx)  -  c/(2Af)  relative  to  a  starting  range  cxQ/2,  with  range 

a 

resolution  cAx/2  «  c/(2NAf). 


It  should  be  noted  that  the  high  resolution  range-profile  produced  at 
the  output  of  the  DFT  extends  from  range  delay  Xq  to  x„  +  1/Af,  where  Xq  Is 
defined  by  eqn.  (2-16),  and  that  these  limits  are  independent  of  the  actual 
time,  t,  of  the  sampling  Instant,  provided  that  Xq  <  t  <  (xq  +  1/Af).  The, 
time,  t,  of  the  sampling  Instant  affects  only  the  degree  of  attenuation  of 


the  target  range  profile  by  the  pulse  shape,  but  has  no  affect  on  the  loca¬ 
tion  of  the  profile  in  the  DFT  output.  This  is  illustrated  in  Fig.  2-3. 

Two  point  scatterers  K  and  B  of  equal  amplitude  are  located  at  range  delays 
x g  and  xg  such  that  Xq  <  x^  <  xg  <  (xq  +  1/Af),  as  Illustrated  In  Fig. 

2 -3a.  If  the  sampling  instant,  t,  is  coincident  with  range  delay  x^  (t  * 
x^)  the  resultant  high  resolution  range  profile  will  appear  as  shown  In 
Fig.  2-3b.  However,  if  the  sampling  time  is  delayed  such  that  t  *  xg,  the 
resulting  high  resolution  profile  will  appear  as  in  Fig.  2-3c.  Note  that 
changing  the  sampling  Instant  from  t  •  xg  to  t  *  xg  has'  only  changed  the 
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RANGE  EXTENT  OF  FFT 
OUTPUT  AX  =  1/Af 


PULSE  RETURN  A/D  SAMPLING  STROBE 


Fig.  2-3  Effect  of  A/D  sampling  time  on  the  high' resolution  range  profile. 

(a)  A  target  consisting  of  two  point-scatterers  located  within  the  same  single- 
pulse  resolution  cell,  and  the  component  pulse  returns  from  each  of  the  two 
scatterers. 

(b)  The  Ugh  resolution  range  profile  at  the  output  of  the  DFT  for  the  target 
configuration  of  Fig.  2-3a  when  the  A/D  sampling  strobe  occurs  at  range  delay  x^« 

(c)  The  high  resolution  ranee  profile  at  the  output  of  the  DFT  for  the  target 
configuration  of  Fig.  2-3a  when  the  A/D  sampling  strobe  occurs  at  range  delay  Xg» 


apparent  amplitudes  of  the  two  scatterers  A  and  B  in  the  high  resolution 
range  profile,  but  has  left  their  locations  in  the  profile  unchanged. 

2 . 2 . 1 . 1  Determination  of  Frequency  Step  Size  and 
A/D  Sampling  Times 

The  choice  of  frequency  step  size,  Af,  and  A/D  sampling  times,  t,  is 
dictated  by  the  need  to  avoid  corruption  of  the  high  resolution  range  pro¬ 
file  of  a  single  range  cell  (i.e.  one  pulse  width)  by  either  aliasing 
(folding)  of  the  returns  from  scatterers  within  that  cell  or  from  scat¬ 
terers  located  in  adjacent  range  cells.  Corruption  cf  the  profile  by 
aliasing  of  the  returns  from  within  that  cell  is  caused  by  choosing  too 
large  a  value  of  Af.  As  pointed  out  in  the  preceding  section,  the  range 
extent  of  the  profile  at  the  output  of  the  DFT  is  equal  to  1/Af.  In  order 
to  avoid  folding  of  the  returns  from  scatterers  within  that  single  pulse 
range  cell,  the  range  extent  of  the  profile,  1/Af,  must  be  made  greater 
than  two  times  the  width  of  a  single  pulse,  i.e.  1/Af  >  2r  -  so  that 
returns  from  scatterers  located  more  than  one  pulse  width  away  from  the 
sampling  time  (the  center  of  the  pulse)  will  be  so  strongly  attenuated  that 
the  amplitudes  of  their  aliased  replicas  become  negligible.'  Thus  the  cri¬ 
terion  on  Af  required  to  prevent  corruption  of  the  profile  by  self-aliasing 
is  Af  <  V2T . 

In  addition,  the  A/D  sampling  time,  ty,  for  the  ath  range  cell  must  be 
chosen  such  that  t^  -  (M  +  1/2 )/Af,  in  order  to  prevent  corruption  of  that 
ceira  high  resolution  range  profile  by  returns  from  adjacent  range  cells. 
Note  that  this  choice  of  t*j  places  tM  in  the  middle  of  the  range  extent  x© 

<  (*o  ♦  VAf)  covered  by  the  output  of  the  DFT.  This  impliec  that  the  *A/D 
sampling  strobes  must  be  synchronized  with  the  variable  frequency  source 


(i.e.,  coherent  frequency  synthesizer)  that  is  used  to  generate  the 
carriers  for  the  frequency  stepped  pulse  train,  and  that  the  frequency 
step  size,  A f ,  oust  be  uniform  over  the  entire  pulse  train.  Recall  that  a 
uniform  frequency  step  size  was  also  a  requirement  for  using  an  DFT  to 
generate  the  high  resolution  range  profiles. 

The  effect  of  aliasing  caused  by  using  too  large  a  frequency  step  size 
is  Illustrated  in  Fig.  2-4.  Two  point  scatterers  A,  B,  of  equal  amplitude 
and  separated  by  less  than  a  pulse  width  are  located  in  adjacent  DFT  range 
intervals,  1/Af,  when  Af  is  chosen  to  be  equal  to  the  reciprocal  of  the 
pulse  width,  -  so  that  the  extent  of  each  range  interval  spanned  by  the 

output  of  the  DFT  is  equal  to  one  pulse  width.  The  range  of  target  B  is 

greater  than  that  of  target  A,  as  shown  in  Fig.  2-4a,  i.e.,  xg  >  x^. 

Although  A  and  B  appear  in  separate  DFT  range  intervals,  1/Af,  their  actual 
physical  separation,  |xg-x^|  is  less  than  a  pulse  width.  For  the  case 
illustrated  in  Fig.  2-4a  |*b“xaI  *  t/2.  Now  if  the  sampling  instant,  t, 
coincides  with  the  midpoint  of  the  first  DFT  range  Interval,  the  resulting 
high  resolution  profile  will  appear  as  illustrated,  In  Fig.  2-4b,  with  the 
return  from  scatterer  B  in  the  next  range  interval  appearing  before  that  of 

scatterer  A.  This  problem  can  be  avoided  by  using  a  smaller  value  of  Af  so 

that  both  scatterers  appear  in  the  same  range  Interval  as  Illustrated  in 
Fig.  2-4c.  Usually  a  value  of  Af  •  1/2*  will  suffice. 

Even  if  Af  is  chosen  as  described  above,  an  ambiguous  rspge  profile 
can  still  result  if  the  sampling  strcbe  is  located  too  far  sway  from  the 
midpoint  of  the  range  delay  Interval  1/Af.  This  effect  is  illustrated  in 
Fig.  2-5.  Again,  two  point-scatterers  A  and  B  are  located  in  adjacent  range 
intervals  as  shown.  Af  is  chosen  such  that  the  range  (delay)  extent,  1/Af, 
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.  WITH  REOUCEO  FPr-fUENCY 
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Tig.  H  Allaalag  of  Clw  high  raaolutloa  rango  prof! la  that  raaulta  f roa  uaa  of 
Coo  largo  a  frequency  atep-alie. 

(a)  Too  poiat-acactarara  aaparated  by  laa*.  th,  -.a  pulaa  width,  but  located  ao 
aa  to  atraddle  an  uaasblguoua  rango  intorval  boundary,  a ad  thoir  corraapoading 
coaponeat  pulaa  raturaa,  uhan  4f«l/i. 

<b)  The  high  raaolutloa  range  profila  obtalnad  at  tha  output  of  tha  DPT  for  tha 
acattarar  eoaf iguratioo  ot  Pig.  2-*a,  ahowing  tha  aliaaad  raapoaaa  of  acattarar 
I  at  xg  “*j-t .  Saapliag  occura  ia  tha  aiddla  of  tha  high  raaolutloa  rango 
Interval,  1/4  f. 

(«)  Tha  high  raaolutloa  ranga  profila  for  tha  acattarar  configuration  of  Pig. 

2 -4a,  ahowlng  how  tha  correct  (unallaaad)  rapponaa  fro*  acattarar  >  la  obtalnad 
by  reducing  tha  frequency  atep-eiac  Af • 1/ (2t >•  Soap ling  ia  again  at  tha  aiddla 
of  tha  rango  internal,  I/A  f. 
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A/D  SAMPLING  STROBE 
AT  tj  »k0*  1  /(2£ft 


Pig.  2-5  Allaalog  of  tha  high  raaolutlon  rang*  proflla  that  raaulca  froa 
lapropar  tlalng  of  tha  A/D  *aap ling  acrobat. 

(a)  Two  polnt-oeactarara  aapararad  by  lata  chan  ona  pulaa  width,  but  locatad  ao 
aa  Co  acraddla  an  unaablpioua  ranga  lncarval  boundary,  and' Chair  corraapondlng 
c'oaponant  pulaa  racurna.  Af-l/(2r),  hut  with  tha  aaapllng  tnacanc  loratad  naar 
tha  ranga  lncarval  boundary. 

(b)  Tha  high  raaolutlon  ranga  proflla  obtaload  at  tha  output  of  tha  DPT  for  tha 

aeattarar  configuration  of  Pig.  2-Ja,  ahowlng  tha  allaaad  raaponaa  of  acattarar 

S  at  xgWg-2r  dua  to  tha  aaapllng  atroba  not  balng  cancarad  to  tha  ranga  tntar- 
val,  1/4  f-2r  . 

(c)  Tha  high  raaolutlon  ranga  proflla  obtalnad  at  tna  output  of  tha’  DPT  for  tha 

acattarar  configuration  of  Pig.  2-5a,  ahowlng  auppraaalon  of  tha  altaaad 

raaponaa  dua  to  acattarar  I  whan  tha  aaapllng  atroba  la  cancarad  In  tha  ranga 
lntarvalt  t#«  «0a  1/ < 2A f )  •  Xg+v. 
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of  each  interval  is  equal  to  2r  in  accordance  vith  the  preceding  criterion. 
H*.**  ver,  if  the  saapling  strobe  is  located  near  the  end  of  the  fir3t  inter¬ 
val  as  shown,  the  resulting  range  profile  will  appear  as  shown  in  Fig.  2-5, 
in  which  the  return  froa  scatterer  B  appears  before  that  fron  A.  However, 

I  ' 

if  the  sampling  instant,  t,  were  aoved  back  to  the  center  of  the  first 
range  interval,  the  "folded”  return  froa  scatterer  B  would  effectively 
disappear  because  of  attenuation  by  the  pulse  shape. 

As  discussed  above,  aabigultle?  and  aliasing  of  the  high  resolution 
range  profiles  can  be  avoided  by  choosing  the  frequency  step  size,  Af,  to 
be  no  greater  than  one-half  the  single  pulse  bandwidth,  and  by  synchro¬ 
nizing  the  tiaing  of  the  A/0  saapling  strobes  to  the  frequency  step  size, 
Af,  in  the  frequency  synthesizer  so  that  all  of  the  strobes  will  appear  . 
exactly  in  the  aiddle  of  each  range  Interval,  as  shown  in  Fig.  2-5c. 

However,  if  the  strobes  are  synchronized  to  Af  and  if  Af  -  V2f,  the 
strobes  illustrated  in  Fig.  2-5  will  be  two  pulse  widths  apart.  As  a 
result,  scatterers  located  in  the  vicinity  of  the  boundary  between  adjacent 
range  intervals,  such  as  the  point  scatterer  A  in  Fig.  2-6,  will  appear  in 
the  range  profile  heavily  attenuated  by  the  skirts  of  the  pulse  shape.  The 
solution  to  this  problem  is  to  place  «  second  set  of  saapling  strobes 
synchronized  to  A {between  the  first  set,  l.e.  at  the  range  Interval  boun¬ 
daries,  In  order  to  generate  a  second  set  of  high  resolution  range  profiles 
that  overlaps  the  original  set.  This  second  set  of  saapling  strobes  is 
indicated  by  the  dashed  lines  with  crosses  in  Fig.  2-6c  and  d.  However, 
the  high  resolution  profiles  derived  froa  this  intermediate  set  of  saapling 
strobes  will  appear  aliased  as  In  the  case  of  the  return  froa  scatterer  B 
in  Fig.  2-bd,  Hence  this  fact  aust  be  taken  into  account  in  sashing  those 
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ALIASEO  RESPONSE 
FROM  SCATTERER  B 
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HIGH  RESOLUTION  RESPONSE  OBTAINED 
FROM  FFT  OF  INTERMEDIATE  SAMPLING 
STROBE  SAMPLES  AMPLITUDE  RESPONSE 
OF  SCATTERS  A  AND  B  IS 
RESTORED  RUT  THAT  OF 
SCATTERER  8  IS  ALIASED 


Fig.  2-b  Illustration  of  the  Mod  for  in  Intermediate  »«c  of  sampling 
strobes  at  the  rang*— lncarval  boundaries,  t§i“k/4f,  to  addition  to  tha 
principal  aat  at  tha  ranga-lntarval  aldpolntas  t#*k/4f  +1/(240  co  provide 


coaplet 
(a)  A 


coverage  between  adjaeant  ranga  Intervale, 
jingle  point-scat tarer  loeatad  naar  tha  edge  of  a  ranga  Interval  and  lea 


pulaa  riturn. 

(b)  Tha  high  resolution  ranga  profile  at  tha  output  of  tha  OFT  In  response  to 
tha  aca:tarar  coaf lgu ration  of  Fig.  2-ba  showing  tha  aavara  attenuation  of  tha 
reeponei  from  aeattarar  A  whan  tha  saapllng  atroba  la  centered  In  tha  range 
Interval;  Af-l/(2r). 

(e)  Two  polnc-acatterara  loeatad  In  adjaeant  ranga  tntervala  near  a  range  Inter¬ 
val  bouidary,  with  Intermediate  sampling  atrobaa  at  tha  ranga  Interval  boun¬ 
darlaa. 

(d)  Tha  high  raaolutlon  range  proflla  at  tha  output  of  tha  DFT  of  tha  lnter- 
nadlata  aaaplaa  for  tha  target  configuration  of  Fig.  2-6c,  showing  the  normal 
response  of  aeattarar  A  and  tha  expected  eltaeed  reaponaa  of  aeattarar  3.  The 
unsllasid  reaponaa  can  be  obtained  bp  Interchanging  tha  flrat  and  seeond  halvea 
of  tha  DFT  Output.  Tha  midpoint  of  tha  raaultlng  proflla  la  at  Xg+i/Af. 
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profiles  with  the  original  ones.  Specifically,  when  processing  the  returns 
sampled  by  the  intermediate  sampling  strobes  located  at  the  range  interval 
boundaries,  k/Af,  the  first  (near)  half  of  the  resulting  high  resolution 
range  profile  must  be  associated  with  the  (far)  range  interval  to  the  right 
of  the  strobe  while  the  second  (far)  half  of  the  high  resolution  profile  is 
associated  with  the  (near)  range  interval  to  the  left  of  the  strobe  loca- 
tion. 

2.2. 1.2  Effect  of  Non-Uniform  Frequency  Spacing  ■' 

In  the  preceding  development,  it  has  been  assumed  that  the  frequency 
spacing,  Af,  between  the  pulses  in  the  stepped  frequency  pulse  train  was 
uniform  (eqn.  2-13).  This  assumption  resulted  in  the  high  resolution  range 
proflle|  q(x,t)|  being  determined  as  the  OFT  of  the  complex  frequency 
samples  (V(f^,t),  eqn.  (2-19).  In  this  section  the  effect  of  non-uniform 
frequency  spacing  on  the  computation  of j  q(x,t)j  will  be  determined. 
Specifically,  the  maximum  allowable  deviation  between  the  actual  frequency 
at  which  the  kth  pulse  is  transmitted,  f*  and  the  frequency  (fo  +  kAf), 
corresponding  to  a  uniform  frequency  spacing,  that  can  be  tolerated  in 
using  eqn.  (2-19)  to  compute)  q(x,t)j  will  be  determined. 

The  assumption  of  uniform  frequency  spacing  between  the  pulses  In  the 
s tepped-f requency  pulse  train  Is  often  an  Idealization,  even  when  the  Indi¬ 
vidual  frequencies  transmitted  are  absolutely  stable.  For  example,  suppose 
that  the  coherent  frequency  synthesizer  is  realized  using  a  number  of  dif¬ 
ferent  crystal  oscillators,  each  of  which  is  absolutely  stable.  Such  a 
system  could  be  realized  either  by  using  a  separate  crystal  to  generate 
each  of  the  N  frequencies,  or  by  using  K  separate  crystals,  where  M«N,  to 
generate  the  N  frequencies  using  a  network  of  frequency  mixers  and 


multipliers.  In  general,  the  spacings  between  adjacent  frequencies 
generated  by  such  a  synthesizer  will  not  be  .precisely  uniform  because  there 
is  an  error  tolerance  associated  with  each  crystal. 

As  previously  discust>od  in  Section  2.1.1,  the  complex  voltage  return, 
V(fk,t),  received  from  a  distributed  target  will  in  general  be  a  function  of 
the  transmitted  frequency,  fk,  in  accordance  with  eqn.  (2-7). 

V(fk,t)  J  p(x)s(t-x)e’"^2lI^t$x 


Let  fk  represent  the  ideal  values  of  transmitted  frequency, 
corresponding  to  a  uniform  frequency  spacing,  Affc,  as  defined  by  eqn  (2-13). 


fk  ■  f0  +  kAf 


(2-13) 


Where  Af  is  the  average  value  of  the  frequency  spacing  over  N  frequencies: 

Af  -  (f n-i  -  f0)/(H— 1) 

Let  the  actual  frequency  transmitted  on  the  kth  pulse  be  gfe,  where  gfc  is 


related  to  fk  by: 


8k  •  *k  +  «k 


(2-21) 


Thus  ck  represents  the  deviation  of  g^  from  the  ideal  value,  fk. 


Assume  that  the  frequencies  fk  and  g^  are  absolutely  stable  such  that 
^>k(t)  =  0.  Then  using  (2-21),  the  expressions  for  V(fktt)  and  V(gk,t)  can 


be  written  as: 


V(fk,t)  -  /  p(x)s(t-x)e"J2,rfltXdx 


(2-22a) 


v<8k,t>  ■  /  p(x)s(t-*)erj2x«tXdx 


(2-22b) 


Now  make  the  transformation: 


u  •  x-t 


and  recall  that  the  pulse  shape  s(u)  la  such  that  s(u)  •  0  for  u  >2t  . 
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Us’ng  the  above  transformation  and  the  relationship  between  and  fk 
given  by  (2-21,  eqns  (2-22a,b)  may  be  rewritten  as  follows: 


V(fk,t)  "  e  J  p(t+u)s(-u)e  ^^^“du 

-2t 


(2-23a) 


v(8k,t)  -  e"J2ir(fk  +  6k)t  /  p(t+a)s(-u)e‘J2irfkVj2lirekUdu  (2-23b) 

-2t 

* 

It  is  reasonable  to  assume  that  the  frequency  deviation,  ek,  between 
gk  and  fk  will  be  a  small  fraction  of  the  nominal  frequency  spacing  Af  i.e. 


Assume: 


|  ekJ  <  0.01  Af 


(2-24) 


Furthermore,  in  accordance  with  Section  2.2. 1.1,  the  upper  bound  on  Af  is 
given  by:  Af  <  -i- 


Combining  the  above  two  criteria  gives: 

|  ek|  t  <  .005  (2-25) 

Consequently,  the  term  e"J2ireku  *  i  in  the  Integrand  of  eqn.  (2-23b) 
so  that  the  two  integrals  in  (2-23a,b)  become  approximately  equal.  As  a 
result,  it  follows  that: 

v<8k,0  -  V(fkft)e‘j2,,ekt  (2-26a) 

Now  even  if  ek  is  small  enough  such  that  (2-25)  is  satisfied,  the 
phase  term  e  J2*6^  i„  (2-26a)  will  generally  not  be  negligible  because 
t>>r .  The  maximum  value  of  t  is  equal  to  the  radar  pri  which  is  usually  at 
least  three  ordern  of  magnitude  greater  than  the  pulse  width  t. 

In  general,  the  complex  returns  V<g|t#t);  received  from  pulses  transmitted 
at  frequencies  gj<  must  be  corrected  for  the  phase  shift,  2*ekt,  due  to  the 
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frequency  spacing  non-uniformity  e^,  before  the  high  resolution  range  pro¬ 
file  can  be  computed  using  eqn  (2-19).  Thus  V(fjt>t)  is  computed  from  the 
measurements  V(gil>t)  as  follows: 

V(fk,t)  -  VCg^Oe^*^  (2-26b) 

Note  that  the  required  phase  shift  correction  is  generally  a  function 
of  both  the  frequency  step  index,  k,  and  the  approximate  range  delay  to  the 
target,  t.  The  effect  of  not  making  this  correction  is  the  same  as  if  the 
ideal  values  V(fk,t)  had  been  measured  and  then  corrupted  with  phase  noise 
of  amplitude  2we{ct  prior  to  computation  of  the  range  profile.  One  of  the 
effects  of  such  phase  noise  is  to  raise  the  sidelobe  level  of  the  resulting 
range  profile.  The  variance  of  the  phase  noise  (zero-mean)  added  to  a  pure 
sine  wave  can  be  shown  to  be  equal  to  the  integral  of  the  relative  noise 
power  in  the  sidelobes  of  the  resulting  power  spectrum.  By  Parseval's 
Theorem  for  a  discrete  N-point  spectrum: 

.  N-l  ,  N-l 

"♦2  -  jf  1  *2<->  ■  4  l  1*00 12  <2-27) 

n-0  k-0 

where  $(k)  is  the  spectrum  of  the  phase  noise  $(n),  and  where  the  phase 
noise  is  assumed  to  be  small  enough  so  that  the  following  approximation 
holds : 

ej*  -  l+j$;  !♦  ]  _<  0.5  rad.  (2-28) 

The  term  N^  in  the  denominator  of  (2-27)  is  the  square  of  the  amplitude  of 
the  spectral  line  in  the  DFT  of  a  unit  amplitude  sine  wave;  •  Hence 

(k) { 2/n2  i8  the  relative  sidelobe  (power)  level  in  DFT  bin  k  relative  to 
the  signal  power  (in  the  bin  containing  the  signal).  It  is  important  to 
note  that  (2-27)  is  a  relation  between  the  phase  noise  level  and  the 
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corresponding  integrated  sldelobe  response.  However,  for  purposes  of  reso¬ 
lution,  the  peak  sidelobe  level  is  usually  of  greater  importance  than  the 
integrated  level.  The  sidelobe  spectrum  is.  often  not  flat,  so  that  the 
highest  levels  usually  occur  in  the  vicinity  of  the  main  lobe.  In  general, 
it  is  not  possible  to  determine  the  peak  sidelobe  level  for  a  given  value 
of  rms  phase  noise  unless  the  spectrum  of  this  noise  is  also  known. 

However,  for  any  given  value  of  rms  phase  noise  the  minimum  value  of  the 
peak  sidelobe  level  is  obtained  when  the  phase  noise  sidelobe  spectrum  is 
flat.  In  this  special  case,  the  relative  sidelobe  level  is  a  constant 
independent  of  k  and  is  given  by: 

rt  m  —  —  m  ■■  1  — 

P  N2  N 

However,  the  assumption  of  a  flat  phase  noise  spectrum  is  probably  unduly 
optimistic,  especially  if  the  frequency  errors  of  the  synthesizer  have  a 
systematic  component  -  i.e. ,  quadratic  or  sinusoidal  as  a  function  of  fre¬ 
quency.  A  more  reasonable  approximation  might  be  to  assume  that  the  shape 
of  phase  noise  power  spectrum  is  a  maximum  at  the  signal  frequency  and 
quadratically  approaches  zero  for  k-N,  as  illustrated  in  Fig.  2*7.  In 
this  case,  the  sidelobe  level  decreases  monotonlcally  from  the  vicinity  of 
the  main  lobe  as  shown  in  Fig.  2-7,  and  the  peak  value  of  the  relative 
sidelobe  level  adjacent  to  the  main  lobe  is  given  by: 


p(max)  -  3 


(2-29) 


Assuming  that  the  naximum  sidelobe  level  is  to  be  no  greater  than  -30  dB 


gives  the  following  bound  on  the  rms  phase  error. 
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fig*  2-7  Phase  noise  spectrum  for  a  sinusoidal  signal,  showing  effect  of  . 
spectrum  shape  on  peak  sidelobe  level.  Spectrum  amplitude  normalized  with 
respect  to  signal  peak  at  signal  frequency. 


39 


a*  <  /  -sSOlN  .  #018  yfj* 


The  corresponding  bound  on  the  rms  frequency  error  is  therefore  given  by: 
2x  lekl  cmax  *  j(  .018/n” 


,  ,  „  .0029/n" 

112 


(2-30) 


where  tMX  is  the  maximum  value  of  the  round  trip  range  delay  to  the 
target  and  N  is  the  number  of  pulses  in  the  train.  In  oust  cases,  tMX 
will  be  equal  to  the  radar  pri.  Note  that  since  t^x  »  r,  if  the  magni¬ 
tude  of  the  frequency  deviations,  |e'£|,  is  bounded  by  eqn.  2-30,  the  basic 
bound  on  of  eqn.  2-24,25,  that  was  assumed  for  the  analysis,  will  auto¬ 
matically  also  be  satisfied. 

2.2.2  Interpretation  of  DFT  Processing  as  Matched  Filtering 
In  the  following,  the  DFT  is  shown  to  be  the  matched  filter  for  a 
frequency-s tepped  pulse  train  with  uniform  frequency  spacing,  provided 
that  there  is  no  relative  radial  motion  between  the  target  and  the  radar. 

Consider  the  return  received  from  a  point  target  located  at  range 
delay  x,  in  response  to  a  pulse  transmitted  at  frequency  f^«  The  phase  of 


the  received  pulse  return  from  the  target  is  then  given  by: 

zjt  *»2irfkX 


(2-31) 


where  the  phase  indicated  by  eqn.  (2-31)  is  measured  relative  to  the  trans¬ 
mitted  carrier  frequency.  In  other  words,  the  coherent  reference  signal 
for  each  pulse  return  is  the  carrier  of  the  pulse  transmitted  at  that  fre¬ 
quency.  The  complex  return  from  this  point  target  may  then  be  written  as: 


V(f|t,f.)  •s(t-x)exp(-j2irficx) 


(2-32) 
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Where  s(t-jc)  is  Che  pulse  shape  of  Che  CransmiCCed  signal.  As  before, 
Che  range  x  may  be  written  written  as  x  *xQ  +Ax,  where  x  ;  defined  by 
eqn.  (2-16).  Then  using  eqn.  (2-13)  and  Che  face  ChaC  AfXg  -M,  eqn.  (2-33) 
can  be  wriCCen  aa: 

V(fk,c)  =*s(t-x)exp(-j2irf0x)exp(-j2irkAf6x)  (2-33) 

The  macched  filCer  for  this  sequence  of  samples,  with  respecc  Co  Che 
sample  index  k,  is  iCs  complex  conjugaCe.  Thus  it  Is  seen  that  the  kernel: 

exp(j2irknAfAx)  (2-34) 

of  the  DFT  represented  by  eqn.  (2-17)  is  the  matched  filter  for  the  signal 
samples  V(fk,t)  when  nAx  *<$x.  Substituting  Che  expression  for  V(fk,t) 
given  by  (2-33)  into  (2-17)  gives: 

N— 1 

|  q(x,  t)j  -  |  •  £  .  8(t-x)  exp(-j2nkAf(nAx-dx))  Af  J  (2-35) 

k-0 

The  maximum  value  of  (2-24)  clearly  occurs  when  nAx  -Ax,  as  i$  required 
for  a  matched  filter.  A  simple  explanation  of  why  a  DFT  is  the  matched 
filter  for  this  pulse  train  waveform  follows:  When  the  returns  from  the 
frequency-stepped  pulse  train  are  coherently  detected  using  each 
transmitted  carrier  as  the  corresponding  coherent  reference,  the  phase  of 
the  detected  returns  from  a  fixed  point  target  will  vary  linearly  with  fre¬ 
quency,  as  previously  indicated  by  eqn*  (2-31).  Similarly  the  phase  of  the 
kernel  of  the  DFT  defined  by- eqn.  (2*17)  also  varies  linearly  with  frequency 
for  a  fixed  value  of  range  delay  x.  Therefore  since  the  phases  of  the 
signal  and  the  filtering  kernel  both  vary  in  the  same  way  as  a  function  of 
sample  index,  k,  the  kernel  is  Che  matched  filter  for  the  signal. 

It  should  be  noted  that  the  fact  that  the  DFT  is  the  matched  filter 
for  the  returns  from  a  train  of  frequency  stepped  pulses,  with  uniform 
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frequency  spacing,  is  not  applicable  to  the  case  of  conventional  digital  pulse 
compression,  where  the  return  from  a  single  swept-f requency  pulse  is 
'  uniformly  sampled  over  an  interval  equal  to  the  duration  of  the  pulse. 

Although  these  two  situations  appear  to  be  similar,  in  that  in  both  cases 
signal  samples  of  the  pulse  returns  are  taken  at  uniformly  spaced  points  in 
carrier  frequency,  there  is  a  subtle  but  important  difference  between  the 
two.  In  the  case  of  the  stepped  frequency  pulse  train,  as  mentioned 
before,  the  detector  reference  signal  for  each  pulse  return  is  the 
transmitted  pulse  carrier  at  that  frequency,  which  changes  from  pulse  to 
pulse.  As  a  result,  the  phases  of  the  detected  returns  from  a  fixed  point 
target  will  vary  linearly  with  frequency.  On  the  other  hand,  in  the  case 
of  the  uniformly  sampled  swept-f requency  pulse,  the  detector  reference 
signal  used  is  a  constant  frequency  waveform.  Consequently,  the  phases  of 
the  detected  signal  samples  will  vary  quadratlcally  with  the  Instantaneous 
frequency  at  each  sampling  point,  rather  than  linearly.  Therefore,  the  DFT 
will  no  longer  be  a  matched  filter  for  the  resulting  signal  samples. 

2.2.3  Effect  of  Radial  Velocity 

In  the  preceding  development  it  was  assumed  that  the  range  between  the 
radar  and  the  target  remained  constant  over  the  duration  of  the  pulse 
train.  The  purpose  of  this  section  is  to  determine  the  effect  of  finite 
relative  radial  velocity  on  the  output  of  the  DFT  matched  filter.  It  does 
not  matter  whether  the  radial  velocity  is  due  to  motion  of  the  target,  the 
radar,  or  both. 

An  uncompensated  nonzero  radial  velocity  is  found  to  have  two  effects 
on  the  range  profile  at  the  output  of  the  DFT.  The  first  is  a  shift  in 
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the  range  profile  of  the  target.  This  is  simply  the  well-known  range- 
doppler  coupling  effect  associated  with  the  use  of  linearly  frequency-swept 
waveforms.  As  will  be  seen,  even  small  values  of  radial  velocity  can  cause 
significant  shifts  in  the  range  profile.  If  the  magnitude  of  the  range 
shift  is  greater  than  tne  range  extent,  c/2Af,  spanned  by  the  output  of  the 
DFT,  the  target  range  profile  will  be  circularly-shifted  and  appear 
"wrapped-a round"  in  the  output  of  the  DFT.  The  second  effect  is  that  of  a 
processing  loss  and  a  degradation  in  resolution  due  to  the  fact  that  DFT  is 
really  only  matched  to  the  signal  samples  for  the  case  of  zero  radial  velo¬ 
city. 

Recall  that  the  complex  return  from  a  single  pulse  transmitted  at  fre¬ 
quency  ffc  is  V(fk,t)  given  by  eqn.  (2-7).  Now  suppose  that  the  relative 
radial  velocity  between  the  radar  and  the  target  is  vr  m/sec.  The  prin¬ 
cipal  effect  of  this  motion  on  the  complex  returns  V(fk,t)  received  from  the 
target  is  to  introduce  a  relative  phase  shift  between  the  returns  received 
from  successive  pulses  in  the  train.  Let  T  represent  the  pulse, repetition 
interval  of  the  pulse  train.  Hence  the  kth  pulse  in  the  train  is 
transmitted  at  frequency  fk  at  time  kT,  relative  to  the  time  of  the  first 
pulse  at  k*0.  The  relative  target  motion  during  this  .interval  is  then  just 
vrkT,  so  that  the  kch  pulse  return  will  undergo  an  additional  phase  shift 
relative  to  that  of  the  first  pulse  of: 

2vrkT 

*k  “  2ir  fk( — J — >  (2-36) 

Therefore  if  the  complex  return  from  the  target,  in  the  absence  of  the 
relative  radial  motion,  Is  given  by  V(fk,t),  the  return  from  the  same 
target  moving  at  velocity  vr  relative  to  the  radar  will  be  given  by: 
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VB(fk,t)  -  V(fk,t)e 


j2*fk(2vrkT/  c) 


(2-37a) 


Clearly,  If  the  radial  velocity  is  known  precisely,  the  returns  Vm(fk,t) 
received  from  the  moving  target  can  be  compensated  for  the  effects  of  rela¬ 
tive  radial  motion  by  removing  this  motion-induced  phase  shift  as  follows: 


V(fk,t)  -  Va(fk,t)e 


-j2*fk(2vrkT/c) 


(2-37b) 


The  resulting  values  of  the  motion  compensated  pulse  returns  V(fk,t) 
are  equivalent  to  those  that  would  have  been  received  from  a  stationary 
target,  and  can  be  processed  to  obtain  the  high-resolution  range  profile 
using  the  procedure  described  previously  in  Section  2.2.1. 

however,  if  the  returns  from  a  moving  target  are  not  motion- 
compensated  before  processing  to  obtain  the  high-resolution  range  profile, 
the  resulting  profile  will  appear  shifted  In  range,  as  well  as  suffering 
some  attenuation  and  dispersion  due  to  the  fact  that  the  DPT  will  no  longer 
be  a  matched  filter  for  the  returns  from  the  moving  target.  To  see  this, 
substitute  eqn.  (2-37a)  for  the  moving  target  returns  into  eqn.  (2-19)  for 
the  DPT.  The  resulting  computed  range  profile  Is  then  given  by: 

i  j  Nr  1  kn 

q  (x.q  -  l  V(fk,t)exp(J2ir(2fkvrkT/c  +  ■=£])  (2-38) 

k-0  ™ 

It  is  now  advantageous  to  define  the  following  dimensionless  parameter: 


2vrNT  vrNT 
cdx  Ar 


(2-39) 


Hence  P  Is  the  ratio  of  the  amount  of  relative  radial  target  motion  over 
the  duration,  NT,  of  the  pulse  train,  to  the  resolution,  Ar,  of  the  high 
resolution  range  profile.  Substituting  this  definition  for  P  into  (2-39) 
and  using  (2-13)  gives: 
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j  q(x,t)j  -  |  l  V(fk,t)exp[ j2*((f0  +  k.Af)PAx  +  n)k/N]  j 
Finally,  using  eqn.  (2-18)  and  defining  the  total  bandwidth  of  the  pulse 
train  as: 

B«,NAf  (2-40) 

gives 

|q(x,t)|-  ^\(fk,t)exp[J2»((^§  +  |)P  +  n>|}  (2-41) 

.  where  x  »  Xq  +  nAx 

Note  that  eqn.  (2-41)  reduces  to  eqn.  (2-19)  for  the  case  of  zero  radial 

velocity  (P«0).  The  principal  effect  of  uncompensated  radial  motion  is  to 

shift  the  apparent  location  of  the  target  in  the'  output  of  the  DFT  by 
fQ  i 

(-g-  +  -j)P  fine  range  resolution  cells  (l.e:,  DFT  bins),  or  equivalently 
a  delay  shift  of  (f0/B  +  l/2)PAx  seconds.  This  is  obvious  from  an  examlna- 
ton  of  the  argument  of  the  complex  exponential  in  eqn.  (2-41)*  The  term 
f0/B  is  the  ratio  of  the  carrier  frequency  of  the  first  pulse  in  the  train 
to  the  pulse  train  bandwidth,  B.  Consequently,  the  term  (f0/B  +V2)  is  the 
ratio  of  the  center  frequency  of  the  pulse  train  to  the  pulse  train  band¬ 
width.  The  range  shift  of  the  target  profile,  in  terms  of  the  number  of 
high  resolution  cells,  Ax,  (or  FFT  bins)  is  therefore  given  by: 

L  -  (f C/B)P  (2-42) 

where  fc  =  f0  +  B/2,  Che  center  frequency  of  the  pulse  train. 

Recall  that  P  represents  the  number  of  fine  range  resolution  cells  that  the 
target  moves  relative  Co  the  radar,  over  the  duration,  NT,  of  the  pulse 
train.  When  f0/3  is  large,  even  small  values  of  P  can  result  in  a  signifi¬ 
cant  range  shift  at  the  output  of  the  DFT  of  the  high  resolution  range 
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profile.  Furthermore,  the  resulting  shift  of  the  profile  is  circular,  modulo 
the  range  extent  represented  by  the  output  of  the  DFT.  For  example,  suppose 
the  DFT  output  represents  a  span  of  64  range  cells.  Then  if  (fc/B)  m  33  and 
P  ■  3,  the  range  profile  will  appear  shifted  by  35  range  cells  ■  (3  x  33  -  64) 
The  result  of  this  shifting  is  to  cause  a  possible  ambiguity  in  the  interpre¬ 
tation  of  the  range  profile  as  previously  explained  in  Section  2.2. 1.1. 

In  addition  to  this  velocity-induced  range  shift,  the  quadratic  phase 
term^  P(k/N)2,  in  the  exponent  of  eqn.  (2-41)  causes  .a  dispersion  of  the 
range  profile.  The  magnitude  of  this  dispersion  can  be  determined  by  eva¬ 
luating  eqn.  (2-41)  for  the  case  of  a  point  target,  at  different  values  of 
P.  The  results  are  illustrated  in  Fig.  2-8  for  the  case  of  Hamming 
weighting  and  the  sampling  Instant,  t,  located  at  the  center  of  the  pulse 
return  from  the  point  target.  The  range  shift  induced  by  the  linear  term 
(fc/B)P  does  not  affect  the  shape  of  the  response  and  thus  is  not  portrayed 
in  Fig.  2-8.  It  is  seen  that  serious  degradation  of  the  response  occurs 
for  values  of  P>3.  At  P  ■, 3,  the  peak  of  the  point  target  response  has 
been  attenuated  by  about  3  dB  and  the  width  of  the  reponse  broadened  by  a 
factor  of  2. 

It  is  noteworthy  that  in  most  practical  cases,  the  degree  of  motion 
compensation  required  to  prevent  circular  shifting  of  the  high  resolution 
profile  is  considerably  more  stringent  than  that  required  to  keep  range 
dispersion  of  the  profile  to  a  tolerable  level.  This  is  illustrated  by  thp 
following  simple  example.  Consider  a  stepped-f requency  MMW  radar  having 
the  following  parameters: 
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P  =  O  (No  Dispersion) 
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RELATIVE  RANGE  RESOLUTION  CELL  OR  FFT  BIN  NUMBER 

Pig.  2-8  Attenuation  and  dispersion  of  the  high  resolution  range  prcflle 
of  a  single  polnt-scatterer  target  due  to  uncoopensated  relative  radial 
velocity. 


35  GHz 


fc  ” 

B  -  500  MHz 
N  •  64  pulses 
T  «50  usee 

Ve  wish  to  determine  the  maximum  values  of  uncompensated  radial  velocity 
error,  vr,  that  can  be  tolerated  in  order  to: 

1.  Prevent  profile  shifting  by  more  than  one  high  resolution  range 
cell  (FFT  bin). 

2.  Range  dispersion  of  the  profile. 

He  find  the  maximum  allowable  value  of  the  dimensionless  parameter  P  in 
each  case  and  then  find  the  values  of  the  velocity  error,  vr,  corresponding 
to  those  values  of  P. 

The  amount  of  profile  shift,  in  high  resolution  range  cells  is  given 
by  eqn.  (2-42).  Thus  the  criterion  for  a  range  shift  of  less  than  one  cell 
is: 

fc 

L  «/P  <  1 

Thus,  for  the  given  radar  parameters,  the  limiting  value  of  P  to  prevent 
profile  shifting  is: 

Pmax  <  J  "  ,0lS 

The  corresponding  value  of  vr  is  found  from  the  definition  of  P  given  by 
eqn.  (2-39).  The  range  resolution  Ar  is  a  function  of  the  waveform  band¬ 
width  B  and  Is: 

C 

Ar  *  7?  •  0.3m  for  the  given  value  of  B. 
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Hence  vr(aax)  •  -  1.4  u/aec 

r  NT  (64X5x10-5), 

This  Is  the  maximum  value  of  uncoapensated  radial  velocity  that  can  be 
tolerated  to  avoid  shifting  of  the  profile.  However  the  profile  dispersion 
caused  by  such  a  small  velocity  error  is  negligible. 

On  the  other  hand,  in  order  to  keep  the  effects  of  range  dispersion  of 
the  profile  to  tolerable  levels,  the  value  of  P  should  probably  be  held  to 
less  than  2.  The  corresponding  value  61  uncoapensated  velocity  error  is 
easily  found  to  be: 


vr(max) 


PAr 

NT 


2(0.3) 

(64)(5xl0-5) 


187  a/sec. 


a  tolerance  that  is  two  orders  of  magnitude  higher  than  that  required  to 
avoid  shifting  of  the  profile. 

In  many  cases,  circular  shifting  of.  the  target  range  profile  may  be 
tolerable,  paitlcularly  if  Af  is  chosen  in  accordance  with  the  criteria  of 
section  2.2. 1.1,  so  that  the  occurrence  of  a  circular  shift  in  the  profile 
can  be  identified.  In  such  cases,  the  accuracy  of  the  required  motion  com¬ 
pensation,  need  only  be  sufficient  to  prevent  excessive  dispersion  of  the 
profile.  As  seen  froa  the  preceding  example,  the  corresponding  value  of 
velocity  error  tolerance  can  be  rather  large,  and  in  those  cases  where  the 
radar  platfora  and/or  target  relative  velocity  is  saaller  than  this 
tolerance,  no  notion  coapensation  at  all  aay  be  required. 


2.3  Noncoherent  Processing 


When  coherent  processing  of  the  measured  coaplex  pulse  returns  cannot 


be  performed  because  the  local  oscillator  phase  drift  functions  fk(t)  are 


too  large  and/or  cannot  be  calibrated  with  sufficient  accuracy,  or  a 
noncoherent  radar  is  used,  it  will  not  be  possible  to  obtain  the  pulse-shape 
weighted  high-resolution  range  profile  of  the  target,  q(x,t).  However,  in 
that  event  it  is  still  possible  to  obtain  the  auto-correlation  Rq(z,t)  of 
the  target's  high-resolution  range  profile  by  noncoherently  processing  the 
magnitudes  of  the  frequency-stepped  pulse  returns  J  V(fk,t)|  .  This  is 
accomplished  as  follows*  The  auto-correlation  of  the  target-range  profile 


is  defined  by: 

8q(z,t)  = 


q*(x,t)q(x  +  z,t)dx 


(2-43) 


As  discussed  previously,  variable  t  in  the  above  expression  represents 
the  time  of  the  sampling  instant,  and  hence  may  be  considered  as  merely  a 
parameter. 

Since  Q(f,t)  is  the  transform  of  q(x,t)  it  follows  from  elementary 
Fourier  transform  theory  that  the  transform  of  Bq(r,t)  is  |  Q(f  ,t)|  2. 

Rq(z,t)<— J  Q(f,t)j  2  (2-44) 

Bjit  by  virtue  of  the  relationship  between  Q(f,t)  and  V(f,t)  given  in 
eqn.  (2-8a)  it  also  follows  that 

|  Q(f,t)|  2  -  |  V(f  ,t)|  2  (2-45) 

Hence,  substituting  (2-45)  into  (2-44)  it  can  be  seen  that  Rq(z,t)  can  be 
determined  from  only  the  magnitude  of  V(f ,t)  as  follows: 

Rq(*,t)  -  j  |  V(f,t)(  2  eJ2irfzdf  (2-46) 

•«  • 

Again,  the  discrete  version  of  this  expression,  as  la  the  case  of  eqn.  (2-11) 
Is  given  by: 
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Rq(z.t) 


(2 -46a) 


-  £  I  v<fk.M  2  eJ2irffcz  Af 

k-0 

The  principal  significance  of  this  result  Is  that  computation  of  Rq(z,t) 
does  not  require  any  knowledge  of  the  phase  of  the  complex  return  V(ffc, t) . 
Therefore,  unlike  the  signal  processing  required  to  determine  q(x,t)  In  eqn. 
(2-11),  the  above  result  is  not  affected  by  either  the  phase  drift,  ffc(t), 
of  the  local  oscillator,  any  differential  phase  delays  as  •  function  af 
frequency  in  the  system,  or  target  radial  velocity.  Furthermore,  since 
only  the  magnitude  of  the  complex  return  V(f,t)  need  be,  measured  to  compute 
Bq(z,t),  a  coherent  radar  is  not  required  in  this  case. 

Although  the  computational  effort  required  to  compute  the  auto¬ 
correlation  function  Rq(z,t)  is  nearly  the  same  as  that  required  to  obtain 
the  underlying  high  resolution  range  profile  q(x,t),  the  former  has  the 
advantage  of  being  more  robust  for  the  following  reasons . 

1.  A  coherent  radar  Is  not  required 

2.  No  motion  compensation  is  required 

3.  The  resulting  range  profile  auto-correlation  function  is  not  sub¬ 
ject  to  arbitrary  velocity-induced  range  shifts  as  in  the  case  of 
the  range  profile. 

4.  When  using  a  DFT,  the  tolerance  on  the  deviation  of  the  frequen¬ 
cies  transmitted  vlth  respect  to  uniform  frequency  spacing  is  two 
to  three  orders  of  magnitude  greater  than  In  the  case  of.  coherent 
processing, 

2.3.1  DFT  Implementation  of  Noncoherent  Processing 
The  Implementation  of  eqn  (2-46a)  using  a  DFT  is  similar  to  that  for 
the  case  of  coherent  processing  previously  described  In  Section  2.2.1. 


Again,  Che  pulse  Crain  frequencies  must  be  uniformly  spaced  in  accordance 
with  eqn.  (2-13).  Similarly,  the  range  delay  separation  variable,  z,  of  the 
auto-correlation  function  can  be  written  as  z  »  nAz,  where  Az  is  the  range 
delay  resolution  corresponding  to  the  waveform  bandwidth,  and  as  before,  is 

given  by:  Az  ■ .  (2-47) 

Therefore,  the  required  processing  can  be  implemented  using  a  DFT,  by 
rewriting  eqn.  (2-46a)  as: 

Rq(nAz,t)  -  V  |  V(fk,t)|  2  eJ2irkn/N  Af  (2-48) 

In  general,  Rq(nAz,t)  will  be  complex.  However,  as  before,  the  function  of 
Interest  is  the  magnitude  of  this  auto-correlation  function.  Furthermore, 
since  |  V(f|(, t)|  2  is  real,  the  magnitude  of  its  transform  will  be  an  even 
function  of  n: 

|  Rq(nAz,t)j  -J  Rq(-nAz,t)| 

or  since  the  OFT  is  circular,  it  then  also  follows  that: 

J  Rq(nAz,t)j  -J  Rq(N-n)Az,t)j 

Consequently,  unless  |  Rq(nAz,t)j  *  0  for  n  >  N/2,  an  aliased  response  will 
be  obtained  as  Illustrated  in  Fig.  2-9.  Note  that  the  extent,  NAz,  covered 
by  the  auto-correlation  function  Rq(nAz,t),  is' equal  to  1/Af,  by  virtue  of 
eqn.  (2-47). 

Because  of  the  attenuation  provided  by  the  envelope  of  the  individual 
pulses,  the  maximum  extent  of  the  target  range  profile  being  processed  is 
assumed  to  be  effectively  limited  to  the  pulse  width,  t .  As  Illustrated  in 
Fig.  2-9c,d  the  extent  of  the  range  profile  auto-correlation  function  for  a 
2-scatterer  target  whose  range  extent  it  r  will  then  be  2t.  In  order  that 
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Fig.  2-9  Illustration  of  th«  Halt  oa  th«  anal  nun  allowable  frequency -step 
size,  At,  required  to  ovoid  olloolng  of  raoge  prof 11a  autocorrelation  fuoeelooo 
obtained  by  a on co he rant  proeaeslng  of  rataraa  f roa  fr«quaney-a tapped  pulse 
tralna. 

(a)  A  two  polat-acattarar  target  having  a  range  extent  of  one  pulaa  width  and 
the  corresponding  coaponent  pulaa  returns. 

(b)  The, high  resolution  range  profile  of  the  target  configuration  of  rig.  2-8a, 
for  a  stapling  tint  half-way  between  the  two  scatterers. 

(c)  Aliased  range  profile  auto-correlation  function  obtained  when  the  frequency 
step  site  of  the  pulse  train  la  too  large  Af>l/(2t). 

(d)  Correct  (unallaaed)  range  profile  auto-correlatlon  function  obtained  when 
the  frequency  step-site  Is  reduced  s*ieh  that  Af<l/ (2t ). 


Che  computed  auto-correlation  function  not  be  corrupted  by  aliasing,  it  is 
necessary  that  the  extent,  NAz,  covered  by  the  auto-correlation  function, 
Rq(z,t)  computed  as  the  transform  (2-45),  be  greater  than  that  of  the 
target's  actual  range  profile  auto-correlation  function,  2t. 

NAz  >  2t 

this  Is  equivalent  to  the  preceding  statement  that: 

Rq(nAz,t)  »  0  for  t  <  nAz  <  y  Az 

Combining  this  constraint  with  the  definition  of  Az  given  by  (2-44)  gives 
the  following  relationship  for  the  constraint  on  the  frequency  spacing  Af. 

Af  <  27  (2-49) 

Note  this  constraint  on  Af  is  the  same  as  that  imposed  earlier  in  the  case 
of  coherent  processing,  although  the  reason  for  the  constraint  in  this  case 
is  somewhat  different.  In  the  present  case,  the  constraint  (2-49)  on  the 
maximum  allowable  value  of  the  frequency  Increment,  Af,  is  Imposed  by  the 
need  to  avoid  corruption  of  the  computed  range  profile  auto-correlation 
function  bj  aliasing.  This  corresponds  to  the  Nyqulst  criterion  that  sta¬ 
tes  that  the  sample  spacing  (in  this  case  in  frequency)  of  a  sampled  func¬ 
tion  (the  pulse  returns)  must  be  less  than  one-half  the  period  of  the 
highest  "frequency"  component  (i.e.  the  pulse  width)  in  order  to  avoid 
aliasing.  In  the  case  of  coherent  processing,  the  corresponding  constraint 
on  Af  required  to  avoid  aliasing  is  I/t,  twice  that  given  by  eqn.  (2-4?). 
However,  in  that  case,  a  further  constaint  on  Af  was  Imposed  in  order  to 
be  able  to  unambiguously  Interpret  the  computed  range  profile.  This  addi¬ 
tional  constraint  is  the  same  as  that  given  by  eqn.  (2-49)  above.  Thus 
eqn  (2*?46)  applies  in  both  eases. 


2.3. 1.1  Effect  r.t  Non-uniform  Frequency  Spacing  for 
Noncoherent  Processing 

The  effect  of  non-uniform  frequency-spacing  on  noncoherent  processing 
using  the  DFT  is  considerably  less  severe  than  in  the  case  of  coherent  pro¬ 
cessing,  previously  described  in  section  2.2. 1.2.  Again,  let  the  actual 
frequencies  transmitted  be  given  by: 

Sk  "  fk  +  ek  (2-21) 

where  fk  are  the  ideal  values  of  the  transmitted  frequencies  corresponding 
to  uniform  spacing  (eqn.  (2-13))  and  ek  are  the  resulting  deviations  about 
the  fk.  The  analysis  of  section  2.2* 1.2  applies  in  this  case  as  well. 
Consequently,  it  follows  that  if  the  bound  on  given  by  eqn.  (2-24)  is 
satistfled:  i.e. 

|  ek  |  <  .01  Af  (2-24) 

then  eqn.  (2-26)  follows: 

V(gk,t)  -  V(fktt)e'j2irekt  (2-26a) 

and  therefore: 

|  V(gkft)  |  2  - j  V(fkft)  I  2  (2-50) 

Therefore,  if  the  magnitude  of  the  frequency  deviations  is  bounded  by 
eqn  (2-24)  or  equivalently  by  (2-25),  then  eqn.  (2-50)  follows,  so  that  the 
range  profile  auto-correlation  function  Rq(s,t)  may  be  computed  from  the 
amplitudes  of  the  returns  j  V(gk}t)  (measured  for  the  actual  transmit  fre¬ 
quencies  gk.  Consequently,  for  noncoherent  processing,  the  applicable  bound 
on  the  allowable  frequency  deviation  is  given  by  eqn.  (2-24) ,  dr  equiva¬ 
lently  by  eqn.  (2-25). 

It  is  interesting  to  compare  this  bound  with  that  derived  previously 
for  the  case  of  coherent  processing,  eqn.  (2-30). 
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eck  |  .1 


.0029  /r 
taax 


(2-30) 


where  tBax  is  the  maximum  range  delay  to  the  target.  The  frequency 
deviation  bound  for  noncoherent  processing  given  by  eqn.  (2-24)  may  also  be 
expressed  by  eqn.  (2-25). 

eikU^  (2-25) 


where  r  is  the  effective  (compressed)  pulse  width  of  the  radar.  A  com¬ 
parison  of  these  two  bounds  gives: 


eik  1.7  cmax 

-TT  “T 


(2-51) 


Since  tgjax  is  usually  eqoal  to  the  radar  pri,  and  N  is  usually  less  than 
128,  the  value  of  the  bound  on  allowable  frequency  deviation  for  non¬ 
coherent  processing,  |  j  ,  will  be  larger  than  that  for  coherent  pro¬ 
cessing,  |  eck  j  ,  by  about  one-tenth  the  ratio  of  the  radar  pri,  Uax  to 
the  (compressed)  radar  pulse  width  t.  In  most  cases  this  ratio  will  be  on 
the  order  of  10^  or  more.  Consequently,  the  frequency  error  tolerances  on 
the  deviations  of  the  pulse  train  frequencies  relative  to  the  values  for 
uniform  frequency  spacing  for  noncoherent  processing  to  obtain  the  range 
profile  auto-correlation  function  ire  considerably  more  relaxed  than  those 
required  for  coherent  processing  to  obtain  the  range  profile. 
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EFFECT  OF  COHERENT  AND  NONCOHERENT  PROCESSING  OF  RETURNS  FROM 
FREQUENCY-AGILE  RADAR  PULSE  TRAINS  ON  TARGET-CLUTTER  CONTRAST  RATIO 
FOR  TARGETS  IMBEDDED  IN  GROUND  CLUTTER  ■ 


A  common  application  of  frequency-agile  radar  pulse  trains  is  the 
detection  and  imaging  of  fixed  targets  located  in  a  ground  clutter 
background.  A  typical  scenario  for  this  application  is  illustrated  in 
Fig.  3-1.  Fixed  target  detection  for  this  scenario  is  based  primarily  upon 
the  contrast  in  amplitude  between  the  return  from  a  target  and  that  from 
the  surrounding  clutter,  based  on  the  fact  that  returns  from  man-made 
metallic  targets  are,  on  the  average,  stronger  than  those  from  clutter 
within  the  same  size  resolution  cell.  Specifically,  the  returns  from  a 
contiguous  patch  of  resolution  cells,  such  as  the  rectangular  array 
illustrated  in  Fig.  3-2,  are  examined,  and  the  amplitude  of  the  return 
from  the  center  cell  of  the  array  is  compared  against  the  average  of  the 
returns  from  the  surrounding  cells.  If  there  is  a  target  in  the  center 
cell  of  the  array,  and  the  remaining  cells  of  die  array  contain  only 
clutter,  a  measure  of  the  contrast  between  the  target  and  the  clutter  is 
the  average  target-to-clutter  power  return  ratio,  T/C.  Strictly  speaking,  , 
If  the  size  of  the  resolution  cells  is  much  larger  than  the  physical  extent 
of  the  target,  the  cell  containing  the  target  will  also  contain  some 
clutter,  so  that  the  contrast  ratio  between  the  tacger.  cell  and  the  clutter, 
cells  will  actually  be  (T+C)/C,  rather  than  just  T/C.  In  any  case,  since 
the  clutter  return  from  a  given  resolution  cell  is  proportional  to  the  pro¬ 
jected  ground  area  of  that  cell,  the  value  of  T/C  will  increase  with 
decreasing  cell  size  (i.e. ,  increasing  resolution),  as  long  as  the  cell 
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size  is  greater  than  the  physical  extent  of  the  target.  This  effect  is 
illustrated  in  Fig.  3-3. 

When  the  resolution  of  the  radar  increases  to  the  point  where  the 
resolution  cell  size  becomes  smaller  than  the  physical  extent  of  the 
target,  the  individual  scattering  centers  of  the  target  start  to  become 
resolved.  In  that  case,  the  value  of  T/C  continues  to  increase,  but  at  a 
less  rapid  rate. 

As  outlined  above,  target-clutter  signal  ratio,  T/C,  is  usually  a  good 
measure  of  the  amplitude  contrast  between  returns  from  targets  and  clutter. 
Since  T/C  is  inversely  proportional  to  the  size  of  the  resolution  cell, 
decreasing  the  size  of  the  resolution  cell  is  one  meads  of  enhancing  the . 
target-clutter  contrast  required  for  reliable  detection. 

The  purpose  of  this  chapter  is  to  determine  the  effects  of  both 
coherent  and  noncoherent  processing  of  the  returns  from  frequency-agile 
pulse  trains  on  the  target-clutter  contrast  in  the  resulting  high- 
resolution  range  profile  and  range  profile  auto-correlation  functions.  The 
effect  of  this  processing  on  target-to-clutter  contrsst  is  represented  in 
terms  of  a  parameter  called  the  contrast  enhancement  ratio.  Eg.  This  para¬ 
meter  is  defined  as  the  ratio  between  the  average  target-clutter  contrast 
in  the  resulting  high  resolution  range  profile  or  range-profile  auto¬ 
correlation  function,  and  the  average  unresolved  target-clutter  contrast 
before  processing  (i.e»  the  average  contrast  between  the  target  and  clutter 
for  the  unresolved  slngle^pulss  returns,  T/C).  In  general,,  this  contrast 
enhancement  ratio  is  an  increasing  function  of  the  resolution  enhancement, 
I,  provided  by  the  pulse  train  bandwidth.  However,  the  contrast  enhan¬ 
cement  ratio  Ex  is  also  a  function  of  the  number,  M,  of  disctete  target 
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TARGET  CLUTTER  SIGNAL  RATIO  =  13.1-0  *  13  1  08 

Fig.  3-3  Effect  of  resolution  cell  site  on  target-clutter  signal  ratio  for 
a  fixed  ground  target. 


scatterers  that  become  resolved  in  the  high  resolution  profile,  and  the 
unresolved  target-to-clutter  signal  ratio,  T/C. 

The  effect  of  coherent  processing  on  the  target-clutter  contrast 
enhancement  provided  by  the  high-resolution  range  profile  is  discussed  in 
section  3.1.  Since  coherent  processing  effectively  divides  the  original 
single-pulse  range  resolution  cell  into  finer  range  subcells,  the  contrast 
between  the  average  target  amplitude  and  the  average  clutter  amplitude  in 
the  resulting  high  resolution  range  profile  increases  nearly  in  proportion 
to  the  corresponding  enhancement  in  resolution. 

The  effect  of  noncoherent  processing  on  the  contrast  enhancement  pro¬ 
vided  by  incoherent,  processing  is  the  subject  of  section  3.2.  In  this 
case,  the  enhancement  in  contrast  between  the  amplitude  of  the  target  peaks 
in  the  range  profile  auto-correlation  function  and  the  clutter  background 
level  of  that  function  is  much  less  than  that  obtained  for  coherent  pro¬ 
cessing.  Indeed,  the  target-clutter  contrast  of  the  range  profile  auto¬ 
correlation  funcrlon  is  often  less  than  the  original  unresolved  target- 
clutter  ratio,  T/C.  In  other  words,  it  is  often  more  difficult  to  detect 
the  presence  of  a  target  by  looking  for  peaks  in  the  high-resolution  range 
profile  auto-correlation  function  than  by  comparing  the  average  amplitudes 
of  the  unresolved  target  and  clutter  returns. 

3. 1  ,  Effect  of  Coherent  Processing  on  Target-to-Clutter  Contrast  Ratio 
As  previously  discussed  in  section  2.2.1,  coherent  processing  of  the 
returns  from  a  frequency-agile  pulse  train  gives  the  true  high-resolution 
range  profile  of  a  target,  with  a  resolution  equal  to  about  Ar  -  c/2B, 
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where  B  “  NAf,  the  pulse  Crain  bandwidth.  Suppose  Che  single  pulse  band¬ 
width  of  the  radar  is  W  ■  1/t,  where  t  is  the  compressed  pulse  width.  The 
resulting  resolution  improvement  can  be  expressed  as  the  ratio  of  the 
single  pulse  range  resolution  AR(*cr/2),  to  that  of  the  high  resolution 
range  profile,  Ar.  Equivalently,  this  resolution  improvement  factor  can 
also  be  expressed  as  the  ratio  of  the  pulse  train  bandwidth  to  that  of  a 
single  pulse. 

I  "  AR/Ar  ■  NAft  or  ~— 

"  (3-1) 

The  parameter  I  is  called  the  resolution  improvement  factor. 

The  effect  of  this  improvement  in  resolution  on  target-to-clutter 
ratio  Is  as  follows.  In  the  simplest  case  where  all  of  the  target  scat¬ 
tering  centers  are  clustered  sufficiently  close  together  so  as  to  remain 
unresolved  at  the  higher  resolution  Ar,  for  example  as  in  the  case  of  a 
single  scatterer  such  as  a  corner  reflector,  the  Improvement  In  target-to- 
clutter  ratio  provided  Is  exactly  equal  to  the  resolution  improvement ,  I. 
The  reason  fcr  this  follows:  in  general,  the  target-to-clutter  ratio  is 


given  by: 


(3-2) 


where  ot  is  the  target  cross*^sectlon,  oQ  is  the  clutter  reflectivity  per 
unit  area,  and  Ac  is  Che  projected  area  of  the  resolution  cell  of  length 
Ah/ cos  <p  and  width  RA8. 

^  •  RAOAR/cos  $  (3-3) 

where  R  is  the  range,  A8  the  effective  beamwldth  of  the  radar,  AR  the  range 
resolution,  and  +  the  depression  angle.  Therefore,  for  constant  values 


of  all  parameters  other  than  the  range  resolution  AR,  T/C  varies  inversely 


with  A R  as  follows: 


,  „  k  at  cos  A 

‘  h  ■■  K  ■  -TiT-1 


(3-4) 


Hence  the  improvement  in  T/C  provided  by  Increased  range  resolution  when 
the  target  remains  unresolved  (i.e.  when  ot  is  unaffected  by  the  change  in 
resolution)  is: 


T/C(Ar)  AR 
T/C(AR)  “  Ar  “  L 


(3-5) 


In  most  practical  cases,  the  high-resolution  subcell  containing  the 
target  will  also  contain  some  clutter  because  of  the  finite  extent  of  that 
cell  in  the  cross-range  dimension.  Consequently,  the  observed  target 
return  will  actually  consist  of  target  plus  clutter.  Therefore,  the 
measure  of  the  target  return  enhancement  provided  by  high  resolution  should 
be. based  on  the  observed  target-plus-clutter  to  clutter  ratio,  and  thus  the 
contrast  enhancement  ratio  for  coherent  processing  is  defined  as  follows: 


_  T/C(Ar)  +  1  I»T/C(AR)  +  1 
1  T/C(AR)  +  1  T/C(AR)  +  1 


(3-6) 


Hence,  in  general,  is  seen  to  be  a  function  of  both  the  resolution 
Improvement  ratio,  I,  and  the  unresolved  target  to, clutter  ratio  T/C(AR). 


For  large  values  of  T/C,  the  above  simplifies  to: 

E1  «  I  for  T/C(AR)  >  4 


(3-6a) 


Thus  at  high  values  of  T/C,  the  contrast  enhancement  ratio  for  coherent 
processing  becomes  independent  of  T/C.  When  the  range-resolution,  Ar, 
obtained  by  coherent  processing  is  fine  enough  to  resolve  the  Individual 
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scatterers  of  the  target,  the  effect  of  Increasing  resolution  on  Ej  is  a 
bit  more  complicated,  but  still  straight  forward. 

Consider  the  case  where  the  target  is  modelled  by  M  distinct  point 
scatterers.  Figure  3-4a  describes  the  range  profile  of  a  target  composed 
of  such  point  scatterers  in  clutter,  within  the  extent  of  the  radar's 
single  pulse  range  resolution,  AR  -  i.e.  the  area  of  the  radar's  single¬ 
pulse  range  resolution  cell  is  assumed  to  be  large  enough  to  encompass  not 
only  the  entire  target  but  also  some  of  the  surrounding  clutter.  Let  the 
range  separations  between  these  M  scatterers  be  given  by  djj,  i»l,...M, 
j  *  1,  ....,  M.  In  general,  the  complex  signal  voltage,  s^,  of  the  unre¬ 
solved  return  from  these  M  scatterers  at  frequency  f^  will  be  given  by 


M  j  i*fk<*lj 

Sk  -  l  Aj  e  c  (3-7) 

j-1 

Air 

where  the  term  — fkdij  represents  the  relative  phase  of  the  return  from  the 
jth  target  scatter  with  respect  to  that  of  the  first  target  scatterer,  and 
Aj  is  the  complex  voltage  of  the  return  from  the  jtl1  scatterer.  The 
received  power  at  frequency  f^  is  therefore  given  by: 


*  M  2  MM  . 

fk  m  sJcsk  *  l  I  A j  |  +  2  l  £  AjAj  j  cos(  — •  f kdij ,  +(frij )  (3-8) 

j-1  j-1  i-j+l  c  j  j 

where  use  has  been  made  of  the  fact  that:  djj  S  dj^-dy 

,  la  [A4A4*] 

*”d  "h"8  ns 5 1“'1 


As  indicated  by  the  second  term  of  eqn.  (3-8),  because  of  interference 
between  the  returns  from  the  individual  scatterers,  the  received  power: 


65 


RANGE  DELAY  TIME 


SAMPLING  TIME  (Range  Delay) 


Fig.  3-4  Formulation  of  the  high  resolution  range  profile  of  a  3-point  scat* 
terer  target  in  clutter.  Scatterer  range  separation  is  less  than  one  pulse 
width. 

(a)  Actual  range  profile  showing  the  overlapping  component  pulse  returns  from 
the  three  point-scatterers. 

(b)  Profile  of  total  power  returns  vs.  A/D  sampling  time  of  the  target-clutter 
profile  of  Fig.  3-4a  for  three  different  illumination  frequencies. 

(c)  The  high  resolution  range  profile  obtained  as  the  DFT  of  the  complex  samples 
of  the  returns  from  a  frequency-stepped  pulse  train  received  from  the  target- 
clutter  profile  of  Fig.  3-4a. 


Pk-|sk|2 

will  vary  as  a  function  of  radar  frequency  fk.  However,  the  average  value 
of  P^,  averaged  over  many  different  frequencies,  will  be  proportional  to 
the  sum  of  the  cross-sections  of  the  individual  scatterers,  i.e. 

r  II  s*|2  5  Il*j|  -  «I«j  <3-») 

fc-1  J“1  j"1 

The  radar  cross-section  <jj  of  each  of  the  j  scatterers  is  proportional  to 
the  reflected  power,  hence  cj  ■  -i|  Aj  |2.  This  situation  is  illustrated  in 
Fig.  3-4b  for  single-pulse  returns  from  the  target-in-clutter  profile  of 
Fig.  3-4a  as  a  function  of  sampling  time.  At  any  given  value  of  sampling 
time  the  received  power  return  from  the  given  resolution  cell  (Fig.  3-4a) 
varies  from  pulse  to  pulse  with  changing  frequency,  as  shown  in  Fig.  3-4b, 
because  of  the  interference  between  the  returns  from  the  individual  scat- 
tererlng  points  (of  the  target  and  clutter)  in  that  cell.'  However,  the 
average  power  return  over  frequency  is  proportional  to  the  combined  pulse- 
shape  weighted  cross-sections  of  the  three  scatterers  plus  the  clutter, 
as  Illustrated  by  the  dashed  curve  in  Fig.  3-4b.  In  this  case  the  target- 
to-clutter  ratio  is  defined  as  the  ratio  of  the  mean  target  return  over 
frequency  to  the  mean  clutter  return  over  frequency,  and  as  in  the  single 
target  scatterer  case,  T/C  is  given  by  eqn.  (3-2),  where  the  total  target 
cross-section  is  the  sum  of  the  individual  scatterer  cross-sections: 

M 

Ot  •  l  (3—10) 

J-l 

Now  suppose  that  the  frequency-agile  pulse  returns  from  the  resolution 
cell  of  Fig.  (3-4a)  are  coherently  processed,  and  the  resulting  resolution. 
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Ar,  is  fine  enough  to  resolve  the  individual  target  scatterers,  cr j.  For 

simplicity,  assume  that  the  M  individual  target  scatterers  are  all  of  equal 

at 

magnitude  aj;  i.e.  oj  ■  —  .  The  resulting  high-resolution  range  profile 
is  illustrated  in  Fig.  (3-4c).  In  effect,  what  has  happened  is  that  the 
extent,  AR,  of  the  single  pulse  range  resolution  cell,  has  now  been  sub¬ 
divided  into  I  smaller  range  cells  each  having  a  finer  resolution  of 
It  ■  AS/1.  Therefore,  as  before,  the  clutter  return  in  each  of  the  fine 
resolution  cells  is  reduced  by  the  resolution  Improvement  factor.  In  other 
words,  if  the  clutter  return  from  the  original  single-pulse  range  resolu¬ 
tion  cell  was  C,  the  clutter  return  in  each  of  the  fine  resolution  cells  is 
now  C/I.  Hence  the  target-clutter  ratio  of  the  returns  from  each  of  the 
resolved  scattering  centers  now  becomes: 

T/C(Ar)  -  i  |  T/C( AR)  (3-11) 

Hence  when  the  target  is  composed  of  M  discrete  scatterers,  all  of  which 
are  of  roughly  equal  amplitude,  and  the  resolution  becomes  fine  enough  to 
individually  resolve  all  of  these  scatterers,  the  contrast  enhancement 


ratio  for  coherent  processing  becomes; 

_  I/M«T/C(AR)+1 

T/C(AR)+1  U 

Note  that  eqn.  (3-12)  reduces  to  (3-6)  for  the  single  resolved-scatterer 


(3-12) 


case  (i.e.,  when  M»l).  Therefore,  the  contrast  enhancement  is  reduced  over 


the  single  scatterer  case  roughly  in  proportion  to  the  number,  H,  of  scat¬ 
terers  that  become  resolved.  This  target-clutter  contrast  enhancement  fac¬ 
tor,  &£,  is  plotted  as  a  function  of  resolution  Improvement  factor,  I,  in. 
Pig.  3-5  for  the  asymptotic  case  of  large  T/C.  At  lower  values  of  T/C,  the 
dependence  of  Ex  on  I  is  similar,  except  that  the  values  of  &x  will  be 


somewhat  smaller  in  accordance  with  eqn.  (3-12).  As  the  resolution  Im¬ 
provement  Increases,  so  does  the  number,  M,  of  scattering  centers  that 
become  resolved.  Hence,  the  overall  contrast  enhancement  ratio  curve 
shifts  towards  the  asymptotic  curves  for  higher  values  of  M  as  I  Increases. 
The  Ideations  of  the  exact  transitions  between  the  curves  for  different  M 
values  will  depend  upon  the  distribution  of  the  locations  of  the  M  scat¬ 
tered  within  the  single-pulse  range  resolution  cell,  the  transitions  bet¬ 
ween  the  curves  for  different  values  of  M  will  not  start  to  tafce  place 
until  I  increases  to  the  point  where  Ar  Is  reduced  to  a  value  equal  to  the 
maximum  scatterer  separation.  Hence,  If  the  scattered  are  widely 
separated  in  range  across  the  cell,  these  transitions  will  start  to  occur 
at  low  values  of  I,  and  conversely. 

It  should  be  noted  that  the  contrast  Improvement  ratios  given. by  eqns. 
(3-6)  and  (3-12)  are  premised  on  the  assumption  that  the  high  resolution 
range  profile  Is  not  corrupted  by  aliasing  -  l.e.  that  no  extraneous 
clutter  returns  are  folded  into  the  profile.  In  other  words,  the  frequency 
step  size  Af  must  be  such  that  Af  1/2**  This  constraint  will  only  be  met 
if  the  number,  N,  of  pulses  in  the  train  Is  such  that  i!l>  21.  In  other 
words,  In  ordtr  eqns  (3-6)  and  (3-10)  to  hold,  the  number  of  pulses  In  the 
train  must  be  at  least  twice  as  great  as  the  Indicated  resolution  improve¬ 
ment  factor  I. 

In  summary,  coherent  processing  of  frequency-agile  returns  to  obtain 
high  resolution  range  profiles  will  Increase  the  contrast  of  the  target 
returns  relative  tc  the  clutter  approximately  in  proportion  to  the  resolu¬ 
tion  Improvement  factor  I,  provided  that  the  frequency  spacing,  Af  1/2* ,  so 
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as  to  prevent  folding  of  extraneous  clutter  into  the  profile.  When  the 
resolution  becomes  fine  enough  to  resolve  the  individual  scattering  centers 
of  the  target,  the  contrast  of  these  scattering  centers  varies  approxima¬ 
tely  in  proportion  to  1/h  where  M  is  the  number  of  resolved  scatterers. 
Hence,  as  resolution  increases  to  the  point  where  the  Individual  target 
scatterers  become  resolved,  target  contrast  continues  to  increase  with  I, 
but  at  a  somewhat  reduced  rate  tmtil  all  of  the  major  scatterers  of  the 
target  have  been  resolved,  as  illustrated  in  Fig.  3-5.  Once  all  of  the 
dominant  scatterers  have  become  resolved,  contrast  enhancement  again 
increases  in  direct  proportion  to  I. 

3.2  Effect  of  incoherent  Processing  on  Target-Clutter  Contrast  Ratio 

As  was  previously  discussed  in  section  2.3,  when  coherent  processing 

of  the  returns  from  a  frequency-agile  pulse  train  is  either  not  possible  or 

not  feasible,  it  is  still  possible  to  obtain  the  auto-correlation  function 

of  the  high  resolution  range  profile,  rather  than  the  range  profile  Itself. 

The  effective  resolution  of  the , resulting  auto-correlation  function  is  the 

c  1 

same  as  that  of  the  associated  range  profile,  Ar,  *  y  .  Hence,  the 

resolution  improvement,  I,  defined  by  eqn  (3-1)  is  the  same  in  both  cases. 

Since  the  resolution  improvement  obtained  from  incoherent  processing 
is  the  same  as  that  obtained  from  coherent  processing,  a  reasonable 
question  is  whether  the  target-clutter  contrast  ratio  of  the  target  "peaks” 
of  the  high-resolution  range-profile  auto-correlation  function  also 
Increases  with  increasing  resolution  in  the  same  way  as  in  the  case  of 
coherent  processing.  The  answer  is  no. 


In  general,  Che  target-clutter  contrast  ratio  between  the  target  peaks 
and  the  clutter  background  level  of  the  range  profile  auto-correlation 
function  will  be  less  than  tnat  for  the  original  unresolved  target  -  i.e. , 
less  chan  T/'C.  However,  this  contrast  ratio  does  increase  with  increasing 
resolution,  but  at  a  such  slower  rate  than  in  the  case  of  coherent 
processing.  As  a  result,  noncoherent  processing  of  returns  froa  frequency- 
agile  pulse  trains  is  not  nearly  as  effective  in  enhancing  the  detec¬ 
tability  of  fixed  targets'  in  ground  clutter  as  coherent  processing  is.  The 
basic  reason  for  this  soaewhat  surprising  result  is  that  the  "clutter** 
background  level  of  the  range-profile  auto-correlation  function  is  propor¬ 
tional  to  the  product  of  the  aaplitude  of  the  target  scatterera  with  that 
of  the  clutter  in  the  resolution  cell  -  hence,  as  the  amplitude  of  the 
individual  target  scatterera  increases,  so  does  that  of  the  clutter 
background  level.  A  acre  quantitative  analysis  of  this  effect  is  given 
below. 

Again,  consider  the  target-ln-clutter  aodel  of  Fig.  3-4a,  and  the 
corresponding  high  resolution  range  profile  of  Fig.  3-4c.  As  previously 
discussed  in  Section  2.3,  the  function  Rq(x)  obtained  by  noncoherent  pro¬ 
cessing  of  the  returns  froa  a  frequency-agile  pulse  train'  is  the  auto¬ 
correlation  of  the high-resolution  range  profile,  q(x),  illustrated  in  Fig. 
3-4c. 

In  general,  obtaining  an  analytical  expression  for  Rq(x)  for  the  case 
of  a  distributed  target  in  clutter  is  difficult.  However,  a  discrete 
approxlaatlon  to  the  auto-correlation  function,  Rq(x),  of  the  range  profile 
illustrated  in  Fig.  3- 4c  can  be  relatively  easily  evaluated,  subject  to  the 
following  assuaptlons. 


1.  The  target  model  consists  of  M  discrete  scatterers  distributed 

in  range  and  having  equal  cross-sections,  a t/M',  but  arbitrary  phases. 

2.  The  separations  between  the  target  scatterers  are  non-equal  and 
unique  for  all  pair  combinations  of  the  M  scatterers. 

3.  The  distribution  of  the  clutter  within  the  single  pulse  range 
resolution  ceil  is  statistically  uniform. 

4.  The  single-pulse  range  resolution  cell  AR  is  subdivided  into  I 
equal  subcells  of  range  extent,  Ar  -  AR/I,  and  the  clutter  return 
from  each  of  these  subcells  is  a  complex  Gaussian  process  that  is 
independent  from  one  subcell  to  the  next.  However,  the 
statistics  of  the  subcell  clutter  distribution  are  the  same  for 
all  subcells.  Note  that  since  each  of  the  subcells  are  implicitly 
resolved,  1  is  equal  to  the  resolution  improvement  ratio,  con¬ 
sistent  with  the  earlier  definition  of  I. 

5.  Those  subcells  containing  the  discrete  target  scatterers  also  con¬ 
tain  clutter., 

6.  The  effect  of  the  relative  attenuation  of  the  returns  from  dif¬ 
ferent  subcells  due  to  the  shape  of  the  radar  pulse  j.s  neglected 
-  In  effect  a  rectangular  pulse  shape  is  assumed. 

7.  The  resolution  Improvement  ratio,  I,  Is  much  greater  than  the  effec¬ 
tive  number  of  discrete  target  scatterers,  M. 

8.  Boundary  effects  In  the  computation  of  the  correlation  function 
are  neglected.  The  maximum  separation  between  target;  scatterers 
is  assumed  small  relative  to  the  range  extent,  AR,  of  the  cell,  so 


that  the  target  profile  Is  not  affected  by  aliasing.  However,  the 
range  auto-correlation  of  the  subcell  clutter  returns  is  assumed 
to  be  circular. 

The  discrete  approximation  to  the  auto-correlation  of  the  range  pro¬ 
file  illustrated  in  Fig.  3-4c  is  then  computed  as  follows.  The  high 
resolution  range  profile  of  Fig.  3-4c  is  subdivided  into  I  discrete  range 
bins,  each  of  width  Ar*  The  return  from  each  bin  is  complex  and  is  due 
either  to  .  'iter,  or  to  one  of  the  M  point  scatterers  of  the  target  plus 
clutter.  In  accordance  with  assumption  4,  the  clutter  returns  from  each 
oubcell  (i.e.  bin)  are  independently  distributed,  with  an  exponential  power 
distribution  that  is  the  same  for  all  1  range  subcells.  The  mean  value  of 
the  power  return  from  the  clutter  in  each  range  subcell  is  thus  C/1,  where 
C  is  the  average  clutter  power  returned  from  the  entire  single  pulse  reso¬ 
lution  cell,  AR.  The  power  return  from  each  of  the  M  target  point  scat¬ 
terers  is  Ot/M,  where  ot  is  the  total  target  cross-section.  In  accordance 
with  assumptions  1  and  4  above,  the  phase  of  the  complex  return  from  each 
of  the  I  .range  subcells,  whether  due  to  the  clutter  or  target  plus  clutter, 
is  random  and  uniformly  distributed  over  the  interval  -w,  +ir. 

Let  the. complex  voltage  return  from  the  clutter  in  any  subcell  be 
represented  by: 

qi  -  at  e^1  i-  I,. ...I  (3-13) 

the  variable  a^  is  Rayleigh  distributed  and  the  variable  ^1  is  uniformly 
distributed  over  -w ,  +r .  as  mentioned  above,  for  all  subcells.  In  the  case 
of  a  target  subcell,  consists  of  the  coherent  sum  of  two  terms;  one  for 


the  target  scatterer  and  the  other  for  the  clutter  In  that  subcell,  as 
follows: 

target  return  clutter  return 
qi  *  ~  e^0t  +  a*  e-i$i  (3-13a) 

The  amplitude  of  the  target  term  is  a  constant  equal  to  o^/M.  The  phase 
of  the  target  term  is  again  uniformly  distributed  over  ~r,  +ir. 

Hence,  in  general,  the  expression  for  nay  be  written  as:' 

qi  -  at  eJ*i  +~  eJ9m  6(i-ja)  (3-13b) 

where  jm  is  the  subcell  location  of  the  mth  target  scatterer  and 

1  i-Jm 

SU-Jm)  -  { 

0  i*ja 

The  discrete  complex  auto-correlation  function  of  the  M-point-scatterer 
target-in-clutter  range  profile  described  by  eqn.  (3-13b)  is  therefore: 

I  * 

RqOO  ■  l  qi  qi+k  (3-i4) 

i-l 

where  is  given  by  eqn.  (3-1 3b)  and  qi-Hc  has  been  circularly  extended: 
qj-qj_I  for  J>1,  in  accordance  with  assumption  7. 

The  auto-correlation  function  described  by  eqn.  (3-14)  is  a  random 
function,  because  of  the  random  nature  of  clutter  returns  used  in  the 
model.  Furthermore,  Rq(k)  is  in  general  complex,  except  for  k-0  where  it 
Is  real  and  positive.  However,  as  stated  previously  in  Chapter  2,  we  are 
primarily  Interested  in  the  magnitude  of  this  profile,  |Rq(k)|.  The 
statistics  of  the  range  auto-correlation  magnitude  profile,  jRq(k)  |  are 
derived  in  Appendix  A.  The  distribution  of  |Rq(k)|  takes  on  three  dif¬ 
ferent  forms,  as  a  function  of  k,  as  described  in  Table  3-1. 
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1. 


2. 


3. 


Note 


Note 


k*0,k* jm 


Table  3-1 

StaClstlcs  of  the  Range  Auto-correlation 
Magnitude  Profile 


Distribution 
Value  of  k  of  |Rq(k) j 


approx. 

normal 


approx. 

Rayleigh 


Mean 

of  lRq(k)| 
qt+C 

/*  l+2ot/C 

'  T  T 


Std.  Deviation 
of  jRq(k)| 


yi+2at/C 


*see  Note 
1  below 


/(4-ir)  l+2qt/C 


k/O.k^Ju, 


approx. 

Ricean 


/l+2at/C  (at/C)2 
/  I  +  M2 


**see  Note  2 
below 


1:  The  result  for  the  standard  deviation  of  Rq(0)  cited  here  is  based 
upon  the  discrete  target-plus-clutter  model  described  earlier  in 
which  it  is  implicitly  assumed  that  the  spacing  between  target 
point  scatterers  is  limited  to  the  integral  multiples  of  Ar,  kAr, 
where  Ar*AR/I  is  the  range  resolution  corresponding  to  the  pulse 
train  bandwidth  and  k<I  is  an  Integer.  It  can  be  shown  that  if 
this  assumption  is  relaxed,  so  that  the  separations  between  the 
target  point  scatterers  can  take  on  a  continuum  of  values  less 
than  AR,  the  standard  deviation  of  Rq(0)  (i.e.,  case  1)  Increases 


to:  C 


/l  +  2(< 


®t/C)+  — additional  term  of 


Ji 


(M-l) 

M 


(ot/ C): 


jr 


represents  the  contribution  of  the  scintillation  between  the 
returns  from  different  target  scatterers  with  changing  frequency. 
This  additional  term  vanishes  either  in  the  case  of  a  single  scat- 
terer  or  where  all  of  the  scatterer-pair  separations  are.  Integral 
multiples,  of  Ar.  Thus,  in  those  cases  the  standard  deviation  of 
Rq(0)  reduces  to  the  value  given  in  Table  3-1.  The  mean  value  of 
Rq(0)  cited  in  Table  3-1  is  not  affected  by  these  considerations. 
However,  it  is  conjectured  that  the  mean  and  standard  deviation  of 
|Rq(k)|  cited  in  Table  3-1  for  cases  2  and  3  will  also  change  in 
proportion  to  the  standard  deviation  of  Rq(0)  When  M>1. 

2  The  mean  value  cited  for  the  Ricean  distribution  of  the  target 
peaks  (case  3)  is  an  approximation  (see  App.  A).  Note  that  it 


reduces  to  within  a  factor  of  v—  ■  1,12  of 


the  mean  value  of  the 
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Rayleigh  distribution  (case  2)  when  at/C“0.  The  error  in  this 
approximation  decreases  as  at/C  increases.  The  standard  deviation 
of  |Rq(k)|  for  case  3  varies  between  a  lower  bound  equal  to  the 

(l+2ot/C)  " 

standard  deviation  for  case  2  (i.e.,  *  /*- - - j- -  )  when 

(<jf/c)  +  0,  ai*d  an  upper  bound  of  B^l.53  as  (ot/C)  -► 


The  range  profile  auto-correlation  function  RqOO  of  a  typical  multi¬ 
point  scatterer  target  in  clutter,  corresponding  to  the  high  resolution 
range  profile  of  Fig.  3-4c,  is  illustrated  in  Fig.  3-6.  It  consists  of  a 
main  peak  at  k*0,  that  rapidly  drops  to  a  background  pedestal  for  k*0,  and 
a  number  of  target  peaks  that  rise  above  this  pedestal  at  values  of  k 
corresponding  to  the  individual  scatterers  pair  separations.  There, Is  one 
target  peak  for  each  unique  scatterer-pair  separation  located  at  k“|jm-jn|, 
where  jm  and  jn  are  the  fine-range  indices  of  the  two  scatterers  comprising 
the  given  scatterer  pair.  Hence  if  the  target  is  composed  of  M  discrete 
scatterers  there  will  be  M(M-l}/2  unique  scatterer-pair  combinations  and 
the  total  number  of  target  peaks  above  the  background  predicted  will  be: 

Np<M(M-l)/2  (3-15) 

The  equality  in  (3-15)  holds  only  when  the  M(M-l)/2  scatterer  pair  separa¬ 
tions  are  all  inique  within  an  amount  equal  to  at  least .Ar,  the  resolution 
of  the  range  profile  auto-correlation  function. 

In  the  case  of  Fig.  3-6,  the  three  target  point  scatterers  are 
labelled  A, 8,  and  C,.as  indicated  in  the  high  resolution  range  profile  of 
Fig.  3-6a.  Hence  there  arc  also  3  peaks  (  Np"  M(M-i)/2  •  (3)(2)/2  )  in  the 
corresponding  range  profile  auto-correlation  function  of  Fig.  3-6b.  These 
peaks  are  labelled  AB,  AC,  and  BC,  corresponding  to  the'  scatterer  pairs  of 
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Fig.  3-6a  with  which  they  are  associated.  Note  that  the  distance  of  each 
peak  of  the  auto-correlation  function  from  the  origin  in  Fig.  3-6b,  is 
equal  to  the  range  delay  separation  between  the  corresponding  pairs  of 
scatterers  in  Fig.  3-6a.  Thus  the  location  of  peak  AC  in  Fig.  3-6b  is 
equal  to  the  range  delay  separation  between  scatterers  A  and  C  in  Fig.  3-6a. 

The  statistics  of  the  central  peak  at  k*0,  the  level  of  the  clutter 
background  pedestal,  and  the  so-called  “target  peaks"  of  Rq(k)  were  pre¬ 
viously  given  in  Table  3-1.  The  issue  at  hand  is  how  good  these  target 
peaks  in  the  range  profile  auto-correlation  function  are  as  an  indicator  of 
the  presence  of  a  target  in  the  corresponding  range  cell.  In  particular, 
how  effective  are  such  peaks  as  an  indicator  of  target  presence  relative  to 
using  the  ratio  of  the  average  amplitudes  of  the  target-plus-clutter  and 
clutter  only  returns  (obtained  from  different  resolution  cells)  as  a 
measure  of  target-to-clutter  contrast  -  i.e.,  relative  to  (at/C  +  1)? 

Note  that  the  contrast  between  the  target  peaks  and  the  clutter 
background  level  of  Rq(k)  is  independent  of  the  absolute  values  of  both 
the  target  and  clutter  return.  As  a  result,  normalizing  Rq(k)  such  that 
Rq(0)al,  will  not  affect  this  contrast. 

The  probability  that  the  presence  of  a  target  peak  in  Rq(k)  can  be 
reliably  detected  is  primarily  dependent  upon  the  contrast  in  amplitude 
between  the  mean  value  of  that  peak  and  the  mean  value  of  the  background  of 
Rq(k).  Let  the  mean  value  of  one  of  the  "target  peaks"  in  Rq(k)  for  an  M- 
point  scatterer  target  in  clutter  be  Pc.  Let  the  mean  value  of  the 
corresponding  background  pedestal  level  by  G.  Then  from  Table  3-1: 


This  targe t-peak-to-backg round  contrast  ratio  for  non-coherent  processing 
corresponds  to  the  resolved  contrast  for  coherent  processing  given  by  eqn. 
(3-12a). 

For  high  values  of  target-to-clutter  ratio,  l.e.,  o^/C  >_  2,  and  for 
(I/M^)  >  4,  {this  reduces  to  approximately:' 


£t 

G 


i 

2  M2  C 


(3-19) 


As  before,  the  target-to-clutter  enhancement  factor  is  based  upon  the 
relative  target  visibility  before  and  after  processing  -  on  the  combined 
return  from  target  plus  clutter,  relative  to  that  from  the  clutter  alone. 
Before  processing,  the  relative  target  visibility  is  (ot  +  C)/C,  whereas 
after  processing  it  is  (P^/G)  (the  amplitude  ratio  of  the  target  peaks  to 
the  background  level).  Therefore/  the  contrast  enhancement  ratio 


for  .noncoherent  processing  Is  given  by: 


Pt/C  /I  »  2  (ot/C)  +  (I/M2)(ot/C)2 
®t/C  +  1  "  (ot/C  +  1)  /^7T /I  +  2(ot/C) 


(3-20) 


As  In  the  coherent  processing  case,  eqn.  (3-12),  this  target  contrast 
enhancement  ratio  for  non-coherent  processing  is  a  function  of  the 
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unresolved  target-to-clutter  ratio  (a^/ C),  the  resolution  improvement  fac¬ 
tor  I,  and  the  number  of  resolved  scatterers,  M. 

The  above  expression  for  the  target  contrast  enhancement  ratio  Ej  is 
plotted  for  several  different  values  of  (ot/ C)  as  a  function  of  resolution 
improvement,  I,  for  the  case  of  two  resolved  scatterers  in  Fig.  3-7. 

The  most  significant  feature  of  these  curves  is  that  in  contrast  to 
the  coherent  processing  case,  the  target-clutter  contrast  enhancement  fac¬ 
tors  are  all  less  than  unity  (have  negative  dB  values)  for  values  of  reso¬ 
lution  improvement,  I,  less  than  about  64.  This  means  that  the  "target 
peaks”  in  the  range  profile  auto-correlation  function  will  generally  be 
less  discernable  from  the  background  clutter  level  of  that  function  than 
the  frequency-averaged  target  plus  ciutier  return  amplitude  would  be  from 
the  mean  value  of  the  surrounding  clutter.  The  visibility  enhancement  fac¬ 
tors  do  increase  as  the  resolution  becomes  finer,  but  at  a  much  lower  rate 
than  in  the  case  of  coherent  processing.  Asymptotically,  the  contrast- 
enhancement  ratio  increases  only  in  proportion  to  /I  rather  than  linearly 
with  I  as  in  the  case  of  coherent  processing. 

Note  that  in  eqn  (3-19,  20)  the  resolution  improvement  factor,  I,  and 
the  number  of  resolved  scatterers,  M,  only  appear  together  in  the  term 
(I/M^).  Hence,  Fig.  3-7  for  the  case  M»2  can  be  used  to  obtain  the  target 
contrast  enhancement  ratio  for  any  value  of  M,  merely  by  defining  an 
equivalent  value  of  t  as  follows: 

I*  -  4(^-)  (3-21) 


Fig.  3-7  Range  profile  auto-correlation  function  Target  Contrast 
Enhancement  Ratio  as  a  function  of  Resolution  Improvement  Factor,  I,  for  a 
two  point-scat terer  target  in  clutter. 


82 


-  .  -w  *  .  -•  \ 


/ 

/ 


where  I  and  M  are  the  actual  resolution  improvement  ratio  and  number  of 
resolved  scatterers,  respectively,  and  I*  is  the  equivalent  value  of  I  to 
use  in  Fig.  3-7.  it  can  be  seen  that,  as  in  the  coherent  case,  the  contrast 
enhancement  ratio  decreases  as  the  number,  M,  of  resolved  scatterers 
increases. 

It  is  also  useful  to  plot  the  contrast  ratio,  P^/G,  of  the  target  peaks 
above  the  background  clutter  level  as  a  function  of  resolution  improvement. 
Such  a  plot  is  illustrated  in  Fig.  3-8  for  the  case  of  a  two-scatterer 
target.  Assume  that  the  target  peaks  have  to  be  at  least  3  dB  above  the 
average  background  level  to  be  considered  visible.  The  value  of 
single-pulse  T/C.  required  to  attain  this  level  of  contrast  at  any  given 
value  of  resolution  Improvement  can  then  be  obtained  from  Fig.  3-8.  For 
example,  in  order  for  the  mean  value  of  the  “target  peak”  of  a  two- 
scatterer  target,  for  which  T/C  **  0  dB,  to  attain  a  level  3  dB  above  the 
background,  a  resolution  Improvement  factor  of  at  least  32  is  required. 
Again,  the  target  peak/back'ground  contrast  ratio  for  values  of  M  other  than 
2  can  be  obtained  from  Fig.  3-8  by  using  the  transformation  defined  by  eqn 
(3-21). 

In  summary,  the  application  of  Incoherent  processing  to  obtain  a  high 
resolution  range-profile  auto-correlation  function  does  not  generally 
enhance  the  contrast  or  visibility  of  a  target  against  the  clutter 
background.  Specifically,  positive  enhancement  of  target-clutter  contrast 
is  only  realized  for  resolution  improvement  ratios  of  32  or  greater*  At 
lower  values  of  resolution  improvement,  target-clutter  contrast  is  actually 
degraded,  with  the  amount  of  degradation  Increasing  with  T/C.  The  reason  . 
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MEAN  TARGET  PEAK  ABOVE  MEAN  BACKGROUND  LEVEL 

(dB) 


RESOLUTION  IMPROVEMENT  FACTOR.  I 

Fig.  3-8  VI sibll icy  of  the  mean  targe t  peak  above  the  mean  background 
level  for  the  range  profile  auto-correlation  function  of  a  two  point- 
scatterer  target  in  clutter. 


for  this  phenomenon  is  that  the. clutter  background  level  of  the  auto¬ 
correlation  function,  against  which  the  target  peaks  are  contrasted,  con¬ 
tains  a  component  proportional  to  the  product  of  the  target  and  clutter 
reflectivities.  Aa  a  result,  this  background  level  also  Increases  with 
target  size.  Although  the  amount  of  contrast  degradation  increases  with 
T/C,  It  Is  generally  not  an  issue  for  strong  targets,  i.e.,  T/C  >  6  dB, 
that  are  strong  enuugh  to  be  visible  anyway,  despite  this  degradation. 

However,  for  weaker  targets,  for  which  T/C  is  between  0  and  6  dB,  the 
resolution  improvement  provided  by  the  pulse  train  bandwidth  should  be  at 
least  32,  in  order  to  keep  the  degradation  in  target-clutter  contrast  to 
less  than  0  dB  (for  a  two-scatterer  target). 

The  principal  advantage  of  incoherent  processing  i3  that  it  does  not 
require  motion  compensation  and  is  therefore  very  robust.  Uncompensated 
target  motion  does  not  affect  the  form  of  the  resulting  range  profile  auto¬ 
correlation  function,  nor  do  the  A/D  range-sampling  strobes  have  to  be 
synchronized  with  the  variable  frequency  source  to  insure  that  the  location 
of  the  target  peaks  in  the  auto-correlation  profile  is  unambiguous, 
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4.0  SIMULATION  RESULTS  FOR  RANGE-PROFILES  AND  RANGE-PROFILE 

AUTO-CORRELATION  FUNCTIONS  OBTAINED  BY  PROCESSING  RETURNS 

FROM  FREQUENCY-STEPPED  PULSE  TRAINS 

Chapter  2  demonstrated  how  high-resolution  range-profiles  and  range- 
profile  auto-correlation  functions  could  be  obtained  by  coherent  and  nonco¬ 
herent  processing,  respectively,  of  the  returns  from  a  frequency-stepped 
pulse  train.  It  was  shown  that  in  both  cases  the  resolution  of  the  resulting 
profile  was  equal  to  the  reciprocal  of  the  pulse  train  bandwidth,  NAf,  and 
that  in  order  to  avoid'  aliasing  or  ambiguities  in  the  resulting  profiles, 
the  Interpulse  frequency  spacing,  Af,  had  to  be  chosen  no  greater  than  one- 
half  the  radar's  single  pulse  bandwidth  i.e.  Af  £  Furthermore,  when 

the  pulse-train  frequency  was  linearly  stepped  (constant  Af),  the  required 
processing  in  both  cases  could  be  implemented  in  the  form  of  a  OFT.  The 
only  difference  between  coherent  and  noncoherent  processing  1.  that  in  the 
former,  the  OFT  of  the  complex  pulse  returns  (amplitude  and  phase)  is  taken 
while  in  the  latter,  the  OFT  of  the  squared-magnitudes  of  those  returns  is 
taken.  Finally,  it  was  shown  that  relative  target  radial  velocity  has  no 
effect  on  the  range  profile  auto-correlation  function  obtained  by  nonco¬ 
herent  processing,  but  that  ir.  the  coherent  case,  unless  motion  compen¬ 
sation  is  applied,  the  effect  of  such  velocity  will  be  to  shift  and  blur 
the  resulting  high  resolution  range  profile. 

In  chapter  3,  the  effect  of  such  processing  on  the  contrast  between 
the  target  and  the  clutter  background  in  the  corresponding  profile  was  exa¬ 
mined.  In  both  cases,  target  visibility  was  shown  to  Increase  with  the 
resolution  improvanwint,  I  »  NAf-r,  provided  by  the  bandwidth  of  the 
frequency-stepped  pulse  train.  Specifically,  In  the  coherent  case,  the 


visibility  of  the  target  in  the  range  profile  was  shown  to  increase 
directly  in  proportion  to  I/M,  where  M  is  the  number  of  resolved  scatterers 
appearing  in  the  profile.  As  a  result,  significant  increases  in  target 
visibility  can  be  realized  by  using  coherent  processing.  In  contrast,  it 
was  shown  that  noncoherent  processing  generally  results  in  a  degradation  of 
target  visibility  in  the  range  profile  auto-correlation  function,  unless 
the  resolution  improvement  factor,  I,  is  j^arge  enough  such  that  the  factor 
(I/M2)  is  greater  than  16. 

The  purpose  of  this  chapter  is  to  illustrate  some  of  the  results 
that  were  derived  in  the  preceding  chapters  in  terms  of  actual  high- 
resolution  range  profiles  and  range  profile  auto-correlation  functions. 

The  profiles  used  to  demonstrate  these  results  are  obtained  by  processing 
simulated  returns  from  frequency-stepped  pulse  trains. 

The  simulated  stepped-frequency  pulse  train  returns  to  be  processed 
are  generated  by  means  of  a  Monte-Carlo  simulation  of  the  radar  returns 
from  a  specified  target  against  a  clutter  background.  The  high  resolution 
range  profile  and  its  auto-correlation  function  are  then  obtained  by  pro¬ 
cessing  these  simulated  returns  using  the  methods  described  in  chapter  2. 

4.1  Description  of  Simulation 

'.'he  computer  simulation  for  generating  the  high-resolution  range  pro¬ 
file  and  its  auto-correlation,  using  the  radar  returns  from  a  frequency- 
stepped  pulse  train,  consists  of  the  following  three  parts. 

1.  Modelling  of  the  target-plus-clutter  radar  reflectivity  distribu¬ 
tion  over  the  extent  of  the  radar's  single-pulse  range-azimuth 


resolution  cell 


2.  Generation  of  the  complex  radar  pulse  returns  from  the 
target-plus-clutter  radar  reflectivity  distribution  defined  in  1. 
at  each  of  the  N  frequencies  used  in  the  pulse  train. 

3.  Generation  of  the  single-pulse-vidth  high  resolution  range  profile 
by  taking  the  FFT  of  the  resulting  N  stepped-frequency  complex 
pulse  returns,  and  the  corresponding  range  profile  auto¬ 
correlation  function  by  taking  the  FFT  of  the  squared  magnitudes 
of  those  pulre  returns. 

The  reflectivity  distribution  of  target  and  clutter  within  the  radar's 
single  pulse  range-azimuth  resolution  cell  is  modelled  as  follows.  The 
single  pulse  resolution  cell  is  subdivided  into  2048  subcells  by  sub¬ 
dividing  the  range  extent  of  the  cell  into  128  subcells  and  its  azimuth 
extent  into  16  subcells,  as  illustrated  in  Fig.  4-1.  The  total  range 
extent  of  128  subcells  actually  corresponds  to  two  6  dB  pulse  widths,  so  as 
to  include  any  returns  from  the  tails  of  the  pulse  shape.  Each  of  these 
2048  subcells  Is  assumed  to  contain  clutter  which  Is  modelled  as  an  inde¬ 
pendent  complex  Gaussian  process  for  each  cell.  The  variance  of  this 
clutter  process  is  uniform  for  all  subcells  and  Is  chosen  such  that  the 
total  clutter  power  received  from  all  2048  subcells  is  unity.  The  target 
is  modelled  by  a  collection  of  point  scatterers.  The  amplitudes  and  sub¬ 
cell  locations  of  these  target  scatterers  are  input  variables  to  the  simu¬ 
lation.  The  phases  of  these  target  point  scatterers  are  randomly  assigned. 
The  complex  reflectivity  of  each  subcell  is  then  weighted  by  the  pulse 
shape  and  beam  shape  profiles  as  a  function  of  the  subcell  location.  For. 
the  simulation  results  to  be  presented,  the  pulse  shape  profile  was  assumed 
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CROSSRANGE  (Azimuth) 


Gaussian  across  Che  range  extent  of  the  cell,  and  the  beam  shape  uniform 
over  the  azimuth  extent  of  the  cell. 

The  t&rget-to-clutter  ratio  is  defined  as  the  ratio  of  the  total 
received  target  power  (the  sums  of  the  squared  weighted  amplitudes  of  the 
target  scatterers)  to  the  weighted  total  clutter  power. 

The  complex  pulse  return  at  each  frequency  is  formed  by  coherently 
summing  the  weighted  complex  returns  from  each  of  the  2048  subcells.  This 
process  is  then  repeated  at  each  of  the  N  frequencies  used  in  the 
frequency-scepped  pulse  train.  The  resulting  set  of  N  frequency-diverse 
pulse  returns  is  then  Hamming-weighted  in  order  to  suppress  the  range  side- 
lobes  in  the  high  resolution  profile  to  be  generated. 

The  high  resolution  range  profile  of  the  cell  is  then  obtained  by  per¬ 
forming  an  N-point  FFT  on  the  N  Hamming-weighted,  frequency-diverse, 
complex  pulse  returns.  Since  the  output  of  the  FFT  is  generally  complex, 
.the  actual  range  picfile  is  obtained  as  the  magnitude  of  the  FFT  output. 

Similarly,  the  range  profile  auto-correlation  faction  is  obtained  by 
performing  an  N-point  FFT  on  the  Hamming  weighted  squared-magnitudes  of 
these  complex  pulse  returns.  Although  the  FFT  input  is  purely  real,  in 
this  case,  the  output  of  this  FFT  will  again  generally  be  complex.  Hence 
the  actual  range-profile  auto-correlation  function  is  again  obtained  as  the 
magnitude  of  the  FFT  output.  Note  that  the  only  difference  between 
generating  the  range  profile  and  its  auto-correlation  is  that  to  obtain  the 
auto-correlation  profile,  the  FFT  is  taken  of  the  squared  magnitudes  of  the 
complex  pulse  returns  rather  than  of  the  complex  returns  themselves. 
Strictly  speaking,  the  range-profile  and  Its  auto-correlation  ars  the 
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complex  outputs  of  the  corresponding  Fourier  transforms.  However,  we  are 
interested  only  in  the  magnitudes  of  those  profiles.  Hence  to  avoid  con¬ 
fusion  with  the  complex  profiles  in  the  sequel,  these  magnitude  Drofiles 
will  be  called  the  range  magnitude  profile  and  the  range  auto-correlation 
magnitude  profiles,  respectively.  It  can  be  shown  that  for  the  case  of  the 
postulated  clutter  random  process,  the  expected  value  of  the  range  auto¬ 
correlation  magnitude  profile  is  equal  to  the  auto-correlation  of  the 
expected  value  of  the  range  magnitude  profile. 

Note  that  because  the  clutter  model  used  in  the  simulation  is  a  random 
process,  the  magnitude  profiles  that  are  output  by  the  simulation  will  also 
be  random  processes.  An  estimate  of  the  mean  value  of  each  magnitude  pro¬ 
file  can  be  computed  as  the  sample  mean  value  of  that  profile  by  averaging 
the  sequence  of  magnitude  profiles  obtained  from  K  independent  trials  of 
the  simulation.  An  estimate  of  the  variance  of  each  profile  can  be  simi¬ 
larly  obtained  as  the  sample  variance  of  the  profile  computed  from  K  inde¬ 
pendent  trials. 

Therefore,  the  above  described  simulation  is  repeated  K  times  and  the 
K  resulting  profiles  are  averaged  to  obtain  the  sample  mean  and  sample 
standard  deviation  of  each  of  the  tifo  magnitude  profile  functions.  For  the 
results  presented, here,  K,  *  30. 

The  resulting  sample  mean  values  of  the  range  magnitude  profile  and  the 
range  auto-correlation  magnitude  profile  are  then  pldtted  as  a  function  of 
tange  and  range  separation  in  pulse  widths,  respectively.  In  addition  to 
the  sample  mean,  these  plots  also  show  the  value  of  the  sample  mean  plus' 
three  times  the.  sample  standard  deviation.  This  latter  quantity  is 
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displayed  in  order  to  portray  the  expected  variation  of  the  profiles  about 
their  mean  value. 

4.2  Simulation  Results 

Simulation  results  for  the  range  magnitude  and  range  auto-correlation 
magnitude  profiles,  computed  using  the  model  described  in  section  4.1,  will 
be  presented  for  the  following  cases. 

1.  Clutter  only. 

2.  A  stationary  two-point-scat terer  target  in  clutter  at  various  ' 
target-to-clutter  ratios. 

3.  A  stationary  three-point-scat terer  target  in  clutter. 

4.  A  stationary  two-polnt-scatterer  in  clutter,  observed  by  a  moving 
radar,  at  various  values  of  uncompensated  radial  velocity. 

5.  A  moving  two-polnt-scatterer  target  in  clutter  for  various  values 
of  uncompensated  radial  velocity. 

Each  of  the  mean-value  and  mean-value-plus-3-sigma  profiles 
illustrated  represents  the  sample  mean  of  30  independent  trials  of  the 
simulation.  The  profiles  are  plotted  in  terms  of  the  following  dimen¬ 
sionless  parameters:  profile  amplitude  in  dB  vs.  fine  range  or  range 
separation  in  pulse  widths. 

4.2.1  Clutter  Profiles 

The  . range  magnitude  and  range  auto-correlation  magnitude  profiles  for 
the  case  where  the  radar  resolution  cell  contains  only  uniform  clutter  are 
Illustrated  in  Figs.  4-2  and  4-3, for  the  case  of  a  64-frequency  pulse  train 
with  a  frequency  separation  between  pulses  of  1/2  the  single  pulse  bandwidth. 
This  frequency-separation  Insures  that  the  resulting  profiles  will  not  be 
corrupted  by  the  effects  of  aliasing.  Hence  the  resolution  Improvement 
factor  (i.e.  pulse  train  bandwidth  ratio)  for  this  case  is  1/2(64)  "  32. 
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64  POINT  FFT,  TW  PRODUCT  =  32.0 
GAUSSIAN  PULSE  SHAPE,  6  dB  WIDTH 
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RELATIVE  RANGE  (Pulse  Widths) 

Fig.  4-2  The  sample-mean  and  the  sampleTiean-plu8-3c  of  the  high  range 

IZ  ctn*  ir°fhile  °f  ClUtter  °nly  VUhta  *  SS. Tsll 

ra«L?!llv  1th  P  flle  shaPes  are  approximately  that  of  the  assumed 
uasian  pulse  shape.  Resolution  Improvement  Factor  1*32. 
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RANGE  SEPARATION  (Pulse  Widths) 

Pig*  4-3  The  sample-mean  and  sample^Bean-plus-3o  of  the  high  resolution 
range  auto-correlation  magnitude  profile  for  clutter  only  within  a  single 
pulse  resolution  cell.  Resolution  Improvement  Factor  £*32. 
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The  high  resolution  range  magnitude  profile  of  the  clutter  only  return 
illustrated  in  Fig.  4-2  reflects  the  Gausnian  pulse  shape  used  in  the  simu¬ 
lation.  It  can  be  seen  chat  the  6  dB  width  of  this  profile  is  approximately 
one  pulse  width  wide  as  expected.  Furthermore,  note  that  the  profile  of 
the  sample-mean-plus-three-times-the-sample-standard-deviation  is  approxi¬ 
mately  8  dB  above  that  of  the  sample  mean  profile.  This  8  dB  differential, 
which  is  uniform  over  the  entire  pulse  shape,  is  a  consequence  of  the  fact 
that  the  clutter  amplitude  is  Rayleigh  distributed.  Since  the  pulse  shape 
profile  is  due  to  the  return  from  clutter,  the  value  of  that  profile  at  any 
point  will  also  be  a  Rayleigh  distributed  random  variable.  It  can  be 
easily  shown  that  the  standard  deviation  and  mean  value  of  the  Rayleigh 
distribution  are  related  by: 


Mz  •  *  • 

Hence  the  ratio  of  the  mean-plus-3-times-the-standard-deviation  to  the  mean 

U:  i^.1+3/S.2.s7 

Uz  11 

Since  the  high  resolution  range  profile  is  an  amplitude,  rather  than  a 
power  profile,  the  value  of  this  ratio  in  dB  is  given  by: 

20  log(2.57)  -  8.19  dB 

which  is  in  good  agreement  with  the  observed  result  In  Fig.  4-2  that  the 
mean-plus-3-sigma  profile  is  uniformly  about  8  dB  above  the  mean  profile. 

.  This  phenomenon  will  also  be  evident  in  the  subsequent  profiles  of  targets 
in  clutter. 

The  corresponding  range  auto-correlation  magnitude  profile  is  shown 
in  Fig.  4-3.  It  consists  of  a  peak  at  zero  lag,  then  a  gradually 


decreasing  pedestal  out  to  a  lag  of  one  pulse  width.  The  height  of  the 
peak  above  the  pedestal  appears  to  be  about  7.0  dB.  This  is  confirmed  by 
the  theory  of  section  3.2  as  follows.  The  height  of  the  peak  at  zero  lag 
corresponds  to  the  total  power  in  the  cell,  which  in  this  case  is  just  C. 
The  height  of  the  clutter  background  pedestal,  G,  is  defined  by  eqn.  3-17, 
which  for  the  case  of  no  target  (ot  ■  0)  reduces  to 

G  -  /77A  ~  (4—1 ) 

Y I 

Hence  the  ratio  of  the  pedestal  height  in  Fig.  4-3  to  the  peak  value  is 


g  /77T 

R(0)  /f 


(4-2) 


Since  I  ■  32  for  the  profile,  of  Fig.  4-3,  the  pedestal  level  should  be 
1//32  of  the  peak,  or  7.5  dB  (10  log/32)  below  the  peak,  which  is  in  good 
agreement  with  the  observed  value  of  7.0  dB.  The  height  of  the  pedestal  of 
the  range  auto-correlation  magnitude  profile  below  the  peak  will  differ 
somewhat  from  the  theoretical  value  obtained  from  eqn.  (4-2)  because  eqn. 

I 

(4-2)  is  premised  on  the  assumption  of  a  uniform  pulse  ehape,  whereas  a. 
Gaussian  pulse  shape  was  used  in  the  simulntion,  and  also  because  of  the 
effect  of  the  Hamming-veighting  used  in  thu  simulation.  In  general 
however,  the  range  auto-correlation  magnitude  profiles  for  clutter  will 
have  the  shape  shown  in  Fig.  4-3,  with  the  relative  height  of  the  pedestal 
decreasing  in  proportion  to  the  square  roo :  of  the  pulse  train  resolution 
improvement  ratio,  as  indicated  by  eqn.  (4-2). 

As  in  the  case  of  the  range  profile,  :he  mean-plus-3-sigma  profile  for 
the  clutter  pedestal  is  seen  to  be  uniformly  about  4.2  dB  above  the  mean 
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profile  at  all  values  of  range  separation  outside  of  the  raainlobe.  This 
is  again  a  consequence  of  the  fact  that  the  probability  distribution  of  the 
magnitude  of  the  range  auto-correlation  function  pedestal  (the  region  out¬ 
side  of  the  mainlobe  and  excluding  any  target  peaks)  is  approximately 
Rayleigh,  as  previously  discussed  in  Section  3.2.  However,'  the  dB  figure 


for  this  offset  is  now  only  half  of  that  for  the  case  of  the  range  profile, 
since  the  auto-._or  relation  function  has  units  of  power  (amplitude-squared ) 
rather  than  amplitude.  Thus  the  fact  that  the  mean-pHs-3-sigma  profile  is 
uniformly  about  4.2  dB  above  the  mean  value  of  the  pedestal,  at  all  points 
of  the  pedestal,  confirms  the  analytical  result  obtained  la  Section  3.2 
that  the  pedestal  amplitude  is  approximately  Rayleigh  distributed. 

4.2.2  Simulation  Results  for  a  Stationary 

Two-Point  Scatterer  Target  in  Clutter 

Consider  a  target  consisting  of  two  equal  amplitude  point  scatterers 
that  are  separated  by  one  quarter  pulse-width  in  range.  Let  this  target  be 
symmetrically  located  about  the  center  of  the  radar's  single  pulse  range 
resolution  cell,  ’.et  the  amplitudes  of  these  two  scatterers  be  such  that 
the  average  unresolved  i:arget-to-clutter  ratio  of  the  combined  return 
(averaged  over  frequency)  from  both  scatterers  in  that  cell  is  about  3.7  dB. 
The  high  resolution  ran;;e  and  range  auto-correlation  magnitude  profiles  of 
this  target-clutter  combination  are  illustrated  in  Fig.  4-4  and  Fig.  4-5 
for  a  resolution  improvement  factor  of  16,  and  in  Fig.  4-6  and  4-7  for  a 
resolution  improvement  factor  of  64, 

The  high  resolution  range  magnitude  profiles  of  Figs.  4-4  and  Fig. 4-6 
clearly  resolve  these  two  target  scatterers  into  the  two  peaks  of  the  pro¬ 
file,  separated  by  one-quarter  pulse  width  and  symmetrically  located  about 


RANGE  SEPARATION  (Pulse  Widths) 

Fig.  4-5  The  sample-mean  and  sample-mean-plus**3o  of  the  high  resolution 
range  auto-correlation  magnitude  profile  for  a  two  point-scatterer  target 
in  clutter.  T/C  »3.97  dB,  Resolution  Improvement  Factor  1-16,  acatterer 
separation  is  0*25  pulse  widths. 
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RELATIVE  RANGE  (Pulse  Widths) 

Fig.  4-6  The  sample-mean  and  sample-mean-plus— 3cr  of  the  high  resolution 
range  magnitude  profile  for  a  two  polnt-scatterer  target  In  clutter. 

T/C  *3.72  dB,  Resolution  Improvement  Factor  1-64,  scatterer  separation  Is 
0.25  pulse  width. 
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Fig.  4-7  The  sample-mean  and  8ample-mean-plu8-3o  pf  the  high  resolution 
range  auto-correlation  magnitude  profile  for  a  two  poiut-scatterer  target 
in  clutter.  T/G  -3.72  dB,  Resolution  Improvement  Factor  1-64,  scatterer 
separation  is  0.23  pulse  width. 
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the  center  of  the  pulse  shape.  The  average  nelght  of  the  two  target  peaks 
above  the  mean  of  the  pulse  shaped  clutter  profile  should  be  equal  to  the 
resolved  target-plus-clutter  ratio:  T/C(Ar)  +  1  •  (I/M)  T/C(AR)  +  1  (see 
eqn.  3-12).  In  the  case  of  Fig.  .4-4  for  I  »  16,  the  peak  heights  of  the 
target  responses  are  seen  to  be  about  11.5  dB  above  the  pulse  shape.  The 
theoretical  value  of  the  resolved  target-plus-clutter  to  clutter  ratio  for 
this  case  is  (I/M)T/C(AR)  +  1  -  (16/2)(2.34)  +  1  -  19.8  -  13  dB.  This 
figure  has  to  be  corrected  for  the  processing  loss  of  1.34  dB  due  to  the 
Hamming-weighting  used  in  the  simulation.  When  this  correction  is  made,  . 
the  resulting  theoretical  value  of  the  relative  amplitude  of  the  profile 
peaks  above  the  clutter  is  found  to  be  11.6  dB,  in  excellent  agreement  with 
the  observed  value  of  11.5  dB. 

Similarly,  in  Fig.  4-6  the  two  target  peaks  are  seen  to  be  about  17.5  dB 
above  the  clutter  profile  for  1*64.  '  In  this  case  the  theoretical  resolved 
target-plus-clutter  to  clutter  ratio  is  18.8  dB.  Again,  subtracting  the 
processing  loss  for  Hamming-weighted  of  1.34  dB,  the  expected  height  of  the 
target  peaks  above  the  clutter  is  17.4  dB,  which  is  in  excellent  agreement 
with  the  observed  value  of  17.5  dB.  - 

In  the  corresponding  range  auto-correlation  magnitude  profiles  of 
Fig.  4-5  and  4-7,  the  pair  of  target-point  scatterers  is  represented  by  the 
single  peak  located  at  a  range  displacement  of  V4  pulse  width,  corresponding 
to  the  range  separation  of  the  two  scatterers.  The  observed  contrast 
enhancement  ratio  is  the  difference  between  the  average  level  of  this  peak 
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above  the  clutter  background  pedestal  and  the  single-pulse  target-plus- 
clutter  to  clutter  ratio.  In  Fig.  4-5,  for  I  *  16,  and  T/C  ■  3.7  dB,  the 
peak  height  above  the  pedestal  is  seen  to  be  about  3.0  dB. 

.'his  is  in  reasonable  agreement  with  the  theoretical  value  of  target 
peak  level  above  the  pedestal  of  about  3.4  dB,  for  1  »  16,  M  ■  2  and  T/C  ■ 

4  dB,  obtained  from  Fig.  3-8.  Similarly  in  Fig.  4-7,  for  I  »  64  and  T/C  ■ 

3.7  dB,  the  observed  value  of  the  average  target  peak  above  the  pedestal  is 

seen  to  be  about  6  dB,  whereas  the  theoretical  value  of  the  target  peak  for 
these  parameters  from  Fig.  3-8  is  about  6.2  dB.  Again  this  is  in  reaso¬ 
nable  agreement  with  the  observed  value  of  6  dB.  The  small  differences  in 

the  target  peak  contrast  between  the  theoretical  values  obtained  from  Fig. 

3-8  and  the  simulation  results  of  Figs.  4-5,  4-7  are  due  to  the  following  2 
effects.  In  contrast  to  the  theoretical  annlysis,  the  simulation  uses 
Hamming-weighting  and  assumes  a  Gaussian  pulse  shape.  The  effect  of 
Hamming-weighting  is  to  lower  the  target  peaks,  while  the  use  of  a 
Gaussian,  rather  than  a  uniform,  pulse  shape  slightly  lowers  the  value  of 
the  pedestal. 

The  final  example  is  Intended  to  illustrate  the  difference  between 
coherent  and  incoherent  processing  in  detecting  targets  at  very  low  target 
to  clutter  ratios.  The  range  magnitude  and  range  auto-correlaton  mangltude 
profiles  are  again  computed  for  the  same  two-scatterer  target  model  as 
before,  but  at  a  relatively  low  target-to-clutter  ratio  of  about  -2  dB,  and 
using  a  64-frequency  pulse  train  having  a  resolution  Improvement  factor  of  32 
.  The  range  magnitude  profile  for  this  case  is  Illustrated  in  Fig.  4-8. 
Although  the  unresolved  target-to-clutter  ratio  is  only  -2  dB,  coherent 


S/C  =  -2  07  dS  BETA  =  1.0  P  =  0.0  30  TRIALS 

64  POINT  FFT.  TW  PRODUCT  =  32.0 
GAUSSIAN  P'JLSE  SHAPE.  6  dB  WIDTH 


RELATIVE  RANGE  (Pulse  Widths) 

Pig.  4-8  The  sample-mean  and  sample-mean-plus-So  of  the  high  resolution 
range  magnitude  profile  for  a  two  point-scatterer  target  in  duttfer. 

T/C  —2.07  dB,  Resolution  Improvement  Factor  1-32,  scetterer  separation  Is 
0.25  pulse  width. 
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RANGE  SEPARATION  (Pulse  Widths) 

Fig,  4-9  The  sample-mean  and  sample-mean-plus-3o  of  the  high  resolution 
range  auto-correlation  magnitude  profile  for  a  two  point-scatterer  target 
in  clutter.  T/C  "-2. 07  dB,  Resolution  Improvement  Factor  1*32,  scatterer 
separation  is  0.25  pulse  width.  Note  that  the  target  peak  at  Ar  -0.25r  is 
barely  discernible  from  the  clutter  background. 
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processing  raises  the  average  value  of  the  target  peaks  to  a  level  of  about 
7.5  dB  above  the  mean  clutter  level.  The  theoretical  value  for  the 
contrast  between  the  resolved  target  peaks  and  the  clutter  background  is 
again  given  by:  (I/M)T/C(AR)  +  1  *  32/2(0.63)  +  1  •  11.08  or  10.4  dB,  less  the 
processing  loss  due  to  Hamming-veighting  of  1.34  dB,  or  9.0  dB,  slightly 
more  than  the  value  obtained  from  the  simulation.  Note  that  although  the 
presence  of  this  small  target  (T/C—2  dB)  would  be  scarcely  visible  against 
the  clutter  in  the  unresolved  case,  the  resolution  Improvement  provided  by 
coherent  processing  really  makes  the  presence  of  the  two  scatterers  stand 
out  above  the  clutter  in  the  resulting  high  resolution  profile. 

The  corresponding  range  1  auto-correlation  profile  for  this  low 
target-clutter  case  is  illustrated  in  Fig.  4-9.  In  contrast  to  the  high 
resolution  range  profile  of  Fig.  4-7,  the  target  peak  located  at  the 
quarter  pulse-width  point  in  the  range  magnitude  auto-correlation  is  barely 
discernable  from  the  clutter  background  pedestal.  Its  level  appears  to  be 
about  1  dB  above  the  clutter  background  pedestal,  which  la  on  the  same  order 
as  the  residual  noise  in  the  sample-mean  amplitude  of  the  pedestal.  This 
la  in  reasonable  agreement  with  the  theoretical  value  of  the  target  peak  ' 
level  of  about  1.7  dB  obtained  from  Fig.  3.8  for  the  indicated  values  of 
I,  M,  and  T/C.  This  example  clearly  shows  the  increase  in  target  contrast 
at  low  value*  of  T/C  provided  by  coherent  processing  relative  to  that  pro¬ 
vided  by  incoherent  processing. 

4.2.3  Simulation  Results  for  a  Stationary  Three-Polnt-Scatterer 
Target  in  Clutter 

Consider  a  target  consisting  of  three  equal  amplitude  point- 
scatterers,  positioned  in  range  such  that  their  range  separation*  are  ' 
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non-equal.  Specifically,  let  the  first  scatterer  be  located  1/8  of  a  pulse 
width  forward  of  the  center  of  the  pulse  return,  the  second  1/16  of  a  pulse 
width  aft  of  the  center,  and  the  third  1/8  of  a  pulse  width  aft  of  the 
center.  Consequently,  the  range  separations  for  the  different  scatterer- 
pair  combinations  will  be  as  follows:  3/16  pulse  width  between  the  first 
and  second  scatterers,  1/4  pulse  width  between  the  first  and  third  scat- 
terers  and  1/16  pulse  width  between  the  second  and  third  scatterers. 

The  corresponding  high  resolution  range  magnitude  profile  of  this 
three-point-scatterer  target  in  clutter  is  illustrated  in  Fig.  4-10  for 
T/C  *  5.84  dB  and  a  resolution  improvement  of  1  ■  64.  The  target  peaks  are 
seen  to  all  be  about  17.5  dB  above  the  pulse-shaped  clutter  background  pro¬ 
file.  The  theoretical  value  for  the  contrast  between  the  resolved  target 
peaks  and  the  clutter  back  in  this  case  is  (64/3)(3.83)+l  •  82.85  or  19.18 
dB.  When  this  is  corrected  for  the  1.34  dB  Hamaing-welghting  loss,  this 
theoretical  result  is  again  in  excellent  agreement  with  the  observed  value 
of  17.5  dB. 

The  corresponding  range  auto-correlation  magnitude  profile  is 
illustrated  in  Fig.  4-11.  There  are  three  peaks  corresponding  to  the  three 
unique  values  of  range  separation  between  any  two  of  the  three-point-scat- 
terers  of  the  target.  Specifically,  since  the  target  scatterers  are 
located  at  ranges  of  -1/8  pulse  width,  +  1/16  pulse  width  and  +1/8  pulse 
width,  relative  to  the  center  of  the  pulse  return,  the  corresponding  range 
separations  between  them  are  3/16,  1/4  and  1/16  of  a  pulse  width.  It  is  at 
these  values  of  range  separation  that  the  target  peaks  of  the  range  .auto¬ 
correlation  profile  shown  in  Fig.  4—11  are  located.  The  height  of  these 
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range  magnitude  profile  for  a  three  polnt-scatterer  target  In  clutter. 
T/C  *5.84  dB,  Resolution  Improvement  Factor  1*64,  scatterer  separations 
are;  .078  pw,  .172  pw,  and  0.25  pw. 
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Pig.  4-11  The  saaple*naeah  and  sampleisean-plus-So  of  the'  high  resolution 
range  auto-correlatlon  magnitude  profile  for  a  three  polnt-scatterer  target 
In  clutter  T/C  -5.84  dB,  Resolution  Improvement  Factor  I*64,  scatterer 
separations  are  0.078  pv,  .172  pv,  and  0.25  pw. 


target  peaks  above  the  clutter  background  pedestal  is  seen  to  be  about  5.2 
dB.  This  is  in  good  agreement  with  the  theoretical  value  of  5.5  dB  obtained 
from  eqn.  (3-23)  for  the  current  parameter  values  of  I  ■  64,  M .»  3  and  T/C  ■ 
5.84  dB.  This  predicted  value  of  the  target  peak  above  the  pedestal  can 
also  be  obtained  from  Fig*  3-7  for  the  two-point  scatterer  case  by  using 
the  equivalent  value  of  I  for  two  scatterers,  eqn.  (3-21),  computed  as  follows: 

**  ■  if?  «>2  ■  7S7T  «>2 ' 28-5 

In  general,  the  number,  Np,  of  discrete  target  peaks  appearing  in  the 
range  auto-correlation  magnitude  profile  will  be  less  than  or  equal 
M(M-l)/2,  where  M  is  the  number  of  discrete  point  scatterers  comprising  the 
target  and  M(M-i)/2  is  the  number  of  different  pair  combinations  that  can 
be  obtained  from  M  different  objects.  The  actual  number  of  target  peaks 
obtained  will  only  be  equal  to  M(M-l)/2  if  the  scatterer  range  separation 
is  unique  for  each  of  the  M(M-l)/2  pairs  of  M  point  scatterers  as  in  the 
preceding  example. 

In  contrast,  consider  the  case  where  the  three-polnt-scatterers  of  the 
target  are  uniformly  spaced  in  range  i/8  of  a  pulse  width  apart  and  sym¬ 
metrically  located  about  the  center  of  the  pulse.  The  high  resolution 
range  amplitude  profile  for  this  case  is  illustrated  in  Fig.  4-12.  As 
before,  the  three  target  peaks  corresponding  to  these  scatterers  are 
clearly  visible  in  this  high  resolution  range  profile. 

The  corresponding  range  auto-correiatlon  amplitude  profile  is  shown  in 
Fig.  4-13.  However,  now  there  are  only  two  target  peaks  rather  than  three, 
one  at  1/8  pulse  width  and  the  other  at  1/4  pulse  width  range  separation. 
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.Fig.  4-12  The  saople-iaean  and  8ample-oean-plaa-3o  of  the  high  resolution 
range  magnitude  profile  for  a  three  point-scatterer  target  in  clutter 
T/C  *5.84  dB,  Resolution  Improvement  Factor  I«64,  equally  spaced  seat— 
terers,  scatterer  separation  0. 125  pw. 
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Also,  the  peak  at  1/8  pulse  width  separation  is  1.5  dB  higher  than  that  at 
1/4  pulse  width.  This  happens  because  the  three  scatterers  are  equally 
spaced  1/8  pulse  width  apart.  Consequently,  two  of  the  possible  three 
scatterer  pairs  both  have  the  sane  range  separation  of  1/8  pulse  width  and 
as  a  result,  the  amplitude  of  the  target  peak  in  the  auto-correlation  pro¬ 
file  located  at  1/8  pulse  width  is  about  1.5  dB  higher  than  that  at  1/4 
pulse  width  range  separation.  The  "double**  peak  at  1/8  pulse  width  is  only 
1.5  dB  higher  than  that  at  1/4  pulse  width,  rather  than  3  dB  higher,  as 
might  have  been  expected.  This  happens  because  the  contribution  due  to 
multiple  scatterer  pairs  having  the  same  range  separation  adds  approxima¬ 
tely  as  the  square-root  of  the  sums  of  the  squares  of  the  cross-sections 
(<jt/M)  of  the  individual  point  scatterers,  in' accordance  with  eqn  (3-20) 
for  the  target  peak  amplitude  Pt.  Hence  the  amplitude  of  a  "double"  peak 
is  approximately  proportional  to: 


'  As  previously  described  in  section  2.2.3,  the  two  effects  of  uncompenr> 
sated  relative  radial  motion  between  the  target  and  the  radar  on  the  high- 
resolution  target  range  profile  are  a  range  shift  proportional  to  the 
radial  velocity,  and  a  dispersion  of  the  profile  because  of  filter  mismatch 
effects  when  the  uncompensated  velocity  becomes  sufficiently  large.  These 
two  effects  can  be  quantified  in  terms  of  the  dimensionless  velocity  para¬ 
meter,  P,  defined  by  eqn.  (2-36). 
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The  range  shift,  Sr,  due  to  uncompensated  radial  velocity,  of  the  high 
resolution  profile  was  previously  given  in  terms  of  the  number  of  high 
resolution  bins  by  eqn.  (2-39): 


LTf 
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(2-31) 


where  f c  is  the  center  frequency  of  and  B  is  the  frequency  spread  across 
the  pulse  train.  Since  the  high  resolution  range  profiles  illustrated  in 
this  chapter  are  presented  in  units  of  pulse  widths,  it  is  convenient  to 
also  express  the  velocity  induced  range  shift  defined  by  eqn.  (2-31)  In 
those  units.  This  can  be  done  by  noting  that  the  relation  between  the 
single  pulse  range  resolution  (i.e.  the  pulse  width)  and  the  high  resolu¬ 
tion  range  subcell  width,  Ar,  is  given  by  eqn.  (3-1).  Therefore,  combining 
eqn.  (2-31)  with  (3-1)  gives  the  following  result  for  the  velocity  induced 
range  shift  in  units  of  pulse  widths. 
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If  one  defines  the  dimensionless  frequency  ratio:  f2  =  fc/B,  then  Srp  can  be' 
expressed  in  terms  of  the  three  dimensionless  parameters t  ft,  P,  and  I.  as: 


srp  *  fzx‘ 


(4-3a) 


Furthermore,  by  using  the  definition  of  P  given  in  eqn.  (2-36)  and  the 
relation  between  the  target  doppler  frequency  and  target  radial  velocity: 
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Srp  can  also  be  expressed  in  terms  of  the  target-doppler  frequency  and  pulse 
repetition  Interval,  T,  as: 

Srp  -  -f  (4-4) 

In  eqn.  (4-4),  the  product,  NT,  Is  the  duration  of  the  pulse  train  and  the 
ratio  (N/I)  is  the  pulse  train  frequency  over-sampling  ratio. 

The  dispersion  of  the  profile  is  strictly  a  function  of  the  parameter 
P,  and  is  generally  negligible  for  values  of  P<3.  These  effects,  of  course, 
disappear  if  the  motion  compensation  phaae  corrections  are  applied  to  the 
pulse  returns  prior  to  computation  of  the  high  resolution  range  profile. 

As  previously  discussed  in  section  2.3,  uncompensated  radial  velocity 
generally  has  no  effect  on  the  range  auto-correlation  profile,  with  one 
exception,  to  be  described  in  the  sequel. 

These  effects  are  illustrated  by  simulation  results  for  a  moving  two- 
point  scatterer  target  in  clutter.  The  model  of  the  two-point  scatterer 
target,  in  these  simulations  is  the  same  as  that  used  previously  in  the  sta¬ 
tionary  target  simulation  of  section  4.2.2. 

Two  distinct  cases  are  described.  In  the  first,,  the  target  is  sta¬ 
tionary  and  the  radar  is  moving.  In  the  second,  the  radar  Is  stationary 
and  the  target  is  moving  radially  through  the  clutter.  The  essential  dif¬ 
ference  between  these  two  cases  is  that  in  the  first,  both  the  target  and 
the  clutter  have  the  same  radial  velocity  relative  to  the  radar,  whereas  in 
the  second  only  the  target  moves  relative  to  the  radar  while  the  clutter  is 
stationary. 


The  high  resolution  range  and  range  auto-correlation  magnitude  profiles 
for  the  moving  radar  case  are  illustrated  in  Fig.  4-14  and  4-15,  respec¬ 
tively  for  P  •  0.5  and  iz  * 67 .  These  profiles  should  be  compared  with 
those  for  the  corresponding  stationary  target  case  illustrated  in  Figs.  4-6 
and  4-7.  It  Is  seen  that  the  high  resolution  range  profile  of  Fig.  4-14  is 
merely  a  circularly  shifted  replica  of  that  for  the  stationary  target  case 
shown  in  Fig.  4-6,  whereas  the  range  auto-correlation  profiles  are  sta¬ 
tistically  the  same  in  both  the  fixed  (Fig.  4-7)  and  the  moving  target 
(Fig.  4-15)  cases,  as  expected.  The  shift  of  the  high  resolution  range 
profile  (Fig.  4-14  relative  to  4-6)  is  0.55  pulse  widths  as  predicted  from 
eqn.  (4-3a)  for  the  given  parameters  P  “0.5,  ft  “  67.  and  I  *  64. 

The  simulation  is  now  repeated  for  the  case  of  a  fixed  radar  with  the 
target  moving  through  the  clutter,  for  the  same  target-clutter  model.  The 
resulting  high-resolution  range  and  range-auto-correlation  profiles  are 
'  illustrated  in  Figs.  4-16  and  4-17.  The  difference  between  the  high 
resolution  range  profile  of  Fig.  4-16  and  that  of  Fig.  4-14  for  the  moving- 
radar  case  is  that  in  Fig.  4-16  only  the  target  response  has  been  shifted, 
whereas  the  high-resolution  pulse-shaped  clutter  profile  remains  fixed. 

This  reflects  the  fact  that  the  target  is  moving  relative  to  the  clutter.  . 
Note  that  the  amount  of  the  target  response  shift  is  the  same  as  that  in 
Fig.  4-14.  Furthermore,  because  of  this  relative  motion  between  the  target 
and  the  clutter,  the  resulting  range  auto-correlation  profile  of,  Fig.  4-17 
1*  somewhat  different  from  those  of  Figs.  4-7  and  4-15.  However,  the  loca¬ 
tion  of  the  target  peak  in  these  auto-correlation  profiles  is  the  same  in 


all  cases 


TARGET  PEAKS 


RELATIVE  RANGE  (Pulse  Widths) 

Pig,  4-14  The  sample-mean  and  saaple-mean-plus— 3a  of  the  high  resolution 
range  magnitude  profile  of  a  two  point-scatterer  target  in  clutter  observed 
by  a  moving  radar.  No  motion  compensation  applied.  As  a  result,  both 
target  and  clutter  profiles  appear  shifted  in  range  by  an  amount  propor¬ 
tional  to  the  radar's  radial  velocity.  Conditions  Identical  to  those  of 
Pig,  4-6. 
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RELATIVE  RANGE  (Piilso  Widths) 

Fig.  4rl6  The,  sample-mean  and  saaple-mean-plus-3o  of  the  high  resolution 
range  aagnitude  profile  for  a  two  polnt-scatterer  target  aovlng  through 
clutter,  observed  by  a  stationary  radar.  Only  the  profile  of  the  two 
target  scatterers  is  shifted  by  an  aoount  proportional  to  the  target  velo¬ 
city.  Conditions  identical  to  those  of  Fig.  4-6. 
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Fig.  4-17  Th«  sample-mean  and  sample— mean— plua-3c  of  the  high  resolution 
range  auto-cor relation  magnitude  profile  for  a  two  polnt-scatterer  target 
moving  through  clutter,  observed  by  a  stationary  radar.  Conditions  Iden¬ 
tical  to  those  of  Fig.  4-6. 


In  summary,  an  uncompensated  radial  velocity  will  first  result  in  a 
shift  of  the  high  resolution  range  profile,  and  if  large  enough,  will  even¬ 
tually  cause  some  dispersion  of  that  profile.  Furthermore,  If  there  is  any 
radial  motion  of  the  target  relative  to  the  ground  clutter,  the  target 
response  will  appear  shifted  with  respect  to  that  of  the  clutter  in  the 
resulting  high  resolution  profile.  The  clutter  response  in  that  profile 
can  usually  be  identified  by  the  fact  that  it  has  the  approximate  shape  of 
the  transmitted  pulse.  On  the  other  hand,  uncompensated  radial  velocity 
appears  to  have  no  effect  on  the  resulting  high-resolution  range  auto¬ 
correlation  profile,  unless  there  is  significant  radial  motion  of  the 
target  relative  to  the  surrounding  clutter.  However,  even  in  that  case,  the 
location  of  the  target  peaks  in  the  profile  appears  unchanged,  and  the  prin¬ 
cipal  effect  is  to  modify  the  shape  of  the  clutter  background  pedestal  of 
the  auto-correlation  profile. 
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Statistics  of  the  Discrete  Range  Profile 
Auto-correlation  Function 


Consider  a  discrete  high-resolution  range  profile ,  made  up  of  I  con¬ 
tiguous  range  subcells.  Each  of  the  I  subcells  contains  either  clutter 
only,  or  a  discrete  point  target  plus  clutter.  There  are  a  total  of  M 
discrete  point  targets  arbitrarily  distributed  over  the  I  subcells;  where  M 
is  less  than  I.  The  complex  return  from  the  clutter  in  the  ith  subcell  is 
given  by: 

qi  “  aA  e3$i  i-1,...  I  (A-l) 

where  is  a  Rayleigh  distributed  random  variable  and  ^  is  a  uniformly 
distributed  phase  angle  over  -»,+*.  Similarly  the  complex  return  from  the 
point  scatterer  in  the  JQth  subcell  is: 

qjB  "  V  e-J®n  ,  m-1,...  M  (A-2) 

where  V  is  the  magnitude  of  the  voltage  return,  from  the  mth  scatterer, 
which  is  assumed  to  be  Constant  aud  equal  for  all  m  sratterers,  6n  is  the 
phase  of  the  return  from  the  mth  scatterer  and  ja  is  the  subcell-location 
1 <1m<I  of  the  mth  scatterer.  The  target  phases,  0a,  are  assumed  to  be 
arbitrary. 

The  discrete  complex  auto-correlation  function  of  the  target-plus- 

„ 

clutter  range  profile  described  by  (A-l)  and  (A-2)  Is  therefore: 


RqOO 
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\  qi+k 


(A-3) 


where  qj  Is  given  either  by  (A-l)  for  the  subcells  containing  clutter  only, 
or  by  the  sum  of  (A-l)  and  (A-2)  for  the  M  subcells  that  contain  a 
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combination  of  a  discrete  point  scatterer  plus  clutter.  When  (A-l)  and 
(A-2)  are  substituted  into  (A-3)  the  following  result  is  obtained: 


M 


M 


Ra(k)  *  I  {<ai  +  v  1  e“Jem  d(i-jn))(ai+k  eJ^i+k  +  V  V  e39n  6(i+k-jn>)} 

H  I— l  m-1  n-  1 

(A- 3a) 


1  i-J 

where:  5(i-j)  =  { 

0  i*j 

Expanding  (A- 3a)  and  making  use  of  the  properties  of  the  5  function  yields: 


I  M 

RqO°  "  iI1  al  ai+k  eJ^i+k-^i)  +  Vj^  ajmHt  ejOa-^-k) 

M  M  M 

+  V  \  aj  eJ  jm+k"flm)  +  V2  \  \  eJ^®n“®a)  6(ja-jn+k)  (A- 4) 

a-l  a-l  n-1 

It  Is  seen  that,  in  general,  Rq(k)  Is  complex.  However,  as  will  be 

seen  shortly,  Rq(0)  is  real  and  positive,  as  it  should  be,  since  the  value 

of  the  auto-correlation  of  a  function  at  zero  lag  is  sinply  the  mean- 

squared  value  of  that  function. 

The  expression  fbr  Rq(k)  given  by  (A-4)  can  be  decomposed  into  three 
distinct  components: 

1.  The  auto-correlation  between  the  returns  due  to  clutter  only. 

This  is  given  by  the  term  £  *i  *i-Ht  e-3  ( ♦i+k“4i) . 

1 

2.  The  cross-correlation  between  the  returns  froa  the  M  target  point 
scatterers  and  clutter  returns.  This  ip  given  by  the  second  and 
third  terms  of  (A-4)  which  may  be  lumped  into  a  single  sum  over 
m— 1 ,  ...  M. 

3.  The  auto-correlation  between  the  target  returns  only.  This  is 
given  by  the  last  term  of  (A-O,  (he  double  sum  over  a  and  n. 


Since  the  clutter  amplitudes  a^  and  the  clutter  and  target  phase 
angles  and  9jm  are  all  random  variables,  Rq(k)  will  be  a  random  function 
of  the  parameter  k.  We  desire  to  determine  the  statistics  of  the  magni¬ 
tude  of  this  auto-correlation  function,  |Rq(k)|.  As  will  be  seen,  the 
distribution  of  |Rq(k)|  and  its  statistics  are  functions  of  the  corre¬ 
lation  lag  k. 

Specifically,  the  distribution  of  Rq(k)  varies  in  accordance  with 
three  different  regions  as  follows: 

1.  k-0:  the  distribution  of  Rq(k)  in  this  region  is  approximately 
non-zero  mean  normal  for  large  I. 

2.  k-jn-jm*0:  the  distribution  of  |Rq(k)|  in  this  region  is  approxi¬ 
mately  Ricean. 

3.  k*0  and  k*|jm-jn|:  the  distribution  of  Rq(k)  in  this  region  is 
approximately  complex-normal  with  zero  mean,  for  large  I. 

Therefore  the  distribution  of  |Rq(k)|  is  approximately  Rayleigh. 

Consider  first  the  case  k*0.  In  this  special  case,  the  expression  for 

Rq (k)  reduces  to: 

*  ▼  M 

R  (0)  -  I  V*  2V  l  a-j  cos<em-t4  )  +  MV2  (A-5) 

■  i-1  m-1  Ja 

Since  the  a*  are  Rayleigh  distributed,  the  a^2  are  exponentially  distri¬ 
buted.  Therefore,  because  of  the  assumed  independence  of  the  clutter 
returns  from  different  subcells,  the  distribution  of  the  first  term  of  ' 
(A-S)  is  chi-squared  with  21  degrees  of  freedom,  which  is  approximately 
normal  for  large  I.  The  mean  value  of  the  a^2  is  merely  the  clutter  power 
returned  from  the  1th  subcell,  which  is  the  total  clutter  power  divided  by 
the  number  of  range  subcells.  Thus, 

E[ai2]  -  C/I.  (A-6a) 
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The  second  term  has  a  complicated  distribution  but  is  zero  mean. 
However,  when  it  is  added  to  the  first  term,  the  resulting  sum  will  still 
be  approximately  normal.  The  third  term  is  merely  a  constant.  Furthermore 
note  that  V2  is  the  power  returned  from  a  single  point  scatterer.  Since 
each  of  the  target  point  scatterers  was  assumed  to  have  equal  magnitude, 
the  power  from  each  scatterer  is  simply  the  total  target  cross-section' 
divided  by  the  number  of  discrete  point  scatterers: 

V2  -  at/M  (A- 6b) 

Hence  the  expected  value  of  Rq(0)  is: 

E[Rq(0)  I  «  I  <-)  +  0  +  M  <^) 

E[Rq(0)]  -  C  +  <rt  (A-7) 

or  as  expected,  E[Rq(0)]  is  simply  the  sum  of  the  mean-square  clutter  and 
target  power  returned  from  the  entire  resolution  cell. 

Tha  variance  of  Rq(0)  can  be  easily  obtained  as  the  sum  of  the  varian¬ 
ces  of  the  individual  terms  in  (A-5)  since  all  the  terms  are  uncorrelated. 

2  1  M 

Thus  on  (o)  *  I  VarU*2)  +  4V2  £  E(aj  2)  E(cos2(8n-<fr.i))  (A-8) 

^  i-1  m-l  m  m 

Since  the  a^2  are  exponentially  distributed,  [AlJ: 

VarCaj2)  *  E2(ai2)  ■  (C/I)2  . . 

Furthermore  for  0m  and  uniformly  distributed  on  it  follows  that: 

E[cos2(9m-<j»jffl)]  -  1/2 

Thus  «Rq2(0)-  I  (f)2  +  2(|%M  (|) 

or  0.  (0)  -  c/(l+2ot/C)/I  (A-9) 
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Hence,  when  k*0,  Rq(0)  is  approximately  normal  with  a  mean  value  given  by 

(A-7)  and  a  standard  deviation  given  by  (A-9). 

Next  consider  the  statistics  of  |Rq(k)|  for  the  case  where:  k*0  and 

k^  j j m- j n  | ,  where  m  and  n  are  the  inc|ices  of  two  distinct  point-scatterers 

of  the  target.  In  this  case,  eqn.  (A-4)  reduces  to: 

I 

Rq(k)  -  a±  aj+fc  e3($i+k"$i) 

+  V  l  (ajm-k  eJOa  “  <frjm-k)  +  »jm+k  eJ^jm+k  "  0m>  (A- 10) 

k*0,  k*  |jn-jQ } 

the  last  term  in  (A-4)  vanishes  because  6(Jm-jQ+k)  -  0  when  k*  |ja-jn| . 

The  first  term  in  (A- 10)  is  the  sum  of  products  of  independent,  zero 
mean  complex  normal  random  variables.  For  large  I,  the  distribution  of 
this  sum  term  will  be  a  zero  mean,  complex,  approximately  normal  random 

variable.  The  second  term  ie  (A-10)  is  the  sum  of  zero  mean,  complex  nor- 

»  •  ' 

mal  random  variables.  This  follows  because  the  a^  are  Rayleigh  distributed 

and  the  arguments  of  the  complex  exponential  are  independent  of  the  a^  and 
uniformly  distributed.  Hence  the  second  term  is  itself  a  zero  mean, 
complex  normal  rv  and  so  the  sum  of  the  two  terms  comprising  Rq(k)  has  a 
distribution  that  is  zero  mean,  complex,  and  approximately  normal  for  large 
I.  A  zero  mean  complex  normal  distribution  is  characterized  by  only  a 
single  parameter  -  its  variance.  This  variance  may  be  easily  found  by 
decomposing  Rq(k)  into  its  real  and  imaginary  parts  and  then  computing  the 
variance  of  either  part.  Denote  this  variance  «s  cx*.  Taking  the  real 
part  of  Rq(k): 

Re[Rq(k)J  •  l  at  an*  cos  (4i+k"ti) 
i-1 

M 

♦  v  I  cos  <®m“*Ja-k>  t  aJa+k  cos  <HB+k-8m>)  U-ll) 

m*l 
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Because  Che  a^  and  the  differential  phase  angles  in  the  above  are  all 
assumed  to  be  Independent  r.v.s,  all  of  the  terms  in  the  above  summation  are 
mutually  independent,  and  therefore  uncorrelated.  In  fact  even  if  the 
were  not  independent,  the  fact  that  the  differential  phase  angles  are  and 
that  the  expected  values  of  their  cosines  are  all  zero  would  be  enough  to 
guarantee  that,  the  individual  terms  in  the  above  summation  are  all  unco r re¬ 
lated.  As  a  result,  the  variance  of  (A- 11)  is  the  sum  of  the  variances  of 
the  individual  terms.  Furthermore,  the  individual  terms  are  all  zero  mean 
and  the  differential  phase  angles  are  independent  of  the  a^,  hence: 

I  « 

°x2  *  I  EIai2  ai+k2)  (1/2)  +  V2  l  1/2  (E[aj  _fc2]  +  E[aj _+k2]) 
i»l  m-1 

-|l  E2[ai2]  +  V2  M  EUi2] 

Finally  since  Ela*2]  -  C/I  and  V2  -  ot/M  the  above  can  be  written  as: 

°rt2  "yC2  (l+2ot/C)/I  .  (A- 12) 


which  is  equal  to  exactly  <t2r<j(0)/2  (see  eqn.  A-9). 

Hence  we  now  know  that  when  k*0  and  k#  Jjm'Jn! »  Rq(k)  is  a  zero  mean, 
complex  normal  random  variable  whose  variance  is  given  by  (A- 12) . 
Consequently  ,lt  follows  that  the  magnitude  |Rq(k) j  is  Rayleigh  distributed 
with  mean  value  and  standard  deviation  ox •  Thus: 


B£  jRqCk)  |  J  -  Q 


(A-13a) 


®  |Rq(k)  j  -y^C  /  (l+2ot/C)/I 


(A- 13b) 
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Finally,  where  k*0  and  k»|jm-jn|,  the  expression  for  Rq(k)  in  (A-4)  reduces 
to  that  of  (A- 10)  plus  the  constant  term  V2  e^®m~®n^.  Thus  if  the 
expression  for  Rq(k)  given  in  (A- 10)  is  designated  as  Q(k),  the  expression 
for  Rq(k)  when  k*0  and  k“|jm-jn|  can  be  written: 

Rq(k)  -  Q(k)  +  V2  e^nfV  (A-14) 

where  Q(k)  is  giyen  by  (A- 10). 

UoweVer,  we  have  previously  established  that  Q(k)  has  a  zero  mean,  complex 

normal  distribution  with  variance  given  by  (A-12).  Thus  the  distribution 

of  Rq(k)  in  this  case  is  complex  normal  with  a  non-zero  mean  value  given  by 

v2  ej<0B-V.  Consequently,  the  distribution  of  |Rq(k)|  in  this  case  is 

Ricean  [A2 ] ,  with  mean  given  by: 

E[|Rq(k)|]  “  /f-oxe-vA/4°x2  {(1  +~2*>  *0  <7^> 

2arx  4ox 


(A- 15) 


where  I{)(x),  IjCx)  are  the  modified  Bessel  functions  of  order  0  and  1 
respectively:  Note  that  the  term  V*/2ax2  can  be  expressed  in  terms  of  the 

target  and  clutter  parameters  as,  follows: 


3  =  V4/2ox2 


(A- 16) 


This  follows  because  ox2  is  the  variance- of  the  complex  normal  distribution 
of  Q(k),  and  is  given  by  (A-12). 

Equation  (A-15)  is  extremely  cumbersome  to  evaluate.  However,  useful 
approximations  can  be  obtained  for  the  case  where  $  is  either  very  large  or 
very  small  [A2]. 
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For  the  case  where  0  >1 ,  a  good  approximation  to  (A-15)  is  given  by: 


E[|Rq(k)|J  S  ot/M  (1  +ig-);  8>1 

On  the  other  hand,  for  0_<1,  (A-15)  is  well  approximated  by: 
EljRqOOll  5^  C  /(l+2ai;/C)/I  (H*/2)  i  8<1 


(A-17) 


(A- 18) 


These  approximations  can  be  put  Into  a  more  convenient  form  by  defining  the 
mean  value  of  the  background  pedestal  -  i.e. ,  the  mean  value  of  |Rq(k)|; 
k*C,  k*  |jm-jn|  as  y: 


M  =  Jfc  J (l+2ot/C)/I 


a  =  '4  u  *  at  ^ 

x  M 

Note  that  in  terms  of  y ,  0  ■  -r  —5 — 

H  vz  - 

and  expressing  the  above  approximations  as  a  function  of  y  and  0 . 
Thus  the  approximation  for  large  0  becomes: 


(A-18a) 


E(|Rq(k)|]  2  y  0(1  +  4*>  S>1 


(A-19) 


Correspondingly,  the  approximation  for  small  0  becomes; 

£{|Rq(k)|]  *  y(l+0/2),  0<1  (A-20) 

Note  that  (A-20)  reduces  to  (A-13a)  when  0*0;  -  in  the  absence  of  a  target 
peak. 

Another  possible  approximation  to  the  mean  value  of  |&q(k)|  at  a  target 
peak  location  is  to  approximate  the  mean  value  of  (Rq(k)|  as  the  squarerroot 
of  the  mean  square  value.  The  mean  square  value  of  |Kq(k)|  At  a  target 
peak,  k«|jm-Jn|,  can  be  found  from  eqn.  (A-14)  as  follows: 

|Rq(k)|2  .  Rq(k)Rq*(k) 

-  |Q(k)  | 2  +  V2  (eJA0Q*(k)  +  e^A0Q(k))  +  V4 
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Hence: 


E|Rq(k)|2  «  E|Q(k)|2  +  V*  (A-21) 

The  expected  values  of  the  middle  term  la  zero  because  Q(k)  la  a  zero-mean 
complex  normal  variable  (see  discussion  of  eqn.  (A-10)  on  p.  A-5). 
Furthermore,  using  eqn.  (A— 12) : 

E|Q(k)|2  -  2ox2  -  C2(i+2ot/C)/I  (A-22) 

Finally  combining  eqns.  (A-21)  and  (A-22)  and  using  the  definition  of  V2  In 
(A-6b)  gives: 


El*q(k) I2  “  2«x2  - 
Thus  E{  |gq(k) | ]  - 


,(k)|2j 


/c2(l+2ot/C)/I  +  (ot/M)2 


(A-23) 


(A-24) 


In  terms  of  the  parameters  u  and  0  this  becomes: 

E(  l®q(k)  |  ]  2  p  /  j  (1+8)  (A-25) 

The  3  approximations  represented  by  (A-19) ,  (A-20)  and  (A-24, 25)  are 
plotted  In  Fig.  A-l.  It  is  seen  that  the  overall  approximation  given  by 
(A-24)  and  (A-25)  is  extremely  good,  being  very  close  to  the  approximation 
(A-17,19)  for  values  of  0>1  and  to  the  approximation  (A-18,20)  for  values 
of  0<1.  Therefore  the  approximation  given  by  A-24  Is  used  In  the  main  body 
of  this  report. 
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Fig.  A- 1  Approximations  to  the  mean  value  of  a  target  peak  relative  to 
that  of  the  clutter  background  pedestal  In  a  range,  auto-correlatlon  magni¬ 
tude  profile  as  a  function  of  the  normalized  polnt-scatterer  target-to- 
clutter  ratio.  The  magnitude  of  the  target  peaks  has  a  Ricean 
distribution,  that  of  the  background  pedestal  is  Rayleigh  distributed. 
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This  report  describes  some  of  the  issues  involved  in  using  stepped -frequency  radar  pulse  trains  to 
obtain  high  resolution  target  range  profiles  and  range  profile  autocorrelation  functions.  Coherent  pro¬ 
cessing  of  the  returns  from  such  pulse  trains  yields  the  target  range  profile,  while  noncoherent  proces¬ 
sing  produces  the  autocorrelation  function  of  that  profile.  The  frequency  stability  requirements  of  the 
source  used  to  generate  the  stepped  frequencies  of  the  pulse  train  are  determined  for  both  the  coherent  and 
noncoher  ent  processing  case.  Finally,  the  target -clutter  contrast  provided  by  the  resulting  range  profiles 
and  range  profile  auto-eorrelaticn  functions  is  determined.  __  ,  ,  -  - 
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